Linear and Nonlinear Thermal Transport in Graphene: Molecular Dynamics Simulations

2011 ◽  
Vol 1347 ◽  
Author(s):  
Bo Qiu ◽  
Yan Wang ◽  
Xiulin Ruan
Keyword(s):  
Molecular Dynamics ◽  
Thermal Transport ◽  
Graphene Nanoribbons ◽  
Ballistic Transport ◽  
Thermal Conductance ◽  
Md Simulations ◽  
Spectral Energy ◽  
Suspended Graphene ◽  
Thermal Rectification ◽  
Dynamics Simulations

AbstractIn this work, we perform molecular dynamics (MD) simulations to study the linear thermal transport in suspended graphene and the nonlinear thermal transport phenomena in graphene nanoribbons (GNR). We use spectral energy density analysis to quantitatively address the relative importance of different types of phonon in thermal transport in suspended graphene. Negative differential thermal conductance (NDTC) and thermal rectification in graphene nanoribbons have been studied using nonequilibrium molecular dyanmics simulations. Ballistic transport regime, sufficient temperature nonlinearity and asymmetry are found to be necessary conditions for the onset of these behaviors.

Necessary conditions for thermal rectification and negative differential thermal conductance in graphene nanoribbons

2011 ◽  
Vol 1347 ◽  
Author(s):  
Yan Wang ◽  
Xiulin Ruan
Keyword(s):  
Molecular Dynamics ◽  
Temperature Gradient ◽  
Necessary Conditions ◽  
Graphene Nanoribbons ◽  
Ballistic Transport ◽  
Thermal Conductance ◽  
Nonequilibrium Molecular Dynamics ◽  
Thermal Rectification ◽  
Preferred Direction ◽  
Dynamics Simulations

ABSTRACTWe have studied negative differential thermal conductance (NDTC) and thermal rectification (TR) in graphene nanoribbons (GNRs) using nonequilibrium molecular dynamics simulations. Strong ballistic transport regime and sufficient temperature gradient are found to be necessary conditions for the onset of both NDTC and TR in GNRs, while the latter also requires asymmetry in structure. Preferred direction of heat transport is also discussed for TR.

Thermal Rectification in Graphene and Carbon Nanotube Systems Using Molecular Dynamics Simulations

2011 ◽  
Author(s):  
Ajit K. Vallabhaneni ◽  
Jiuning Hu ◽  
Yong P. Chen ◽  
Xiulin Ruan
Keyword(s):  
Molecular Dynamics ◽  
Carbon Nanotube ◽  
Temperature Difference ◽  
Graphene Nanoribbons ◽  
Md Simulations ◽  
System Size ◽  
Heat Current ◽  
Thermal Rectification ◽  
Mean Temperature ◽  
Rectification Factor

We investigate the thermal rectification phenomena in asymmetric graphene and carbon nanotube systems using molecular dynamics (MD) simulations. The effects of various parameters, including mean temperature, temperature difference, and system size on rectification factor have been studied. In homogenous triangular graphene nanoribbons (T-GNR), the heat current is normally higher from wide to narrow end than that in the opposite direction, resulting in a positive rectification factor. The rectification factor increases further for a double layered T-GNR. It is also found that varying the parameters like mean temperature can result in reverse of the sign of thermal rectification factor. In the case of carbon nanotube (CNT)–silicon system, the heat current is higher when heat flows from CNT to silicon. The thermal rectification factor is almost independent of the diameter of CNT. In both cases, the rectification factor increases with the imposed temperature difference.

Role of Edge Chirality and Isotope Doping in Thermal Transport and Thermal Rectification in Graphene Nanoribbons

2011 ◽  
Author(s):  
Yan Wang ◽  
Xiulin Ruan
Keyword(s):  
Thermal Conductivity ◽  
Thermal Transport ◽  
Graphene Nanoribbons ◽  
Mass Density ◽  
Transport Processes ◽  
Phonon Modes ◽  
Cross Sectional ◽  
Thermal Rectification ◽  
The Difference ◽  
Dynamics Simulations

