-
Temperature Dependence of Phonon Energies and Lifetimes in Single- and Few-layered Graphene
Authors:
Markos Poulos,
Konstantinos Papagelis,
Emmenuel N. Koukaras,
George Kalosakas,
Giorgia Fugallo,
Konstantinos Termentzidis
Abstract:
In this work, we have studied the phonon properties of multi-layered graphene with the use of Molecular Dynamics (MD) simulations and the k-space Autocorrelation Sequence (k-VACS) method. We calculate the phonon dispersion curves, densities of states and lifetimes $τ$ of few-layered graphene of 1-5 layers and graphite. $Γ$-point phonon energies and lifetimes are investigated for different temperat…
▽ More
In this work, we have studied the phonon properties of multi-layered graphene with the use of Molecular Dynamics (MD) simulations and the k-space Autocorrelation Sequence (k-VACS) method. We calculate the phonon dispersion curves, densities of states and lifetimes $τ$ of few-layered graphene of 1-5 layers and graphite. $Γ$-point phonon energies and lifetimes are investigated for different temperatures ranging from 80 K to 1000 K. The study focuses on the impact of the interlayer interaction and temperature on the energies and lifetimes of the $Γ$-point phonons, as well as the type of interlayer potential used. For the later we used the Kolmogorov-Crespi (KC) and the Lennard-Jones (LJ) potentials. We have found that the number of layers $N$ has little effect on the intra-layer (ZO and G) mode energies and greater effect on the inter-layer (Layer Shearing and Layer Breathing) modes, while $τ$ is generally affected by $N$ for all modes, except for the Layer Shear mode. The trend of $N$ on the lifetimes was also found to independent of the type of potential used. For the Raman-active G phonon, our calculations show that the lifetime increase with $N$ and that this increase is directly connected to the strength of the interlayer coupling and how this is modelled.
△ Less
Submitted 27 June, 2024;
originally announced June 2024.
-
Bubble relaxation dynamics in homopolymer DNA sequences
Authors:
Malcolm Hillebrand,
George Kalosakas,
Alan R Bishop,
Charalampos Skokos
Abstract:
Understanding the inherent timescales of large bubbles in DNA is critical to a thorough comprehension of its physicochemical characteristics, as well as their potential role on helix opening and biological function. In this work we employ the coarse-grained Peyrard-Bishop-Dauxois model of DNA to study relaxation dynamics of large bubbles in homopolymer DNA, using simulations up to the microsecond…
▽ More
Understanding the inherent timescales of large bubbles in DNA is critical to a thorough comprehension of its physicochemical characteristics, as well as their potential role on helix opening and biological function. In this work we employ the coarse-grained Peyrard-Bishop-Dauxois model of DNA to study relaxation dynamics of large bubbles in homopolymer DNA, using simulations up to the microsecond time scale. By studying energy autocorrelation functions of relatively large bubbles inserted into thermalised DNA molecules, we extract characteristic relaxation times from the equilibration process for both adenine-thymine (AT) and guanine-cytosine (GC) homopolymers. Bubbles of different amplitudes and widths are investigated through extensive statistics and appropriate fittings of their relaxation. Characteristic relaxation times increase with bubble height and width. We show that, within the model, relaxation times are two orders of magnitude longer in GC sequences than in AT sequences. Overall, our results confirm that large bubbles leave a lasting impact on the molecule's dynamics, for times between 0.4-200ns depending on the homopolymer type and bubble shape, thus clearly affecting long-time evolutions of the molecule.
△ Less
Submitted 11 August, 2022;
originally announced August 2022.
-
Distributions of Bubble Lifetimes and Bubble Lengths in DNA
Authors:
Malcolm Hillebrand,
George Kalosakas,
Haris Skokos,
Alan Bishop
Abstract:
We investigate the distribution of bubble lifetimes and bubble lengths in DNA at physiological temperature, by performing extensive molecular dynamics simulations with the Peyrard-Bishop-Dauxois (PBD) model, as well as an extended version (ePBD) having a sequence-dependent stacking interaction, emphasizing the effect of the sequences' guanine-cytosine (GC)/adenine-thymine (AT) content on these dis…
▽ More
We investigate the distribution of bubble lifetimes and bubble lengths in DNA at physiological temperature, by performing extensive molecular dynamics simulations with the Peyrard-Bishop-Dauxois (PBD) model, as well as an extended version (ePBD) having a sequence-dependent stacking interaction, emphasizing the effect of the sequences' guanine-cytosine (GC)/adenine-thymine (AT) content on these distributions. For both models we find that base pair-dependent (GC vs AT) thresholds for considering complementary nucleotides to be separated are able to reproduce the observed dependence of the melting temperature on the GC content of the DNA sequence. Using these thresholds for base pair openings, we obtain bubble lifetime distributions for bubbles of lengths up to ten base pairs as the GC content of the sequences is varied, which are accurately fitted with stretched exponential functions. We find that for both models the average bubble lifetime decreases with increasing either the bubble length or the GC content. In addition, the obtained bubble length distributions are also fitted by appropriate stretched exponential functions and our results show that short bubbles have similar likelihoods for any GC content, but longer ones are substantially more likely to occur in AT-rich sequences. We also show that the ePBD model permits more, longer-lived, bubbles than the PBD system.
△ Less
Submitted 19 November, 2020; v1 submitted 20 August, 2020;
originally announced August 2020.
-
Chaotic dynamics of graphene and graphene nanoribbons
Authors:
Malcolm Hillebrand,
Bertin Many Manda,
George Kalosakas,
Enrico Gerlach,
Charalampos Skokos
Abstract:
We study the chaotic dynamics of graphene structures, considering both a periodic, defect free, graphene sheet and graphene nanoribbons (GNRs) of various widths. By numerically calculating the maximum Lyapunov exponent, we quantify the chaoticity for a spectrum of energies in both systems. We find that for all cases, the chaotic strength increases with the energy density, and that the onset of cha…
▽ More
We study the chaotic dynamics of graphene structures, considering both a periodic, defect free, graphene sheet and graphene nanoribbons (GNRs) of various widths. By numerically calculating the maximum Lyapunov exponent, we quantify the chaoticity for a spectrum of energies in both systems. We find that for all cases, the chaotic strength increases with the energy density, and that the onset of chaos in graphene is slow, becoming evident after more than $10^4$ natural oscillations of the system. For the GNRs, we also investigate the impact of the width and chirality (armchair or zigzag edges) on their chaotic behavior. Our results suggest that due to the free edges the chaoticity of GNRs is stronger than the periodic graphene sheet, and decreases by increasing width, tending asymptotically to the bulk value. In addition, the chaotic strength of armchair GNRs is higher than a zigzag ribbon of the same width. Further, we show that the composition of ${}^{12}C$ and ${}^{13}C$ carbon isotopes in graphene has a minor impact on its chaotic strength.
△ Less
Submitted 12 June, 2020; v1 submitted 20 March, 2020;
originally announced March 2020.
-
Heterogeneity and chaos in the Peyrard-Bishop-Dauxois DNA model
Authors:
Malcolm Hillebrand,
George Kalosakas,
Adrian Schwellnus,
Charalampos Skokos
Abstract:
We discuss the effect of heterogeneity on the chaotic properties of the Peyrard-Bishop-Dauxois nonlinear model of DNA. Results are presented for the maximum Lyapunov exponent and the deviation vector distribution. Different compositions of adenine-thymine (AT) and guanine-cytosine (GC) base pairs are examined for various energies up to the melting point of the corresponding sequence. We also consi…
▽ More
We discuss the effect of heterogeneity on the chaotic properties of the Peyrard-Bishop-Dauxois nonlinear model of DNA. Results are presented for the maximum Lyapunov exponent and the deviation vector distribution. Different compositions of adenine-thymine (AT) and guanine-cytosine (GC) base pairs are examined for various energies up to the melting point of the corresponding sequence. We also consider the effect of the alternation index, which measures the heterogeneity of the DNA chain through the number of alternations between different types (AT or GC) of base pairs, on the chaotic behavior of the system. Biological gene promoter sequences have been also investigated, showing no distinct behavior of the maximum Lyapunov exponent.
△ Less
Submitted 13 February, 2019; v1 submitted 16 November, 2018;
originally announced November 2018.
-
Uniaxial Compression of Suspended Single and Multilayer graphenes
Authors:
A. P. Sgouros,
G. Kalosakas,
C. Galiotis,
K. Papagelis
Abstract:
The mechanical response of single and multiple graphene sheets under uniaxial compressive loads was studied with molecular dynamics simulations, using different semi-empirical force fields at different boundary conditions or constrains. Compressive stress-strain curves were obtained and the critical stress/strain values were derived. For single layer graphenes, the critical stress/strain for buckl…
▽ More
The mechanical response of single and multiple graphene sheets under uniaxial compressive loads was studied with molecular dynamics simulations, using different semi-empirical force fields at different boundary conditions or constrains. Compressive stress-strain curves were obtained and the critical stress/strain values were derived. For single layer graphenes, the critical stress/strain for buckling was found to scale to the inverse length square. For multilayer graphenes the critical buckling stress also decreased with increasing length, though at a slower rate than expected from elastic buckling analysis. The molecular dynamics results are compared to the linear elasticity continuum theory for loaded slabs. Qualitatively similar behavior is observed between the theory and numerical simulations for single layer graphene, while discrepancies were noted for multilayers may be due to their discrete nature and stacking interactions.
△ Less
Submitted 28 April, 2018;
originally announced April 2018.
-
Atomistic potential for graphene and other sp$^2$ carbon systems
Authors:
Zacharias G. Fthenakis,
George Kalosakas,
Georgios D. Chatzidakis,
Costas Galiotis,
Konstantinos Papagelis,
Nektarios N. Lathiotakis
Abstract:
We introduce a torsional force field for sp$^2$ carbon to augment an in-plane atomistic potential of a previous work (Kalosakas et al, J. Appl. Phys. {\bf 113}, 134307 (2013)) so that it is applicable to out-of-plane deformations of graphene and related carbon materials. The introduced force field is fit to reproduce DFT calculation data of appropriately chosen structures. The aim is to create a f…
▽ More
We introduce a torsional force field for sp$^2$ carbon to augment an in-plane atomistic potential of a previous work (Kalosakas et al, J. Appl. Phys. {\bf 113}, 134307 (2013)) so that it is applicable to out-of-plane deformations of graphene and related carbon materials. The introduced force field is fit to reproduce DFT calculation data of appropriately chosen structures. The aim is to create a force field that is as simple as possible so it can be efficient for large scale atomistic simulations of various sp$^2$ carbon structures without significant loss of accuracy. We show that the complete proposed potential reproduces characteristic properties of fullerenes and carbon nanotubes. In addition, it reproduces very accurately the out-of-plane ZA and ZO modes of graphene's phonon dispersion as well as all phonons with frequencies up to 1000~cm$^{-1}$.
△ Less
Submitted 27 November, 2017; v1 submitted 1 November, 2017;
originally announced November 2017.
-
A Torsional potential for graphene derived from fitting to DFT results
Authors:
Georgios D. Chatzidakis,
George Kalosakas,
Zacharias G. Fthenakis,
Nektarios N. Lathiotakis
Abstract:
We present a simple torsional potential for graphene to accurately describe its out-of-plane deformations. The parameters of the potential are derived through appropriate fitting with suitable DFT calculations regarding the deformation energy of graphene sheets folded around two different folding axes, along an armchair or along a zig-zag direction. Removing the energetic contribution of bending a…
▽ More
We present a simple torsional potential for graphene to accurately describe its out-of-plane deformations. The parameters of the potential are derived through appropriate fitting with suitable DFT calculations regarding the deformation energy of graphene sheets folded around two different folding axes, along an armchair or along a zig-zag direction. Removing the energetic contribution of bending angles, using a previously introduced angle bending potential, we isolate the purely torsional deformation energy, which is then fitted to simple torsional force fields. The presented out-of-plane torsional potential can accurately fit the deformation energy for relatively large torsional angles up to 0.5 rad. To test our proposed potential, we apply it to the problem of the vertical displacement of a single carbon atom out of the graphene plane and compare the obtained deformation energy with corresponding DFT calculations. The dependence of the deformation energy on the vertical displacement of the pulled carbon atom is indistinguishable in these two cases, for displacements up to about 0.5 $Å$. The presented potential is applicable to other sp$^2$ carbon structures.
△ Less
Submitted 27 July, 2017;
originally announced July 2017.
-
Graphenes flakes under controlled biaxial deformation
Authors:
Charalampos Androulidakis,
Emmanuel N. Koukaras,
John Parthenios,
George Kalosakas,
Konstantinos Papagelis,
Costas Galiotis
Abstract:
Thin membranes, such as monolayer graphene of monoatomic thickness, are bound to exhibit lateral buckling under uniaxial tensile loading that impairs its mechanical behaviour. In this work, we have developed an experimental device to subject 2D materials to controlled equibiaxial strain on supported beams that can be flexed up or down to subject the material to either compression or tension, respe…
▽ More
Thin membranes, such as monolayer graphene of monoatomic thickness, are bound to exhibit lateral buckling under uniaxial tensile loading that impairs its mechanical behaviour. In this work, we have developed an experimental device to subject 2D materials to controlled equibiaxial strain on supported beams that can be flexed up or down to subject the material to either compression or tension, respectively. Using strain gauges in tandem with Raman spectroscopy measurements, we monitor the G and 2D phonon properties of graphene under biaxial strain and thus extract important information about the uptake of stress under these conditions. The experimental shift over strain for the G and 2D Raman peaks were found to be in the range of 62.3 \pm 5 cm^-1/%, and 148.2 \pm 6 cm^-1/%, respectively, for monolayer but also bilayer graphenes. The corresponding Gruneisen parameters for the G and 2D peaks were found to be between 1.97 \pm 0.15 and 2.86 \pm 0.12, respectively. These values agree reasonable well with those obtained from small-strain bubble-type experiments. The results presented are also backed up by classical and ab initio molecular dynamics simulations and excellent agreement of Γ-E2g shifts with strains and the Gruneisen parameter was observed.
△ Less
Submitted 13 June, 2016;
originally announced June 2016.
-
Raman spectroscopy of graphene at high pressure: substrate and pressure transmitting media effects
Authors:
K. Filintoglou,
N. Papadopoulos,
J. Arvanitidis,
D. Christofilos,
O. Frank,
M. Kalbac,
J. Parthenios,
G. Kalosakas,
C. Galiotis,
K. Papagelis
Abstract:
The pressure evolution of the Raman spectrum of graphene grown by chemical vapour deposition on polycrystalline copper is investigated with the use of a polar and a non-polar pressure transmitting medium (PTM). The G and 2D Raman bands exhibit similar pressure slopes for both PTM irrespectively of any unintentional initial doping and/or strain of the samples. Our analysis suggests that any pressur…
▽ More
The pressure evolution of the Raman spectrum of graphene grown by chemical vapour deposition on polycrystalline copper is investigated with the use of a polar and a non-polar pressure transmitting medium (PTM). The G and 2D Raman bands exhibit similar pressure slopes for both PTM irrespectively of any unintentional initial doping and/or strain of the samples. Our analysis suggests that any pressure-induced charge transfer effects are negligibly small to influence the pressure response of graphene. This is determined by the mechanical stress due to the pressure-induced substrate contraction and the transfer efficiency of the latter to the graphene layer, as well as the PTM-graphene interaction. For the non-polar PTM, a peculiar pressure behavior of graphene is observed in the PTM solidification regime, resembling that of free standing graphene.
△ Less
Submitted 30 April, 2013;
originally announced April 2013.
-
In-plane force fields and elastic properties of graphene
Authors:
G. Kalosakas,
N. N. Lathiotakis,
C. Galiotis,
K. Papagelis
Abstract:
Bond stretching and angle bending force fields, appropriate to describe in-plane motion of graphene sheets, are derived using first principles' methods. The obtained force fields are fitted by analytical anharmonic energy potential functions, providing efficient means of calculations in molecular mechanics simulations. Numerical results regarding the mechanical behavior of graphene monolayers unde…
▽ More
Bond stretching and angle bending force fields, appropriate to describe in-plane motion of graphene sheets, are derived using first principles' methods. The obtained force fields are fitted by analytical anharmonic energy potential functions, providing efficient means of calculations in molecular mechanics simulations. Numerical results regarding the mechanical behavior of graphene monolayers under various loads, like uniaxial tension, hydrostatic tension, and shear stress, are presented, using both molecular dynamics simulations and first principles' methods. Stress-strain curves and elastic constants, such as, Young modulus, Poisson ratio, bulk modulus, and shear modulus, are calculated. Our results are compared with corresponding theoretical calculations as well as with available experimental estimates. Finally, the effect of the anharmonicity of the extracted potentials on the mechanical properties of graphene are discussed.
△ Less
Submitted 21 November, 2012;
originally announced November 2012.
-
Dependence on temperature and GC content of bubble length distributions in DNA
Authors:
G. Kalosakas,
S. Ares
Abstract:
We present numerical results on the temperature dependence of the distribution of bubble lengths in DNA segments of various guanine-cytosine (GC) concentrations. Base-pair openings are described by the Peyrard-Bishop-Dauxois model and the corresponding thermal equilibrium distributions of bubbles are obtained through Monte Carlo calculations for bubble sizes up to the order of a hundred base pai…
▽ More
We present numerical results on the temperature dependence of the distribution of bubble lengths in DNA segments of various guanine-cytosine (GC) concentrations. Base-pair openings are described by the Peyrard-Bishop-Dauxois model and the corresponding thermal equilibrium distributions of bubbles are obtained through Monte Carlo calculations for bubble sizes up to the order of a hundred base pairs. The dependence of the parameters of bubble length distribution on temperature and the GC content is investigated. We provide simple expressions which approximately describe these relations. The variation of the average bubble length is also presented. We find a temperature dependence of the exponent c that appears in the distribution of bubble lengths. If an analogous dependence exists in the loop entropy exponent of real DNA, it may be relevant to understand overstretching in force-extension experiments.
△ Less
Submitted 25 June, 2009;
originally announced June 2009.
-
Probability distribution of bubble lengths in DNA
Authors:
S. Ares,
G. Kalosakas
Abstract:
The distribution of bubble lengths in double-stranded DNA is presented for segments of varying guanine-cytosine (GC) content, obtained with Monte Carlo simulations using the Peyrard-Bishop-Dauxois model at 310 K. An analytical description of the obtained distribution in the whole regime investigated, i.e., up to bubble widths of the order of tens of nanometers, is available. We find that the dec…
▽ More
The distribution of bubble lengths in double-stranded DNA is presented for segments of varying guanine-cytosine (GC) content, obtained with Monte Carlo simulations using the Peyrard-Bishop-Dauxois model at 310 K. An analytical description of the obtained distribution in the whole regime investigated, i.e., up to bubble widths of the order of tens of nanometers, is available. We find that the decay lengths and characteristic exponents of this distribution show two distinct regimes as a function of GC content. The observed distribution is attributed to the anharmonic interactions within base pairs. The results are discussed in the framework of the Poland-Scheraga and the Peyrard-Bishop (with linear instead of nonlinear stacking interaction) models.
△ Less
Submitted 16 February, 2007; v1 submitted 5 July, 2006;
originally announced July 2006.
-
Nonlocal interactions in doped cuprates: correlated motion of Zhang-Rice polarons
Authors:
L. Hozoi,
S. Nishimoto,
G. Kalosakas,
D. B. Bodea,
S. Burdin
Abstract:
In-plane, inter-carrier correlations in hole doped cuprates are investigated by ab initio multiconfiguration calculations. The dressed carriers display features that are reminiscent of both Zhang-Rice (ZR) CuO4 singlet states and Jahn-Teller polarons. The interaction between these quasiparticles is repulsive. At doping levels that are high enough, the interplay between long-range unscreened Coul…
▽ More
In-plane, inter-carrier correlations in hole doped cuprates are investigated by ab initio multiconfiguration calculations. The dressed carriers display features that are reminiscent of both Zhang-Rice (ZR) CuO4 singlet states and Jahn-Teller polarons. The interaction between these quasiparticles is repulsive. At doping levels that are high enough, the interplay between long-range unscreened Coulomb interactions and long-range phase coherence among the O-ion half-breathing vibrations on the ZR plaquettes may lead to a strong reduction of the effective adiabatic energy barrier associated to each polaronic state. Tunneling effects cannot be neglected for a relatively flat, multi-well energy landscape. We suggest that the coherent, superconducting quantum state is the result of such coherent quantum lattice fluctuations involving the in-plane O ions. Our findings appear to support models where the superconductivity is related to a lowering of the in-plane kinetic energy.
△ Less
Submitted 23 May, 2006;
originally announced May 2006.
-
Comment on "Can one predict DNA Transcription Start Sites by Studying Bubbles?"
Authors:
C. H. Choi,
A. Usheva,
G. Kalosakas,
K. O. Rasmussen,
A. R. Bishop
Abstract:
Comment on T.S. van Erp, S. Cuesta-Lopez, J.-G. Hagmann, and M. Peyrard, Phys. Rev. Lett. 95, 218104 (2005) [arXiv: physics/0508094].
Comment on T.S. van Erp, S. Cuesta-Lopez, J.-G. Hagmann, and M. Peyrard, Phys. Rev. Lett. 95, 218104 (2005) [arXiv: physics/0508094].
△ Less
Submitted 9 December, 2005;
originally announced December 2005.
-
Structurally specific thermal fluctuations identify functional sites for DNA transcription
Authors:
G. Kalosakas,
K. Ø. Rasmussen,
A. R. Bishop,
C. H. Choi,
A. Usheva
Abstract:
We report results showing that thermally-induced openings of double stranded DNA coincide with the location of functionally relevant sites for transcription. Investigating both viral and bacterial DNA gene promoter segments, we found that the most probable opening occurs at the transcription start site. Minor openings appear to be related to other regulatory sites. Our results suggest that coher…
▽ More
We report results showing that thermally-induced openings of double stranded DNA coincide with the location of functionally relevant sites for transcription. Investigating both viral and bacterial DNA gene promoter segments, we found that the most probable opening occurs at the transcription start site. Minor openings appear to be related to other regulatory sites. Our results suggest that coherent thermal fluctuations play an important role in the initiation of transcription. Essential elements of the dynamics, in addition to sequence specificity, are nonlinearity and entropy, provided by local base-pair constraints.
△ Less
Submitted 8 September, 2003; v1 submitted 5 September, 2003;
originally announced September 2003.