A method that utilises atomic trajectories and velocities from molecular dynamics simulations has been suitably adapted and employed for the implicit calculation of the phonon dispersion curves of graphene. Classical potentials widely used in the literature were employed. Their performance was assessed for each individual phonon branch and the overall phonon dispersion, using available inelastic x-ray scattering data. The method is promising for systems with large scale periodicity, accounts for anharmonic effects and non-bonding interactions with a general environment and it is applicable under finite temperatures. The temperature dependence of the phonon dispersion curves has been examined with emphasis on the doubly degenerate Raman active Γ-E2g phonon at the zone centre, where experimental results are available. The potentials used show diverse behaviour. The Tersoff-2010 potential exhibits the most systematic and physically sound behaviour in this regard and gives a first-order temperature coefficient of χ = −0.05 cm−1/K for the Γ-E2g shift in agreement with reported experimental values.
It has long been known that double-stranded DNA is subject to temporary, localized openings of its two strands. Particular regions along a DNA polymer are destabilized structurally by available thermal energy in the system. The localized sequence of DNA determines the physical properties of a stretch of DNA, and that in turn determines the opening profile of that DNA fragment. We show that the Peyrard-Bishop nonlinear dynamical model of DNA, which has been used to simulate denaturation of short DNA fragments, gives an accurate representation of the instability profile of a defined sequence of DNA, as verified using S1 nuclease cleavage assays. By comparing results for a non-promoter DNA fragment, the adenovirus major late promoter, the adeno-associated viral P5 promoter and a known P5 mutant promoter that is inactive for transcription, we show that the predicted openings correlate almost exactly with the promoter transcriptional start sites and major regulatory sites. Physicists have speculated that localized melting of DNA might play a role in gene transcription and other processes. Our data link sequence-dependent opening behavior in DNA to transcriptional activity for the first time.
We investigate polaron properties in the semiclassical Holstein model in one, two, and three dimensions, using two methods: a simple and efficient numerical scheme and a variational approach. We obtain accurate information on the energy and the existence regimes of the polaron state. We study small oscillations of the polaron through normal-mode analysis, investigate their stability, and evaluate the density of states of the polaron eigenfrequencies for different system parameters. The normal-mode analysis shows in the one-dimensional case a pinning-breathing mode crossing in the transition region of the small to large polaron. Similar crossings do not exist in higher dimensions.
Bond stretching and angle bending force fields, appropriate to describe in-plane properties of graphene sheets, are derived using first principles' methods. The obtained force fields are fitted by analytical anharmonic potential energy functions, providing efficient means of calculations in molecular mechanics simulations. Using both molecular dynamics simulations and first principles' methods, numerical results regarding the mechanical behavior of graphene monolayers under various loads, like uniaxial tension in different directions or hydrostatic tension, are presented and compared. Graphene's response in shear stress is also investigated using molecular dynamics, where a noticeable asymmetric mechanical behavior is found. Stress-strain curves and elastic constants, such as, Young modulus, Poisson's ratio, bulk modulus, and shear modulus, are calculated. Our results are compared with available experimental estimates, as well as, with corresponding theoretical calculations. Finally, the effects of the anharmonicity of the extracted bond stretching and angle bending potentials on the mechanical properties of graphene are discussed.
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 ± 5 cm–1/% and 148.2 ± 6 cm–1/%, respectively, for monolayer but also bilayer graphenes. The corresponding Grüneisen parameters for the G and 2D peaks were found to be between 1.97 ± 0.15 and 2.86 ± 0.12, respectively. These values agree reasonably 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 Grüneisen parameter was observed.