Thermal transport processes in graphene nanoribbons (GNRs) within and beyond the linear response regime has been studied using classical molecular dynamics simulations. Zigzag-edged GNRs have higher thermal conductivities than armchair-edged ones, and the difference diminishes with increasing width. Analysis on the cross-sectional distribution of heat flux reveals that edge atoms cannot transport thermal energy as efficiently as interior ones. Edge localization of phonon modes reduces thermal transport through edge carbon atoms, especially on armchair edges, which results in a lower thermal conductivity. Isotope (13C) doping can reduce the thermal conductivity of GNRs by 30%–40% by an addition of only ∼20% isotope atoms. The significant reduction in thermal conductivity is partially attributed to phonon localization induced by isotope defects, which is confirmed by phonon mode participation ratio analysis. We also demonstrate that a GNR asymmetric in edge chirality or mass density can generate considerable thermal rectification, which is essential for developing GNR-based thermal management devices.

scholarly journals Atomistic Nodal Approach: a General Molecular Dynamics Method For Computing Local Thermal Conductance at Atomistic Resolution

2021 ◽  
Author(s):  
Mingxuan Jiang ◽  
Juan D. Olarte-Plata ◽  
Fernando Bresme
Keyword(s):  
Molecular Dynamics ◽  
Thermal Transport ◽  
Temperature Drop ◽  
Thermal Conductance ◽  
Fundamental Property ◽  
Lower Surface ◽  
Coordination Numbers ◽  
Interfacial Thermal Conductance ◽  
Non Equilibrium Molecular Dynamics ◽  
Dynamics Simulations

The Interfacial Thermal Conductance (ITC) is a fundamental property of mate- rials and has particular relevance at the nanoscale. The ITC quanti�es the thermal resistance between materials of dierent compositions or between uids in contact with materials. Furthermore, the ITC determines the rate of cooling/heating of the materi- als and the temperature drop across the interface. Here we propose a method to com- pute local ITCs and temperature drops of nanoparticle- uid interfaces. Our approach resolves the ITC at the atomic level using the atomic coordinates of the nanomaterial as nodes to compute local thermal transport properties. We obtain high-resolution descriptions of the interfacial thermal transport by combining the atomistic nodal ap- proach, computational geometry techniques and \computational farming" using Non- Equilibrium Molecular Dynamics simulations. We illustrate our method by analyzing various nanoparticles as a function of their size and geometry, targeting experimentally relevant structures like capped octagonal rods, cuboctahedrons, decahedrons, rhombic dodecahedrons, cubes, icosahedrons, truncated octahedrons, octahedrons and spheres. We show that the ITC of these very dierent geometries can be accurately described in terms of the local coordination number of the atoms in the nanoparticle surface. Nanoparticle geometries with lower surface coordination numbers feature higher ITCs, and the ITC generally increases with decreasing particle size.

Position Sensitivity Study in Molecular Dynamics Simulations of Self-Organized Development of 3D Nanostructures

2017 ◽  
Vol 885 ◽  
pp. 216-221
Author(s):  
Dávid Fülep ◽  
Ibolya Zsoldos ◽  
István László
Keyword(s):  
Molecular Dynamics ◽  
Graphene Nanoribbons ◽  
Md Simulations ◽  
Sensitivity Study ◽  
Basic Difference ◽  
Self Organized ◽  
3D Nanostructures ◽  
Position Sensitivity ◽  
Dynamics Simulations ◽  
Zigzag Type

The sensitivity of defect free fusion of straight carbon nanotubes from graphene nanoribbons to the position of the nanoribbon edge positions has been investigated. A basic difference between the behavior of armchair and zigzag type nanoribbons was observed. When placing armchair type graphene nanoribbons above each other identical, fitting positions are obtained automatically. Zigzag type graphene nanoribbons, however, must not be placed above each other in identical positions. From the viewpoint of defect-free fusion, according to the MD simulations symmetric on nearly symmetric positions of the ribbon edges are favorable.

scholarly journals Molecular Junctions Enhancing Thermal Transport within Graphene Polymer Nanocomposite: A Molecular Dynamics Study

Nanomaterials ◽  
2021 ◽  
Vol 11 (10) ◽  
pp. 2480
Author(s):  
Alessandro Di Pierro ◽  
Bohayra Mortazavi ◽  
Alberto Fina
Keyword(s):  
Molecular Dynamics ◽  
Polymer Matrix ◽  
Thermal Transport ◽  
Strong Dependence ◽  
Thermal Conductance ◽  
Polymer Nanocomposite ◽  
Molecular Junctions ◽  
Tight Contact ◽  
Thermally Conductive ◽  
Dynamics Simulations

