1. Introduction
Helical structures formed at the molecular level have attracted the attention of researchers for a long time. Many molecular structures in living organisms exist in the form of helices, in particular DNA and RNA molecules and various peptides. Such molecules are natural polymers. The successive development of methods of polymer organic chemistry has led to the possibility of the synthesis of artificial stereoregular polymers, the threads of which are also ordered in a spiral manner (helical polymers)
[1][2][3]. The spiral topology opens up great opportunities for various technological applications regarding helical polymers. In particular, the circularly polarized luminescence (CPL) of cis-cisoid polyene-based polymers
[4], the magneto-optical properties of helical poly-3-(alkylsulfone)thiophene
[5], various stimuli-responsive devices as sensors
[6], asymmetric catalysis and chiral recognition
[2], and chirality-induced spin selectivity (CISS)
[7] should be mentioned.
The development of experimental methods for the synthesis of nanostructures that are much more extended in one dimension than in the other two (for brevity, researchers call such objects one-dimensional, 1D) made it possible to obtain nano objects of various morphologies and structures, including helical ones
[8]. Among the most frequently used synthetic methods that allow for 1D helical objects to be obtained are the following: (i) the screw-dislocation-driven growth of 1D crystals
[8] and (ii) various methods based on the addition of a chiral component to the reaction mixture (template-assisted synthesis)
[9]. The 1D nanostructures with intrinsic chirality can be noted separately. The existence of such 1D crystals is due to the fact that the symmetry of the initial 3D crystals of these compounds is described by the Sohncke symmetry groups
[10]. Among such structures, recently synthesized ultrathin tellurium nanorods can be mentioned
[11].
Helical 1D nanostructures, here referred to as
nanohelicenes [12], with both helical polymers and nanorods with intrinsic chirality can be considered. In molecular [
n]helicenes, the number
n denotes the number of ortho-fused benzene rings
[13]. If
n tends to infinity, then a polymeric [∞]helicene will be obtained, which can be considered as a ribbon infinitely turned along the helical axis (
Figure 1). By considering [∞]helicene as the basic nanohelicene, it is possible to construct many modifications that will have different properties. At the same time, all nanohelicenes will have two main properties—a continuous π-system, and a helical topology.
Figure 1. Molecular [n]helicenes and nanohelicene. From left to right: [6]helicene, [7]helicene, [16]helicene, [26]helicene, and [∞]helicene = nanohelicene. Gray and white spheres are carbon and hydrogen atoms, respectively.
Further, the consideration of nanohelicenes as “graphite” nanorods with a screw dislocation
[14][15] makes it possible to understand these objects as structures formed by a graphene ribbon of a special morphology, which allows it to helically turn around a screw axis.
Therefore, other names are also used for nanohelicenes: “graphene spiral”
[16], “graphene helicoid”
[17], “helical graphene nanoribbon”
[18], and “helical graphene”
[19]. In addition, conventional graphene nanoribbons can be twisted into “helicoidal graphene nanoribbons”
[20], which are qualitatively different from the objects considered here.
2. Theoretical Consideration of Molecular Helicenes and Their Properties
Molecular helicenes have a long story. In 1956, the [6]helicene was synthesized
[21], which can be considered as one coil of the [∞]helicene. The synthesis of helicenes with two coils began in 1967 with the synthesis of the [7]helicene
[22]. In 2015, the [14]helicene was obtained
[23], i.e., helicene with three incomplete coils, and, to date, helicene with maximum
n (
Figure 1) remains to be obtained experimentally. It can be assumed that a further increase in
n is a difficult synthetic problem, which is due to significant steric hindrance existing for inner atoms. However, there are some experimental data
[24] that enable one to speculate on the presence of helicenes in interstellar space up to the [26]helicene (with more than four coils), which can be explained on the basis of the principles of minimum energy and polarizability (
Figure 1)
[25][26][27].
Thus, at the moment, the properties of helicenes with a further increase in
n, and the approximation of properties to
n → ∞ can only be estimated based on quantum chemical calculations. The first quantum chemical calculations (at the level of semi-empirical CNDO/S, PM3, and the HF method) of the helicenes appeared quite a long time ago and dealt with the structure of helicenes, racemization barriers, and chiroptical properties
[28][29][30]. In
[31], the effect of the distortion of the π-system due to twisting was theoretically investigated by comparing the energy characteristics of [
n]helicenes and their planar counterparts, [
n]phenacenes. The calculations were carried out at the level of HF/6-31G* and B3LYP/6-31G* at
n = 6–10 and
n = 16 and in the more powerful 6-311G** basis set for
n = 6–10.
The dependence of the difference, Δ
E, between the energies of helicenes and phenacenes with increasing
n was considered
[31]:
The results show a linear increase in the ΔE between helicenes and phenacenes with increasing n. In other words, an increase in n does not lead to any stabilization of helicene. At n = 16 (largest experimentally synthesized helicene), the difference reaches 83.4 kcal/mol (the result of B3LYP/6-31G* calculations). This makes it possible to understand the observed difficulties in obtaining helicenes with a further increase in n.
The difference in magnetic susceptibilities between helicenes and phenacenes was also investigated in
[31]. This quantity enables us to determine the aromaticity of the studied molecules. According to the data of
[31], helical distortions in the π-system expectedly lead to a decrease in the aromaticity of helicenes in comparison with phenacenes. However, this decrease is not so large—the loss of aromaticity in helicenes is approximately 10–13% compared to phenacenes.
Similar results were obtained by the authors of
[32], who performed DFT calculations using the functional B3LYP with the TZVP basis set and the empirical Grimme correction to account for the dispersion forces
[33] and applying the conductor-like screening model (COSMO)
[34] to take into account the solvent (acetonitrile). Incremental Gibbs free energies, Δ
Gn, were obtained, which show the difference between the homodesmic addition of naphthalene (C
10H
8) to helicene or phenacene:
Here,
G([
n]) and
G([
n − 1]) are the Gibbs free energies of [
n]helicene or [
n]phenacene and [
n − 1]helicene or [
n − 1]phenacene, respectively. The Δ
Gn reaches a constant value of ~10 kJ/mol at
n ~ 6 for helicenes, which qualitatively corresponds to the result of
[31] on inherent instability and explains the difficulty of the synthesis of helicenes with a large
n. A similar result was obtained at the M06-2X/6-311+G(3d,p) level of calculations
[35].
Researchers also note that the small difference between the aromaticity of helicenes and phenacenes noted earlier in
[31] is confirmed by the results of subsequent works
[32][35].
For the first time, the possibility of using helicene in terms of springs for application in nanotechnologies was discussed in
[36]. From this point of view, at the level of the semi-empirical PM3 Hamiltonian, the force constants of axial tension for neutral helicenes and their cationic and anionic forms were estimated. Calculations were made for helicenes with
n = 12 and 18, which corresponds to the number of coils of two and three, respectively. According to
[36], the force constants are in the range 9.8–14.6 kcal·mol
−1·Å
−2 for small “springs” and approximately in the range 4.9–5.4 kcal·mol
−1·Å
−2 for large “springs”. The differences in hardness depend little on the electronic state of helicenes (i.e., neutral, cationic, or anionic states).
In
[32], a somewhat larger value was obtained for [14]helicene, equal to 23.6 kcal·mol
−1·Å
−2. It can be assumed that the difference is due to a more accurate calculation method, which was carried out within the DFT (B3LYP functional) in the TZVP basis set, taking into account the solvent and dispersion corrections. A further theoretical study of the influence of the axial deformations of [∞]helicene and its non-hexagonal modifications on mechanical stability, deformability, and fracture processes was carried out in
[37] at the PBE/PAW level of theory. It was found that the presence of pentagons or heptagons instead of hexagons had a significant influence on the mechanical properties, and some helicenes could achieve very large reversible tensile strains (more than 200%). The roughly estimated effective tensile strength of pristine [∞]helicene is equal to 117 GPa.
The authors of
[32] also made an attempt to computationally evaluate the convergence of the physicochemical properties of [
n]helicenes with increasing
n. The results of TD DFT calculations of excitation energies
[32] agree well with experiments and demonstrate the convergence of the results at
n ~ 14. In particular, for [14]helicene, the calculated wavelength for the
S0 →
S1 vertical transition is 476 nm (2.60 eV).
The convergence of the gap between the HOMO (highest occupied molecular orbital) and LUMO (lowest unoccupied molecular orbital) was studied at the B3LYP/6-31G* computational level
[38]. Calculations were carried out for helicenes and phenacenes (and their analogs with sulfur) with an
n of up to 30 and under PBCs (periodic boundary conditions) in order to obtain the change in the electronic band gap of the structures with
n → ∞. The data obtained for [14]helicene in
[38] correlate with the results in
[32], taking into account the fact that the HOMO–LUMO gap is only an upper estimate of the energies of the
S0 →
S1 vertical transition obtained by the TD DFT.
The obtained extrapolations are very close to the results of the calculations under PBCs (2.90 eV and 3.59 eV for [∞]helicene and [∞]phenacene, respectively
[38]). The value 2.90 eV correlates with the 2.49 eV value obtained for [∞]helicene by the authors of earlier work
[39].
Ref.
[39] is the first work in which the infinite [∞]helicene was also considered (in the periodic model at the tight-binding level). In addition, researchers note that in
[39], not only simple [∞]helicene but also its other topological variants (anthra-helicene and benzo-helicene) were considered for the first time. [∞]benzo-helicene has been shown to be potentially metallic, with a zero band gap, unlike simple [∞]helicene, which is a semiconductor.
The band structure of these two 1D structures shows that the [∞]benzo-helicene conduction band is “inside” the band gap of [∞]helicene, which allowed the authors of Ref.
[39] to propose a variant of a quantum dot constructed as a hybrid of semi-conductor helicene and metallic benzo-helicene. The
I/
V characteristics (i.e., current-voltage) were calculated for finite size
[21]benzo-helicene, [24]helicene, and their hybrid. As expected from the band structure, the
I/
V curve of the hybrid is similar to that of semiconductor helicene. This can also be understood in terms of the molecular orbitals (MO) of the hybrid. It can be seen from the figure that the MOs, which correspond to the energy levels within the band gap of the semiconductor part, will be localized on the metal part of the hybrid. If the voltage exceeds the band gap value of the semiconductor part, then the conductivity is restored, and MOs are delocalized throughout the hybrid. Hybrids that combine different types of helicenes have not been theoretically studied anywhere except in
[39].
Researchers also note that in
[39], for the first time, such issues as the Peierls distortion of metallic helicene were raised, and a decrease in the band gap with “expanding” helicene (i.e., a transition from simple helicene to anthra-helicene) was also shown. These issues were repeatedly analyzed in subsequent theoretical works on nanohelicenes.
In
[40], the stability and properties of some other classes of helicenes were theoretically studied at the B3LYP/6-31G* level. Some ways for constructing a spiral from various combinations of hexagons, pentagons, and tetragons from carbon atoms were considered, including the “lateral-extending” of the helicene’s π-system, which is the object of study in subsequent theoretical works.
The presence of a continuous π-system and axial chirality makes helicenes and their derivatives promising for achieving non-linear optical (NLO) properties. At the level of the semi-empirical AM1 Hamiltonian, within the framework of the time-dependent Hartree-Fock (TD HF) theory, the change in the first hyperpolarizability (
) and its projection on the dipole moment (
) were theoretically studied with a successive increase in
n (up to
n = 19) for helicenes and their antiaromatic analogs, phenylenes
[41]. It was found that
increases monotonically with the size of the system. At the same time, in helicenes,
is positive and presents quasi-periodic oscillations with the helix.
It was theoretically shown at the PBEPBE/6-31G(d,p) level that the modification of the helicene backbone by fusing the azulene element leads to a significant enhancement of the second-order NLO properties due to the polar charge distribution in the pentagon end (a component of azulene)
[42]. NLO properties have also been shown to exhibit only slight changes under an external strain.
Further works in this area concerned third-order NLO properties, two-photon circular dichroism, and so on. A detailed review of the calculations and experiments on the NLO properties of helicenes and their derivatives is presented in a recent book
[43].
Theoretical studies of electron transport through helicenes were the subject of works
[44][45][46][47][48]. The common result for all the works mentioned is the conclusion about the presence of a strong dependence of conductivity and other properties on the applied mechanical stress. This is a fairly expected result based on the understanding of helicenes as molecular springs.
In particular, helicene-based molecular junctions constructed as (diaza)helicenes placed between two Au(111) electrodes were studied in
[44]. Helicene compression leads to a significant increase in conductance, whereas stretching leads to a decrease. In the range of low voltages, conductance falls off linearly on a logarithmic scale as a function of strain. The thermopower shows a sign change, which implies a change in the transport mechanism of the junction. The dependence of the thermoelectric
figure of merit (FOM),
zT, on strain for increasing
n values shows that the maximal values of FOM were obtained for the longest helicene. The maximal
zT values can be as large as ~0.6, which far exceeds the values for currently known molecular junctions.
In
[45], a wider tension-compression interval for [12]helicene was theoretically investigated, and the results on
I/V behavior correlate with the results of
[44] in the corresponding interval. However, the use of a wider interval made it possible to detect an increase in conductivity in the region of large stretches. Thus, current dependence on pitch has a
U-shape. The first part of the curve (the distance between coils is less than 4 Å) corresponds to the transmission pathways between helicene coils (through space). The ascending part of the curve corresponds to a change in the conduction mechanism, and the conduction is due to transmission along the helicene ribbon (through bond).
At PBE/DZP computational level, the authors of
[46] investigated the electronic transport properties of spiral-shaped molecules coupled to metal electrodes. Three molecules were composed of [7]helicene and two thiol-[7]helicenes with benzenethiolates on both sides. The data on local transmission pathways demonstrate that there is no spiral current flowing through [7]helicene, and the total current is only due to electron tunneling in the inter-layer. In the case of thiol-[7]helicenes, the electron distribution in frontier orbitals will change, and the spiral current will appear, which makes it possible to generate a magnetic field. The negative differential resistance behavior of
I/
V curves was also found in
[46]. In other words, in a certain interval of bias voltage, the current decreases with bias increases.
A similar result about the absence of a helical current for helicenes was obtained in
[47]. It was also shown in
[48] that the conductivity of the oligomers of [
n]benzo-helicene (see
[39]) is significantly higher than that of conventional [
n]helicenes. In the limit, this corresponds to the metallic conductivity of [∞]benzo-helicene that was shown in
[39].
In addition to electrical conductivity studies, the theoretical calculations were also concerned with the potential application of helicene molecules to control spin transport. Such control has many advantages over conventional electronic transport and can potentially be used in nanodevices that are smaller and more robust than modern nanodevices. The use of helicene derivatives as a spin filter exhibiting the CISS effect was shown experimentally
[49]. According to the general theory, the driving force behind CISS is spin-orbit interaction, SOC. Despite the fact that for carbon atoms, the SOC value is very small (about 6 meV), the CISS for these molecules can be significant
[50], which opens up great opportunities for the use of organic chiral molecules, and in particular, helicenes.
The first theoretical consideration of spin-polarized electron transport in helicene was carried out in
[51]. The [
n]helicenes (
n = 12, 16, 20, and 24) combined with semi-infinite graphene nanoribbons were considered. The spin polarization,
Ps, is determined through spin-up
G↑ and spin-down
G↓ conductance, according to Equation (3):
Using the tight-binding Hamiltonian with the addition of the SOC Hamiltonian, it was shown that the
Ps of [20]helicene could reach 27%, which is comparable to the experimentally found values in
[49]. Interestingly,
Ps exhibits significant oscillations by varying the Fermi energy,
E, because of quantum interference effects. In addition, the influence of axial deformations on
Ps was studied. Axial deformations significantly change the overall shape of the
Ps(
E) curve; however, the oscillation behavior remains in the energy spectrum regardless of the applied mechanical stress, and the spin filtration efficiency remains quite large. The
Ps(
E) curves for other helicenes differ significantly in terms of peak heights and their location, but the overall picture remains the same as in the case of [20]helicene.
Since the CISS effect relates primarily to excited electron states, it should be observed not only in the case of flowing current but also in the case of a photocurrent promoted by photo-excited electrons. In
[52], an experimental and theoretical study of [7]helicene molecules adsorbed on metal surfaces (Cu(332), Ag(110), and Au(111)) was carried out. An excess of 6–8% of longitudinal spin polarization,
PZ, was demonstrated:
Here,
I↑,↓ denotes the count rates in the two detectors, and
Seff is the effective Sherman function. It is important to note that circularly polarized UV/Vis light served as the source of spin polarization. Thus, the obtained experimental results connect two promising phenomena—CISS and CPL, which is possible in this case. In the computational part of Ref.
[52], the DFT method was used to study the [7]helicene molecules adsorbed on different surfaces. Then, using the model Hamiltonian to calculate the CISS,
PZ was obtained for helicenes with different numbers of coils. The results showed a significant increase in
PZ for the system with the maximum number of coils (three in the case under study). Interestingly, the authors of
[53] explored the CISS associated with inter-system crossing (ISC) in carbon-sulfur [
n]helicenes. It has been suggested that the efficient ISC of these molecules might be beneficial for an efficient CISS.
Spin transport through helicenes bonded to a carbon zigzag-terminated nanoribbon was theoretically studied in
[54]. Some peripheral modifications of helicene with one coil were considered, starting from the basic [6]helicene. The helicene periphery was successively expanded by adding from one to three carbon hexagons. The calculations were carried out at the DFT+NEGF (non-equilibrium Green function) level with the DZP basis set.
It was shown that in the case of basic helicene, both spin states on the central spiral moiety reveal transport via inter-layer tunneling without the generation of spiral current, which is similar to the results of
[45][46]. Spin-polarized transport is present but is relatively small. However, modification by adding carbon hexagons resulted in a robust spin-polarized interaction, especially in the case of the appearance of trigonal sections. According to the Lieb theorem
[55] (regarding the total spin of the Hubbard model in bipartite sublattices
A and
B), these parts have extra atoms that break the balance between the number of atoms in sublattices, i.e.,
NA ≠
NB. Consequently, the emergence of magnetism is induced by edge imbalance. Additionally, depending on the location of additional hexagons, helicene-based devices exhibit an on-off spin-polarization phenomenon, which manifests itself in the presence/absence of spin-splitting that is dependent on bias. Interestingly, the location of additional hexagons also affects the type of electron tunneling. In the case of trigonal modification, it was possible to redirect the current completely in a spiral, i.e., intra-layer tunneling is a major contribution to current, which is in contrast to the basic helicene case.
Despite the fact that the theoretical works considered mainly refer to molecular helicenes, the manifestation of many of the effects found can also be expected for nanohelicenes. Researchers especially note the importance of structural modification, which consists of the addition of carbon hexagons to the carbon skeleton of the initial [n]helicene. From the results of theoretical studies, it can be clearly seen that modification leads to a significant change in the optical, electronic, and magnetic properties, and taking into account the spring topology of helicenes makes such modifications extremely promising for creating corresponding opto-, electro-, and magneto-mechanical devices at the nano level.