1
|
Mahalingam V, Harursampath D. Effect of the Base-Pair Sequence on B-DNA Thermal Conductivity. J Phys Chem B 2021; 125:10652-10656. [PMID: 34533965 DOI: 10.1021/acs.jpcb.1c04318] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The thermal conductivity of double-stranded (ds) B-DNA was systematically investigated using classical molecular dynamics (MD) simulations. The effect of changing the base pair (bp) on the thermal conductivity of dsDNA needed investigation at a molecular level. Hence, four sequences, viz., poly(A), poly(G), poly(CG), and poly(AT), were initially analyzed in this work. First, the length of these sequences was varied from 4 to 40 bp at 300 K, and the respective thermal conductivity (κ) was computed. Second, the temperature-dependent thermal conductivities between 100 and 400 K were obtained in 50 K steps at 28 bp length. The Müller-Plathe reverse nonequilibrium molecular dynamics (RNEMD) was employed to set a thermal gradient and obtain all thermal conductivities in this work. Moreover, mixed sequences using AT and CG sequences, namely, A(CG)nT (n = 3-7), ACGC(AT)mGCGT (m = 0-5), and ACGC(AT)nAGCGT (n = 1-4), were investigated based on the hypothesis that these sequences could be better thermoelectrics. One-dimensional lattices are said to have diverging thermal conductivities at longer lengths, which violate the Fourier law. These follow the power law, where κ ∝ Lβ. At longer lengths, the exponent β needs to satisfy the condition β > 1/3 for divergent thermal conductivity. We find no such significant Fourier law violation through divergence of thermal conductivities at 80 bp lengths or 40 bp lengths. Also, in the case of the second study, the presence of short (m ≤ 2) encapsulated AT sequences within CG sequences shows an increasing trend. These results are important for engineering DNA-based thermal devices.
Collapse
Affiliation(s)
- Vignesh Mahalingam
- Department of Aerospace Engineering, Indian Institute of Science, Bengaluru 560012, India
| | | |
Collapse
|
2
|
Elenewski JE, Wójtowicz G, Rams MM, Zwolak M. Performance of reservoir discretizations in quantum transport simulations. J Chem Phys 2021; 155:124117. [PMID: 34598565 DOI: 10.1063/5.0065799] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Abstract
Quantum transport simulations often use explicit, yet finite, electronic reservoirs. These should converge to the correct continuum limit, albeit with a trade-off between discretization and computational cost. Here, we study this interplay for extended reservoir simulations, where relaxation maintains a bias or temperature drop across the system. Our analysis begins in the non-interacting limit, where we parameterize different discretizations to compare them on an even footing. For many-body systems, we develop a method to estimate the relaxation that best approximates the continuum by controlling virtual transitions in Kramers turnover for the current. While some discretizations are more efficient for calculating currents, there is little benefit with regard to the overall state of the system. Any gains become marginal for many-body, tensor network simulations, where the relative performance of discretizations varies when sweeping other numerical controls. These results indicate that typical reservoir discretizations have little impact on numerical costs for certain computational tools. The choice of a relaxation parameter is nonetheless crucial, and the method we develop provides a reliable estimate of the optimal relaxation for finite reservoirs.
Collapse
Affiliation(s)
- Justin E Elenewski
- Biophysical and Biomedical Measurement Group, Microsystems and Nanotechnology Division, Physical Measurement Laboratory, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, USA
| | - Gabriela Wójtowicz
- Jagiellonian University, Institute of Theoretical Physics, Łojasiewicza 11, 30-348 Kraków, Poland
| | - Marek M Rams
- Jagiellonian University, Institute of Theoretical Physics, Łojasiewicza 11, 30-348 Kraków, Poland
| | - Michael Zwolak
- Biophysical and Biomedical Measurement Group, Microsystems and Nanotechnology Division, Physical Measurement Laboratory, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, USA
| |
Collapse
|
3
|
Abstract
The thermal conductivity of B-form double-stranded DNA (dsDNA) of the Drew-Dickerson sequence d(CGCGAATTCGCG) is computed using classical molecular dynamics (MD) simulations. In contrast to previous studies, which focus on a simplified 1D model or a coarse-grained model of DNA to reduce simulation times, full atomistic simulations are employed to understand the thermal conduction in B-DNA. Thermal conductivities at different temperatures from 100 to 400 K are investigated using the Einstein-Green-Kubo equilibrium and Müller-Plathe non-equilibrium formalisms. The thermal conductivity of B-DNA at room temperature is found to be 1.5 W/m·K in equilibrium and 1.225 W/m·K in the non-equilibrium approach. In addition, the denaturation regime of B-DNA is obtained from the variation of thermal conductivity with temperature. It is in agreement with previous studies using the Peyrard-Bishop-Dauxois model at a temperature of around 350 K. The quantum heat capacity (Cvq) has given additional clues regarding the Debye and denaturation temperature of 12-bp B-DNA.
Collapse
Affiliation(s)
- Vignesh Mahalingam
- Department of Aerospace Engineering, Indian Institute of Science, Bengaluru 560012, India
| | | |
Collapse
|
4
|
Zwolak M. Analytic expressions for the steady-state current with finite extended reservoirs. J Chem Phys 2020; 153:224107. [PMID: 33317280 PMCID: PMC8356363 DOI: 10.1063/5.0029223] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Open-system simulations of quantum transport provide a platform for the study of true steady states, Floquet states, and the role of temperature, time dynamics, and fluctuations, among other physical processes. They are rapidly gaining traction, especially techniques that revolve around "extended reservoirs," a collection of a finite number of degrees of freedom with relaxation that maintains a bias or temperature gradient, and have appeared under various guises (e.g., the extended or mesoscopic reservoir, auxiliary master equation, and driven Liouville-von Neumann approaches). Yet, there are still a number of open questions regarding the behavior and convergence of these techniques. Here, we derive general analytical solutions, and associated asymptotic analyses, for the steady-state current driven by finite reservoirs with proportional coupling to the system/junction. In doing so, we present a simplified and unified derivation of the non-interacting and many-body steady-state currents through arbitrary junctions, including outside of proportional coupling. We conjecture that the analytic solution for proportional coupling is the most general of its form for isomodal relaxation (i.e., relaxing proportional coupling will remove the ability to find compact, general analytical expressions for finite reservoirs). These results should be of broad utility in diagnosing the behavior and implementation of extended reservoir and related approaches, including the convergence to the Landauer limit (for non-interacting systems) and the Meir-Wingreen formula (for many-body systems).
Collapse
Affiliation(s)
- Michael Zwolak
- Biophysical and Biomedical Measurement Group, Microsystems and Nanotechnology Division, Physical Measurement Laboratory, National Institute of Standards and Technology, Gaithersburg, MD 20899, USA
| |
Collapse
|
5
|
Wójtowicz G, Elenewski JE, Rams MM, Zwolak M. Open System Tensor Networks and Kramers' Crossover for Quantum Transport. PHYSICAL REVIEW. A 2020; 101:10.1103/PhysRevA.101.050301. [PMID: 33367191 PMCID: PMC7754794 DOI: 10.1103/physreva.101.050301] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Tensor networks are a powerful tool for many-body ground states with limited entanglement. These methods can nonetheless fail for certain time-dependent processes-such as quantum transport or quenches-where entanglement growth is linear in time. Matrix-product-state decompositions of the resulting out-of-equilibrium states require a bond dimension that grows exponentially, imposing a hard limit on simulation timescales. However, in the case of transport, if the reservoir modes of a closed system are arranged according to their scattering structure, the entanglement growth can be made logarithmic. Here, we apply this ansatz to open systems via extended reservoirs that have explicit relaxation. This enables transport calculations that can access steady states, time dynamics and noise, and periodic driving (e.g., Floquet states). We demonstrate the approach by calculating the transport characteristics of an open, interacting system. These results open a path to scalable and numerically systematic many-body transport calculations with tensor networks.
Collapse
Affiliation(s)
- Gabriela Wójtowicz
- Jagiellonian University, Institute of Theoretical Physics, Lojasiewicza 11, 30-348 Kraków, Poland
| | - Justin E. Elenewski
- Biophysics Group, Microsystems and Nanotechnology Division, Physical Measurement Laboratory, National Institute of Standards and Technology, Gaithersburg, MD, USA
- Institute for Research in Electronics and Applied Physics, University of Maryland, College Park, MD, USA
| | - Marek M. Rams
- Jagiellonian University, Institute of Theoretical Physics, Lojasiewicza 11, 30-348 Kraków, Poland
| | - Michael Zwolak
- Biophysics Group, Microsystems and Nanotechnology Division, Physical Measurement Laboratory, National Institute of Standards and Technology, Gaithersburg, MD, USA
| |
Collapse
|
6
|
Topology, landscapes, and biomolecular energy transport. Nat Commun 2019; 10:4662. [PMID: 31604949 PMCID: PMC6789131 DOI: 10.1038/s41467-019-12700-w] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2019] [Accepted: 09/12/2019] [Indexed: 12/20/2022] Open
Abstract
While ubiquitous, energy redistribution remains a poorly understood facet of the nonequilibrium thermodynamics of biomolecules. At the molecular level, finite-size effects, pronounced nonlinearities, and ballistic processes produce behavior that diverges from the macroscale. Here, we show that transient thermal transport reflects macromolecular energy landscape architecture through the topological characteristics of molecular contacts and the nonlinear processes that mediate dynamics. While the former determines transport pathways via pairwise interactions, the latter reflects frustration within the landscape for local conformational rearrangements. Unlike transport through small-molecule systems, such as alkanes, nonlinearity dominates over coherent processes at even quite short time- and length-scales. Our exhaustive all-atom simulations and novel local-in-time and space analysis, applicable to both theory and experiment, permit dissection of energy migration in biomolecules. The approach demonstrates that vibrational energy transport can probe otherwise inaccessible aspects of macromolecular dynamics and interactions that underly biological function. Understanding vibrational energy transfer in macromolecules has been challenging to both theory and experiment. Here the authors use non-equilibrium molecular dynamics to reveal the relationship between heat transport in a model peptide, emergent nonlinearity, and the underlying free energy landscape.
Collapse
|
7
|
Dugar P, Chien CC. Geometry-induced local thermal current from cold to hot in a classical harmonic system. Phys Rev E 2019; 99:022131. [PMID: 30934240 DOI: 10.1103/physreve.99.022131] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2018] [Indexed: 06/09/2023]
Abstract
The second law of thermodynamics requires the overall thermal current to flow from hot to cold. However, it does not forbid a local thermal current from flowing from cold to hot. By coupling a harmonic system of three masses connected by a few springs to two Langevin reservoirs at different temperatures, a local atypical thermal current is found to flow from cold to hot in the steady state while the overall thermal current is still from hot to cold. The direction of the local thermal current can be tuned by the mass, spring constant, and system-reservoir coupling. The local thermal current can vanish if the parameters are tuned to proper values. We also consider nonlinear effect from the system-substrate coupling and find that the local atypical thermal current survives in the presence of the nonlinear potential. Moreover, the local atypical thermal current is robust against asymmetry of the system-reservoir coupling, inhomogeneity of the nonlinear potential, and additions of more masses and springs. In molecular or nanomechanical systems where the setup may find its realization, the direction of the local thermal current may be controlled by mechanical or electromagnetic means, which may lead to applications in information storage.
Collapse
Affiliation(s)
- Palak Dugar
- School of Natural Sciences, University of California, Merced, California 95343, USA
| | - Chih-Chun Chien
- School of Natural Sciences, University of California, Merced, California 95343, USA
| |
Collapse
|
8
|
Marty-Roda M, Dahlen O, van Erp TS, Cuesta-López S. Improving the mesoscopic modeling of DNA denaturation dynamics. Phys Biol 2018; 15:066001. [PMID: 29775183 DOI: 10.1088/1478-3975/aac61c] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
Although previously developed mesoscopic DNA models have successfully reproduced thermodynamic denaturation data, recent studies show that these overestimate the rate of denaturation by orders of magnitude. Using adapted Peyrard-Bishop-Dauxois (PBD) models, we have calculated the denaturation rates of several DNA hairpins and made comparison with experimental data. We show that the addition of a barrier at the onsite potential of the PBD model gives a more accurate description of the unzipping dynamics of short DNA sequences. The new models provide a refined theoretical insight on the dynamical mechanisms of unzipping which can have implications for the understanding of transcription and replication processes. Still, this class of adapted PBD models seems to have a fundamental limitation which implies that it is not possible to get agreement with available experimental results on the dynamics of long DNA sequences and at the same time maintain the good agreement regarding its thermodynamics. The reason for this is that the denaturation rate of long DNA chains is not dramatically lowered by the additional barrier-as the base-pairs that open are more likely to remain open, facilitating the opening of the full DNA molecule. Some care has to be taken, since experimental techniques suitable to the study of denaturation rates of long sequences seem not to agree with other experimental data on short DNA sequences. Further research, both theoretical and experimental, is therefore needed to resolve these inconsistencies-which will be a starting point for new minimalistic models that are able to describe both thermodynamics and dynamics at a predictive level.
Collapse
Affiliation(s)
- Marta Marty-Roda
- International Research Center in Critical Raw Materials-ICCRAM, University of Burgos, Plaza Misael Bañuelos s/n, 09001 Burgos, Spain. Advanced Materials, Nuclear Technology and Applied Bio/Nanotechnology, Consolidated Research Unit UIC-154, Castilla y Leon, Spain
| | | | | | | |
Collapse
|
9
|
|
10
|
Chien CC, Velizhanin KA, Dubi Y, Ilic BR, Zwolak M. Topological quantization of energy transport in micro- and nano-mechanical lattices. PHYSICAL REVIEW. B 2018; 97:10.1103/PhysRevB.97.125425. [PMID: 30997441 PMCID: PMC6463522 DOI: 10.1103/physrevb.97.125425] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/09/2023]
Abstract
Topological effects typically discussed in the context of quantum physics are emerging as one of the central paradigms of physics. Here, we demonstrate the role of topology in energy transport through dimerized micro- and nano-mechanical lattices in the classical regime, i.e., essentially "masses and springs". We show that the thermal conductance factorizes into topological and nontopological components. The former takes on three discrete values and arises due to the appearance of edge modes that prevent good contact between the heat reservoirs and the bulk, giving a length-independent reduction of the conductance. In essence, energy input at the boundary mostly stays there, an effect robust against disorder and nonlinearity. These results bridge two seemingly disconnected disciplines of physics, namely topology and thermal transport, and suggest ways to engineer thermal contacts, opening a direction to explore the ramifications of topological properties on nanoscale technology.
Collapse
Affiliation(s)
- Chih-Chun Chien
- School of Natural Sciences, University of California, Merced, CA 95343, USA
| | | | - Yonatan Dubi
- Department of Chemistry and the Ilse Katz Institute for Nanoscale Science and Technology, Ben-Gurion University of the Negev, Beer-Sheva 84105, Israel
| | - B. Robert Ilic
- Center for Nanoscale Science and Technology, National Institute of Standards and Technology, Gaithersburg, MD 20899, USA
| | - Michael Zwolak
- Center for Nanoscale Science and Technology, National Institute of Standards and Technology, Gaithersburg, MD 20899, USA
| |
Collapse
|
11
|
Chien CC, Kouachi S, Velizhanin KA, Dubi Y, Zwolak M. Thermal transport in dimerized harmonic lattices: Exact solution, crossover behavior, and extended reservoirs. Phys Rev E 2017; 95:012137. [PMID: 28208409 PMCID: PMC5473347 DOI: 10.1103/physreve.95.012137] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2016] [Indexed: 11/07/2022]
Abstract
We present a method for calculating analytically the thermal conductance of a classical harmonic lattice with both alternating masses and nearest-neighbor couplings when placed between individual Langevin reservoirs at different temperatures. The method utilizes recent advances in analytic diagonalization techniques for certain classes of tridiagonal matrices. It recovers the results from a previous method that was applicable for alternating on-site parameters only, and extends the applicability to realistic systems in which masses and couplings alternate simultaneously. With this analytic result in hand, we show that the thermal conductance is highly sensitive to the modulation of the couplings. This is due to the existence of topologically induced edge modes at the lattice-reservoir interface and is also a reflection of the symmetries of the lattice. We make a connection to a recent work that demonstrates thermal transport is analogous to chemical reaction rates in solution given by Kramers' theory [Velizhanin et al., Sci. Rep. 5, 17506 (2015)]2045-232210.1038/srep17506. In particular, we show that the turnover behavior in the presence of edge modes prevents calculations based on single-site reservoirs from coming close to the natural-or intrinsic-conductance of the lattice. Obtaining the correct value of the intrinsic conductance through simulation of even a small lattice where ballistic effects are important requires quite large extended reservoir regions. Our results thus offer a route for both the design and proper simulation of thermal conductance of nanoscale devices.
Collapse
Affiliation(s)
- Chih-Chun Chien
- School of Natural Sciences, University of California, Merced, California 95343, USA
| | - Said Kouachi
- Department of Mathematics, Qassim University, Al-Gassim, Buraydah 51452, Saudi Arabia
| | - Kirill A Velizhanin
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - Yonatan Dubi
- Department of Chemistry and the Ilse Katz Institute for Nanoscale Science and Technology, Ben-Gurion University of the Negev, Beer-Sheva 84105, Israel
| | - Michael Zwolak
- Center for Nanoscale Science and Technology, National Institute of Standards and Technology, Gaithersburg, Maryland 20899, USA
| |
Collapse
|
12
|
Velizhanin KA, Sahu S, Chien CC, Dubi Y, Zwolak M. Crossover behavior of the thermal conductance and Kramers' transition rate theory. Sci Rep 2015; 5:17506. [PMID: 26634333 PMCID: PMC4669443 DOI: 10.1038/srep17506] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2015] [Accepted: 10/23/2015] [Indexed: 11/26/2022] Open
Abstract
Kramers' theory frames chemical reaction rates in solution as reactants overcoming a barrier in the presence of friction and noise. For weak coupling to the solution, the reaction rate is limited by the rate at which the solution can restore equilibrium after a subset of reactants have surmounted the barrier to become products. For strong coupling, there are always sufficiently energetic reactants. However, the solution returns many of the intermediate states back to the reactants before the product fully forms. Here, we demonstrate that the thermal conductance displays an analogous physical response to the friction and noise that drive the heat current through a material or structure. A crossover behavior emerges where the thermal reservoirs dominate the conductance at the extremes and only in the intermediate region are the intrinsic properties of the lattice manifest. Not only does this shed new light on Kramers' classic turnover problem, this result is significant for the design of devices for thermal management and other applications, as well as the proper simulation of transport at the nanoscale.
Collapse
Affiliation(s)
| | - Subin Sahu
- Center for Nanoscale Science and Technology, National Institute of Standards and Technology, Gaithersburg, MD 20899, USA
- Maryland Nanocenter, University of Maryland, Collage Park, MD 20742, USA
- Department of Physics, Oregon State University, Corvallis, OR 97331, USA
| | - Chih-Chun Chien
- School of Natural Sciences, University of California, Merced, CA 95343, USA
| | - Yonatan Dubi
- Department of Chemistry and the Ilse Katz Institute for Nanoscale Science and Technology, Ben-Gurion University of the Negev, Beer-Sheva 84105, Israel
| | - Michael Zwolak
- Center for Nanoscale Science and Technology, National Institute of Standards and Technology, Gaithersburg, MD 20899, USA
- Department of Physics, Oregon State University, Corvallis, OR 97331, USA
| |
Collapse
|
13
|
Maximiano RV, Weber G. Deoxyinosine mismatch parameters calculated with a mesoscopic model result in uniform hydrogen bonding and strongly variable stacking interactions. Chem Phys Lett 2015. [DOI: 10.1016/j.cplett.2015.04.045] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
14
|
Li B, Han W, Byun M, Zhu L, Zou Q, Lin Z. Macroscopic highly aligned DNA nanowires created by controlled evaporative self-assembly. ACS NANO 2013; 7:4326-33. [PMID: 23550794 DOI: 10.1021/nn400840y] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/21/2023]
Abstract
By subjecting DNA aqueous solution to evaporate in a curve-on-flat geometry that was composed of either a spherical lens or a cylindrical lens situated on a flat substrate, a set of highly aligned DNA nanowires in the forms of spokes and parallel stripes over a macroscopic area (i.e., millimeter scale) were successfully created. The DNA molecules were stretched and aligned on polymer-coated substrate by the receding meniscus. The imposed curve-on-flat geometry provided a unique environment for controlling the flow within the evaporating solution by eliminating temperature gradient and possible convective instability and, thus, regulated the formation of DNA nanowires. Such controlled evaporative self-assembly is remarkably easy to implement and opens up a new avenue for crafting large-scale DNA-based nanostructures in a simple and cost-effective manner, dispensing with the need for lithography techniques.
Collapse
Affiliation(s)
- Bo Li
- School of Materials Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, United States
| | | | | | | | | | | |
Collapse
|