Thermal conductivity of polymer-based (nano)composites is typically limited by thermal resistances occurring at the interfaces between the polymer matrix and the conductive particles as well as between particles themselves. In this work, the adoption of molecular junctions between thermally conductive graphene foils is addressed, aiming at the reduction of the thermal boundary resistance and eventually lead to an efficient percolation network within the polymer nanocomposite. This system was computationally investigated at the atomistic scale, using classical Molecular Dynamics, applied the first time to the investigation of heat transfer trough molecular junctions within a realistic environment for a polymer nanocomposite. A series of Molecular Dynamics simulations were conducted to investigate the thermal transport efficiency of molecular junctions in polymer tight contact, to quantify the contribution of molecular junctions when graphene and the molecular junctions are surrounded by polydimethylsiloxane (PDMS) molecules. A strong dependence of the thermal conductance was found in PDMS/graphene model, with best performances obtained with short and conformationally rigid molecular junctions. Furthermore, the adoption of the molecular linkers was found to contribute additionally to the thermal transport provided by the surrounding polymer matrix, demonstrating the possibility of exploiting molecular junctions in composite materials.

Thermal Rectification of Silicene Nanosheets With Triangular Cavities by Molecular Dynamics Simulations

2017 ◽  
Vol 139 (5) ◽  
Author(s):  
Yuan Feng ◽  
Xingang Liang
Keyword(s):  
Molecular Dynamics ◽  
Thermal Transport ◽  
Simulation Method ◽  
Nonequilibrium Molecular Dynamics ◽  
Dimensional Structure ◽  
Honeycomb Lattice ◽  
Rectification Effect ◽  
Thermal Rectification ◽  
Dynamics Simulations ◽  
Silicon Based

Silicene, the silicon-based two-dimensional structure with honeycomb lattice, has been discovered and expected to have tremendous application potential in fundamental industries. However, its thermal transport mechanism and thermal properties of silicene have not been fully explained. We report a possible way to control the thermal transport and thermal rectification in silicene nanosheets by distributing triangular cavities, which are arranged in a staggered way. The nonequilibrium molecular dynamics (NEMD) simulation method is used. The influences of the size, number, and distribution of cavities are investigated. The simulation results show that reflections of phonon at the vertex and the base of the triangular cavities are quite different. The heat flux is higher when heat flow is from the vertex to the base of cavities, resulting in thermal rectification effect. The thermal rectification effect is strengthened with increasing cavity size and number. A maximum of thermal rectification with varying distance between columns of cavities is observed.

scholarly journals Atomistic Nodal Approach: a General Molecular Dynamics Method For Computing Local Thermal Conductance at Atomistic Resolution

2021 ◽  
Author(s):  
Mingxuan Jiang ◽  
Juan D. Olarte-Plata ◽  
Fernando Bresme
Keyword(s):  
Molecular Dynamics ◽  
Thermal Transport ◽  
Temperature Drop ◽  
Thermal Conductance ◽  
Fundamental Property ◽  
Lower Surface ◽  
Fluid Interfaces ◽  
Coordination Numbers ◽  
Non Equilibrium Molecular Dynamics ◽  
Dynamics Simulations

The Interfacial Thermal Conductance (ITC) is a fundamental property of materials and has particular relevance at the nanoscale. The ITC quantifies the thermal resistance between materials of different compositions or between fluids in contact with materials. Furthermore, the ITC determines the rate of cooling/heating of the materials and the temperature drop across the interface. Here we propose a method to compute local ITCs and temperature drops of nanoparticle-fluid interfaces. Our approach resolves the ITC at the atomic level using the atomic coordinates of the nanomaterial as nodes to compute local thermal transport properties. We obtain high-resolution descriptions of the interfacial thermal transport by combining the atomistic nodal approach, computational geometry techniques and "computational farming'' using Non-Equilibrium Molecular Dynamics simulations. We illustrate our method by analyzing various nanoparticles as a function of their size and geometry, targeting experimentally relevant structures like capped octagonal rods, cuboctahedrons, decahedrons, rhombic dodecahedrons, cubes, icosahedrons, truncated octahedrons, octahedrons and spheres. We show that the ITC of these very different geometries can be accurately described in terms of the local coordination number of the atoms in the nanoparticle surface. Nanoparticle geometries with lower surface coordination numbers feature higher ITCs, and the ITC generally increases with decreasing particle size.

Assessing the Antimalarial Potentials of Phytochemicals: Virtual Screening, Molecular Dynamics and In-Vitro Investigations

2019 ◽  
Vol 16 (3) ◽  
pp. 291-300
Author(s):  
Saumya K. Patel ◽  
Mohd Athar ◽  
Prakash C. Jha ◽  
Vijay M. Khedkar ◽  
Yogesh Jasrai ◽  
...  
Keyword(s):  
Molecular Dynamics ◽  
Structural Integrity ◽  
Md Simulations ◽  
Tylophora Indica ◽  
Multidrug Resistance Protein 1 ◽  
Receptor Complexes ◽  
Resistance Protein ◽  
Dynamics Simulations ◽  
Stable Trajectory

Background: Combined in-silico and in-vitro approaches were adopted to investigate the antiplasmodial activity of Catharanthus roseus and Tylophora indica plant extracts as well as their isolated components (vinblastine, vincristine and tylophorine). </P><P> Methods: We employed molecular docking to prioritize phytochemicals from a library of 26 compounds against Plasmodium falciparum multidrug-resistance protein 1 (PfMDR1). Furthermore, Molecular Dynamics (MD) simulations were performed for a duration of 10 ns to estimate the dynamical structural integrity of ligand-receptor complexes. </P><P> Results: The retrieved bioactive compounds viz. tylophorine, vinblastin and vincristine were found to exhibit significant interacting behaviour; as validated by in-vitro studies on chloroquine sensitive (3D7) as well as chloroquine resistant (RKL9) strain. Moreover, they also displayed stable trajectory (RMSD, RMSF) and molecular properties with consistent interaction profile in molecular dynamics simulations. </P><P> Conclusion: We anticipate that the retrieved phytochemicals can serve as the potential hits and presented findings would be helpful for the designing of malarial therapeutics.

scholarly journals Unveiling Carbon Dioxide and Ethanol Diffusion in Carbonated Water-Ethanol Mixtures by Molecular Dynamics Simulations

Molecules ◽  
2021 ◽  
Vol 26 (6) ◽  
pp. 1711
Author(s):  
Mohamed Ahmed Khaireh ◽  
Marie Angot ◽  
Clara Cilindre ◽  
Gérard Liger-Belair ◽  
David A. Bonhommeau
Keyword(s):  
Carbon Dioxide ◽  
Molecular Dynamics ◽  
Diffusion Coefficients ◽  
Hydrodynamic Radius ◽  
Md Simulations ◽  
Molecular Models ◽  
Experimental Values ◽  
Carbon Dioxide Co2 ◽  
Dynamics Simulations ◽  
Apparent Disagreement

The diffusion of carbon dioxide (CO2) and ethanol (EtOH) is a fundamental transport process behind the formation and growth of CO2 bubbles in sparkling beverages and the release of organoleptic compounds at the liquid free surface. In the present study, CO2 and EtOH diffusion coefficients are computed from molecular dynamics (MD) simulations and compared with experimental values derived from the Stokes-Einstein (SE) relation on the basis of viscometry experiments and hydrodynamic radii deduced from former nuclear magnetic resonance (NMR) measurements. These diffusion coefficients steadily increase with temperature and decrease as the concentration of ethanol rises. The agreement between theory and experiment is suitable for CO2. Theoretical EtOH diffusion coefficients tend to overestimate slightly experimental values, although the agreement can be improved by changing the hydrodynamic radius used to evaluate experimental diffusion coefficients. This apparent disagreement should not rely on limitations of the MD simulations nor on the approximations made to evaluate theoretical diffusion coefficients. Improvement of the molecular models, as well as additional NMR measurements on sparkling beverages at several temperatures and ethanol concentrations, would help solve this issue.

Sign in / Sign up

Export Citation Format

Share Document