Density Functional Molecular Dynamics Simulations Research Articles

Mechanical and electronic behaviour of TMDC nanotubes and monolayers: molecular simulations

ABSTRACT Transition metal dichalcogenides (TMDCs) are among the most exciting materials in nano-size particles. In this study, MX2 (M = Mo and W; X = S, Se, and Te) nanotubes and monolayers were evaluated to show the structural, mechanical, and electronic properties. The results showed that the stability of MX2 monolayers changes similarly to nanotubes with the sequence of MS2 > MSe2 > MTe2 (M = Mo and W). Besides, MX2 nanotubes and monolayers are semiconductors, and the band gap energy increases with a similar sequence with stability properties. The structural stability and band gap energy of WX2 structures are more than MoX2 structures. Young's modulus of nanotubes was studied with density functional theory (DFT) and molecular dynamics (MD) simulations with Stillinger-Weber (SW) potential methods. The results are similar to the formation and band gap energy showing that Young's modulus of the MX2 structures is in order of MS2 > MSe2 > MTe2 (M = Mo and W) and WX2 > MoX2 (X = S, Se or Te). Structural defects in monolayers, which were M-vacancy (VM), X2-vacancy (VX2), M-X2-vacancy (VM-X2), and M-X4-vacancy (VM-X4), were studied and it showed that the mechanical properties could be decreased. However, the defects do not affect electronic properties.

Biodegradability assessment of HDPE-based biocomposites: Influence of starch and fiber composition

This research aims to analyze the biodegradation dynamics of a tertiary composite blend, including High-Density Polyethylene (HDPE), starch and linen fiber, and their combined effect on decay processes in authentic environmental settings. It investigates the relationship between fiber content and decomposition rates, details the biodegradation mechanisms, and evaluates the reactive profiles of the involved constituents. Decay kinetics and the biodegradation mechanism of three formulations: HDPE60S40, HDPE60S20F20, and HDPE60S30F10, representing composites with 60% HDPE, complemented by 40%, 20% starch and 20%, 10% linen fiber, respectively, are examined. HDPE60S30F10 is noted for its superior biodegradation rates, showing a 1.2% weight loss in soil and 9.89% in marine conditions and an increased resistance to shearing forces, whereas HDPE60S40 recorded a weight loss of 0,63% in soil and 2.59% in seawater against 1,7% and 6.64% in soil and seawtaer, respectively recorded with HDPE60S20F20. Density Functional Theory (DFT) and Molecular Dynamics (MD) simulations complement these findings, presenting HDPE60S40 as the most rigid, HDPE60S20F20 as the most ductile with a bulk modulus of 13.34GPa, and HDPE60S30F10 exhibiting the best shear resistance with a shear modulus of 12.48GPa. Scanning Electron Microscopy (SEM) and Fourier-Transform Infrared Spectroscopy (FTIR) analyses confirm microbial involvement and significant surface erosion, particularly indicating particularly starch degradation. The results suggest that integrating linen fiber into the composites enhances biodegradation.

Biodegradability assessment of HDPE-based biocomposites: Influence of starch and fiber composition

This research aims to analyze the biodegradation dynamics of a tertiary composite blend, including High-Density Polyethylene (HDPE), starch and linen fiber, and their combined effect on decay processes in authentic environmental settings. It investigates the relationship between fiber content and decomposition rates, details the biodegradation mechanisms, and evaluates the reactive profiles of the involved constituents. Decay kinetics and the biodegradation mechanism of three formulations: HDPE60S40, HDPE60S20F20, and HDPE60S30F10, representing composites with 60% HDPE, complemented by 40%, 20% starch and 20%, 10% linen fiber, respectively, are examined. HDPE60S30F10 is noted for its superior biodegradation rates, showing a 1.2% weight loss in soil and 9.89% in marine conditions and an increased resistance to shearing forces, whereas HDPE60S40 recorded a weight loss of 0,63% in soil and 2.59% in seawater against 1,7% and 6.64% in soil and seawtaer, respectively recorded with HDPE60S20F20. Density Functional Theory (DFT) and Molecular Dynamics (MD) simulations complement these findings, presenting HDPE60S40 as the most rigid, HDPE60S20F20 as the most ductile with a bulk modulus of 13.34GPa, and HDPE60S30F10 exhibiting the best shear resistance with a shear modulus of 12.48GPa. Scanning Electron Microscopy (SEM) and Fourier-Transform Infrared Spectroscopy (FTIR) analyses confirm microbial involvement and significant surface erosion, particularly indicating particularly starch degradation. The results suggest that integrating linen fiber into the composites enhances biodegradation.

Zea mays bracts extract as an eco-friendly corrosion inhibitor for steel in HCl pickling solution: Experimental and simulation studies

Active phytochemicals with unique electron structures can adsorb on metal surfaces and effectively mitigate metal corrosion. Recently, plant extracts have garnered significant attention as green and renewable corrosion inhibitors. In this study, the abundant agricultural byproducts, namely Zea mays bracts, were used as the extract raw materials, and Zea mays bracts extract (ZMBE) was served as an eco-friendly corrosion inhibitor for mild steel (MS) in 1 M HCl pickling solution. The organic substances in ZMBE were analyzed via ultra-high performance liquid chromatography- quadrupole time-of-flight mass spectrometry (UPLC-QTFMS). The interaction between the active functional groups and MS was investigated using fourier transform infrared (FT-IR) and ultraviolet–visible (UV–visible) spectroscopy and X-ray photoelectron spectroscopy (XPS). The corrosion inhibition behavior was scrutinized through potentiodynamic polarization (PDP) curves and electrochemical impedance spectroscopy (EIS). The results show that functional groups such as aromatic rings, double bonds, and hydroxyls successfully adsorb on the MS surface to form the protective films. ZMBE exhibits a mixed-type corrosion inhibition with a “geometric coverage” effect. The corrosion inhibition efficiency (η) increases with ZMBE concentration and solution temperature, the optimum η with 5.0 g·L−1 ZMBE at 318 K is 96.2 %. The adsorption of active ingredients follows the Langmuir and El-Awady isotherms. Density functional theory (DFT) and molecular dynamics (MD) simulations prove the synergistic effect of complex components in ZMBE and the corrosion inhibition mechanism. ZMBE is a promising acid pickling corrosion inhibitor with broad application prospects.

Coordination structures and stabilities of Am(III) adsorption complexes at kaolinite(0 0 1)-water interface

Understanding the structures and stabilities of actinide species is crucial for comprehending the sorption mechanisms at mineral–water interfaces. However, the relevant molecular-level fundamental aspects of actinide coordination chemistry at the mineral–water interface remain unclear. We herein conduct a systematic theoretical study of coordination structures and adsorption energies of Am3+ and Am(OH)2+ sorption complexes at the kaolinite(001)-water interface by integrating static density functional theory (DFT) calculations and ab initio molecular dynamics (AIMD) simulations. Computational results show that stable outer- and inner-sphere adsorption complexes are respectively formed through hydrogen bonding and coordinate bonding with surface aluminol groups. The most stable inner-sphere species are predicted to be tridentate surface complexes with an eight-folded coordination structure in the first coordination sphere of Am(III). This work deepens the understanding of Am(III) adsorption and complexation behaviors at the mineral–water interface, providing essential insights for actinide migration in natural environments and nuclear waste disposal. The synergy between static DFT calculations for initial screening and AIMD simulations for comprehensive dynamic analysis in this work enhances the efficiency and accuracy of theoretical predictions for actinide environmental chemistry and can be helpful for the broader research fields of geochemistry and environmental science.

Unveiling a New 2D Semiconductor: Biphenylene-Based InN.

The two-dimensional (2D) materials class earned a boost in 2021 with biphenylene synthesis, which is structurally formed by the fusion of four-, six-, and eight-membered carbon rings, usually named 4-6-8-biphenylene network (BPN). This research proposes a detailed study of electronic, structural, dynamic, and mechanical properties to demonstrate the potential of the novel biphenylene-like indium nitride (BPN-InN) via density functional theory and molecular dynamics simulations. The BPN-InN has a direct band gap energy transition of 2.02 eV, making it promising for optoelectronic applications. This structure exhibits maximum and minimum Young modulus of 22.716 and 22.063 N/m, Poisson ratio of 0.018 and -0.008, and Shear modulus of 11.448 and 10.860 N/m, respectively. To understand the BPN-InN behavior when subjected to mechanical deformations, biaxial and uniaxial strains in armchair and zigzag directions from -8 to 8% were applied, achieving a band gap energy modulation of 1.36 eV over tensile deformations. Our findings are expected to motivate both theorists and experimentalists to study and obtain these new 2D inorganic materials that exhibit promising semiconductor properties.

Probing the interaction mechanism of choline chloride-based deep eutectic solvents with flavonoid and polyphenol systems

Choline chloride-based deep eutectic solvents (DES) have been demonstrated as promising green solvents for the extraction of flavonoids and polyphenols. This study investigates the interaction between six different hydrogen bond donors (HBD) choline chloride-based DES and three types of polyphenolic flavonoid substrates. Density functional theory (DFT) and molecular dynamics (MD) simulations suggest that Cl− dominates the interaction primarily through electrostatic forces. The extraction effect of alcohol DES on rutin was optimized obviously, and the maximum interaction could reach 737.6 kJ/mol. Acid DES has a better extraction effect on anthocyanins, with a maximum interaction of 376.4 kJ/mol. With the increase of the carbon chain of HBA and the increase of the number of functional groups, the interaction between DES and the active substance decreases. As HBD, ethanolamine with both hydroxyl and amino groups showed the best effect in the extraction process of DES, forming the maximum interaction energy of 1575.8 kJ/mol. Additionally, multivariate analysis is employed to establish the relationship between different variables and interaction energy. This study provides micro-theoretical support for the extraction of flavonoids and polyphenols.

X-functionalized molybdenene monolayers (X = O, F, Al, Si, Cl)

Molybdenene synthesized experimentally is the sole Dirac material with metallic properties so far [Sahu et al., Nat. Nanotechnol. 18, 1430 (2023)], exhibiting a wide array of unique and outstanding properties with potentially extensive applications. However, the free-standing molybdenene structure is not so stable, which highly limits its further exploitation. In this work, we employ density functional theory calculations and ab initio molecular dynamics simulations to investigate molybdenene monolayers functionalized with various elements (X = H, Li, Be, B, C, N, O, F, Na, Mg, Al, Si, P, S, Cl) at both basal surfaces, aiming at achieving the stable free-standing molybdenene derivatives. By evaluating the energetic, mechanical, dynamical, and thermodynamic properties, we confirm some stable monolayer structures of X-functionalized molybdenene (X = O, F, Al, Si, Cl), which are named as X-Molybdenene. Therefore, our findings stabilize the molybdenene via surface functionalization, which is crucial for future experimental validations and applications.

Formation of Supernarrow Borophene Nanoribbons

AbstractBorophenes have sparked considerable interest owing to their fascinating physical characteristics and diverse polymorphism. However, borophene nanoribbons (BNRs) with widths less than 2 nm have not been achieved. Herein, we report the experimental realization of supernarrow BNRs. Combining scanning tunneling microscopy imaging with density functional theory modeling and ab initio molecular dynamics simulations, we demonstrate that, under the applied growth conditions, boron atoms can penetrate the outermost layer of Au(111) and form BNRs composed of a pair of zigzag (2,2) boron rows. The BNRs have a width self‐contained to ∼1 nm and dipoles at the edges to keep them separated. They are embedded in the outermost Au layer and shielded on top by the evacuated Au atoms, free of the need for post‐passivation. Scanning tunneling spectroscopy reveals distinct edge states, primarily attributed to the localized spin at the BNRs’ zigzag edges. This work adds a new member to the boron material family and introduces a new physical feature to borophenes.

Formation of Supernarrow Borophene Nanoribbons.

Borophenes have sparked considerable interest owing to their fascinating physical characteristics and diverse polymorphism. However, borophene nanoribbons (BNRs) with widths less than 2 nm have not been achieved. Herein, we report the experimental realization of supernarrow BNRs. Combining scanning tunneling microscopy imaging with density functional theory modeling and ab initio molecular dynamics simulations, we demonstrate that, under the applied growth conditions, boron atoms can penetrate the outermost layer of Au(111) and form BNRs composed of a pair of zigzag (2,2) boron rows. The BNRs have a width self-contained to ∼1 nm and dipoles at the edges to keep them separated. They are embedded in the outermost Au layer and shielded on top by the evacuated Au atoms, free of the need for post-passivation. Scanning tunneling spectroscopy reveals distinct edge states, primarily attributed to the localized spin at the BNRs' zigzag edges. This work adds a new member to the boron material family and introduces a new physical feature to borophenes.

Exploring the mechanism of nerve agent (Tabun and Sarin) adsorption on carbon nanocones: Computational insights

Carbon NanoCones (CNCs), distinguished by their conical shape and large surface area, are evaluated for their potential application in nerve gas [tabun (GA) and sarin (GB)] adsorption. Density Functional Theory (DFT) calculations, Molecular Dynamics (MD), and Monte Carlo (MC) simulations are employed to investigate this phenomenon comprehensively. The study findings indicate that tabun (GA) and sarin (GB) undergo spontaneous adsorption, characterized by significant negative energy values. Density Functional Theory (DFT) calculations indicate electron donation tendencies of nerve agents to CNC surfaces. Computational analysis shows higher adsorption efficiency in the inner portion (P1) of CNCs. Molecular Dynamics (MD) and Monte Carlo (MC) simulations confirm planar adsorption configurations parallel to CNC surfaces. The DFT calculations show that GA has a far greater adsorption energy on a Carbon Nanocone compared to GB (−44.96 kcal/mol for GA, whereas GB has an adsorption energy of −35.95 kcal/mol). Furthermore, the analysis reveals significant differences between Sarin and Tabun nerve agents. Sarin exhibits higher global hardness, indicating increased stability, while Tabun shows higher overall softness, suggesting greater reactivity. Despite similar electronegativities, Sarin’s greater ionization energy and electron affinity imply enhanced stability and reactivity. Differences in HOMO and LUMO energies highlight unique reactivity profiles, with Sarin being more susceptible to nucleophilic attacks and Tabun to electrophilic attack. Sarin’s larger energy gap between HOMO and LUMO orbitals signifies its superior stability under electronic disturbances.

Computational Insights into Chromene/pyran Derivatives: Molecular Docking, ADMET Studies, DFT Calculations, and MD Simulations as Promising Candidates for Parkinson's Disease.

Parkinson's disease (PD) is a neurodegenerative condition characterized by both motor and non-motor symptoms. Although PD is commonly associated with a decline of dopaminergic neurons in the substantia nigra, other diagnostic criteria and biomarkers also exist. In the search for novel therapeutic agents, chromene and pyran derivatives have shown potential due to their diverse pharmacological activities. This study utilizes a comprehensive computational approach to investigate the viability of chromene/pyran compounds as potential treatments for PD. The drug-likeness characteristics of these molecules were analyzed using ADMET (Absorption, Distribution, Metabolism, Excretion, and Toxicity) studies. Molecular docking was performed against PDB ID: 2V5Z. The best three molecules chosen were compound 7, compound 24, and compound 67 have a binding energy of -6.7, -8.6, and -10.9 kcal/mol. Molecules demonstrating positive blood-brain barrier permeability, good solubility, and favorable binding affinity were further evaluated using Density Functional Theory (DFT) calculations and Molecular Dynamics (MD) simulations to assess their electronic structure and stability. DFT calculations indicated that molecule 82 has a dipole moment of 15.70 D. RMSD and RMSF results confirmed the stability of the complexes over a 100 ns simulation, with a maximum of 3 hydrogen bonds formed.

The cage effect of electron beam irradiation damage in cryo-electron microscopy

Electron beam irradiation can cause damage to biological and organic samples, as determined via transmission electron microscopy (TEM). Cryo-electron microscopy (cryo-EM) significantly reduces such damage by quickly freezing the environmental water around organic molecules. However, there are multiple hypotheses about the mechanism of cryo-protection in cryo-EM. A lower temperature can cause less molecular dissociation in the first stage, or frozen water can have a “cage” effect by preventing the dissociated fragments from flying away. In this work, we use real-time time-dependent density functional theory molecular dynamics(rt-TDDFT-MD) simulations to study the related dynamics. We use our recently developed natural orbital branching (NOB) algorithm to describe the molecular dissociation process after the molecule is ionized. We find that despite the difference in surrounding water molecules at different temperatures, the initial dissociation process is similar. On the other hand, the dissociated fragments fly away at room temperature, while they remain in the same cage when frozen water is used. Our results provide direct support for the cage effect mechanism.

Oligo-Adenine and Cyanuric Acid Supramolecular DNA-Based Hydrogels Exhibiting Acid-Resistance and Physiological pH-Responsiveness.

Expanding the functions and applications of DNA by integrating noncanonical bases and structures into biopolymers is a continuous scientific effort. An adenine-rich strand (A-strand) is introduced as functional scaffold revealing, in the presence of the low-molecular-weight cofactor cyanuric acid (CA, pKa 6.9), supramolecular hydrogel-forming efficacies demonstrating multiple pH-responsiveness. At pH 1.2, the A-strand transforms into a parallel A-motif duplex hydrogel cross-linked by AH+-H+A units due to the protonation of adenine (pKa 3.5). At pH 5.2, and in the presence of coadded CA, a helicene-like configuration is formed between adenine and protonated CA, generating a parallel A-CA triplex cross-linked hydrogel. At pH 8.0, the hydrogel undergoes transition into a liquid state by deprotonation of CA cofactor units and disassembly of A-CA triplex into its constituent components. Density functional theory calculations and molecular dynamics simulations, supporting the structural reconfigurations of A-strand in the presence of CA, are performed. The sequential pH-stimulated hydrogel states are rheometrically characterized. The hydrogel framework is loaded with fluorescein-labeled insulin, and the pH-stimulated release of insulin from the hydrogel across the pH barriers present in the gastrointestinal tract is demonstrated. The results provide principles for future application of the hydrogel for oral insulin administration for diabetes.

Ab initio exploration of modified carbon nanotubes as potential corrosion inhibitors

In order to develop novel unexplored potential corrosion inhibitors, covalently modified single-walled carbon nanotubes (SWCNT) via benzoic (–PhCOOH) and aniline (–PhNH2) groups are being investigated as corrosion inhibitors for the first time. Utilizing a comprehensive approach, this study employed density functional theory (DFT), Monte Carlo (MC), and molecular dynamics simulations (MD) to assess the adsorption behavior of modified nanotubes as corrosion inhibitors on the Cu(111) surface within a simulated aqueous HCl corrosion medium. The results provided molecular information on the adsorption capability, geometry adsorption centers, and adsorption energies (Eads) of carbon nanotubes on the surface of Cu(111). The adsorption energy values unveiled robust interactions between SWCNT–PhCOOH and SWCNT–PhNH2 inhibitors and the Cu(111) surface, suggesting a highly effective corrosion protection mechanism. The calculated Eads values exhibited notable ranges, spanning from –260.82 to –308.18 kcal/mol for SWCNT–PhCOOH and –220.92 to –261.01 kcal/mol for SWCNT–PhNH2 with the maximum probability values, representing the most favorable adsorption scenarios, determined to be –292.96 and –229.39 kcal/mol, respectively. A key insight from Monte Carlo simulations underscored the inherent spontaneity of the adsorption process, corroborated by the consistently negative Eads values. These findings collectively underscore the substantial affinity of the inhibitors to the copper surface, contributing to a deeper comprehension of their corrosion inhibition capabilities.

The inhibitor activity of some azo compound derivatives using density functional theory and molecular dynamics simulations

Corrosion poses a significant challenge for metals and alloys, carrying considerable economic and safety consequences. Various manufacturing processes, including surface preparation techniques, further exacerbate it. While corrosion inhibitors are employed to safeguard these materials from corrosion, they can also induce oxidative damage. This study aims to determine the corrosion inhibition efficiency of the specified materials. Density Functional Theory (DFT) and molecular dynamics simulations (MDS) were used to investigate the electronic characteristics of diverse compounds. The study included key parameters including highest occupied molecular orbital energy (EHOMO), lowest unoccupied molecular orbital energy (ELUMO), energy gap (Eg) between ELUMO and EHOMO, dipole moment, global hardness, softness (σ), ionization energy (I), electron affinity (A), electronegativity (χ), back-donation energy Eb-d, global electrophilicity (ω), electron transfer, global nucleophilicity (ε), and total energy (sum of electronic and zero-point energies). The calculations were executed utilizing the DFT method, employing a 6-311G++(d,p) basis set at the B3LYP level. Our research aimed to rank materials based on their effectiveness in inhibiting corrosion by considering EHOMO, ELUMO, and Eg. The order of effectiveness was found to be as follows: S4 > S3 > S1 > Mad1 > Mad3 > Mad4 > P-azo > Bad4. Additionally, our Investigation delved into alterations in electron affinity (A) and ionization potential (I) within different solvent phases, highlighting their importance in determining the chemical reactivity and stability of the studied compounds. The reactive sites were selected using Fukui indices, and the molecular electrostatic potential (MESP) surface was characterized—a Monte Carlo simulation to further understand these chemicals' adsorption behavior on Fe (110) surfaces. Our findings revealed that these compounds exhibit corrosion inhibition potential due to their high EHOMO, A, σ, ΔN, ΔE-back-donation, and low Eg, ELUMO, I, and η. Moreover, our MESP surface analysis demonstrated their ability to donate electrons to the metal's d-orbitals and receive electrons through back-donation. The compounds investigated in this study exhibit the potential to be employed as corrosion inhibitors in various industrial sectors. As a result, through the use and determination of quantum computational chemistry parameters, the level of corrosion resistance was as follows S3 > S1 > Bad4 > Mad4 > Mad3 > S4 > Mad1 > P-azo.

Valorizing the potential of saffron petals extract for aluminium corrosion control: An integrated approach involving extraction, experimental and computational analysis

Our study focuses on enhancing the antioxidant potential of saffron petals extract (SPE) as a sustainable corrosion inhibitor for aluminium (Al) in a 3.5% NaCl solution. The phytochemical screening and antioxidant activity are identified. Various investigative methodologies, including experimental techniques such as Electrochemical Impedance Spectroscopy (EIS) and Potentiodynamic Polarization (PDP), were employed to evaluate the effectiveness of this compound in inhibiting corrosion. Computational analysis, involving Density Functional Theory (DFT) and Molecular Dynamics (MD) simulations, were utilized to carry out a Theoretical study on the main constituents of SPE. The surface characteristics and composition of the corroded Al were examined using a combination of Scanning Electron Microscopy (SEM), Energy-Dispersive X-ray Spectroscopy (EDX) and Fourier Transform Infrared Spectroscopy (FT-IR). The Electrochemical results suggest that the effectiveness of the inhibitor was dependent on its concentration, reaching 95.11% at 500 ppm of SPE. The PDP indicates that mixed inhibition control effectively slowed down the corrosion of Al. Furthermore, the adsorption behaviour of SPE onto Al follows to the Langmuir isotherm model. SEM, EDX and FT-IR analyses further verified significant alterations in the surface morphology and roughness of Al and confirming the successful formation of a protective layer on the Al surface. Ultraviolet-Visible spectroscopy (UV-Vis) reveals evidence of interaction between Al and the SPE molecules. Computational studies substantiated these observations, demonstrating the reactivity and adsorption patterns of the evaluated major bioactive molecules of SPE on the Al surface.

Experimental and theoretical insights into Artemisia Stems aqueous extract as a sustainable and eco-friendly corrosion inhibitor for mild steel in 1M HCl environment.

The present research demonstrates an innovative investigation of environmentally friendly mild steel (M-steel) corrosion inhibition using the artemisia stems aqueous extract (ASAEx) as an inhibitor in hydrochloric acid 1M. The standard extraction technique of hydrodistillation was used for producing the aqueous solutions of ASAEx. To assess the ratios of the chemical components, phytochemical screening was used to identify the stems of this plant. We used a variety of methods and techniques in our research on corrosion inhibition, including weight loss measures, surface analysis methods like XPS and SEM/EDS, electrochemical testing like PDP and EIS, as well as computational lead compound evaluation. Maximum inhibitory efficacy was achieved with 400mg/L ASAEx in 1M HCl at 303K, i.e. 90%. The PDP investigation verified the mixed-kind inhibitor status of the ASAEx extract. To describe the surface of M-steel, fitting and synthetic data were used to identify a constant phase element (CPE). SEM surface analysis was also used to detect the ASAEx effect on the surface of M-steel. X-ray photoelectron spectroscopy (XPS) analysis shows the presence of trace molecules of ASAEx on M-steel surface characterizing the bands in Maj-ASAEx (major compound of ASAEx). Density functional theory (DFT) and molecular dynamics simulations (MDs) were used in computational chemistry to clarify the adsorption mechanism and inhibitory impact.

Modulating Thermal Conductivity via Targeted Phonon Excitation.

Thermal conductivity is a critical material property in numerous applications, such as those related to thermoelectric devices and heat dissipation. Effectively modulating thermal conductivity has become a great concern in the field of heat conduction. Here, a quantum modulation strategy is proposed to modulate the thermal conductivity/heat flux by exciting targeted phonons. It shows that the thermal conductivity of graphene can be tailored in the range of 1559 W m-1 K-1 (decreased to 49%) to 4093 W m-1 K-1 (increased to 128%), compared with the intrinsic value of 3189 W m-1 K-1. The effects are also observed for graphene nanoribbons and bulk silicon. The results are obtained through both density functional theory calculations and molecular dynamics simulations. This novel modulation strategy may pave the way for quantum heat conduction.

Alkyl chain length affecting the structure and interactions of choline-based amino acid ionic liquids: Insights from a computational study

The influences of the alkyl chain length of amino acid anions on the structure and interactions of choline-based amino acid ([Cho][AA]) ionic liquids were investigated based on density functional theory (DFT) and molecular dynamics (MD) simulations. The most stable configuration of the individual [Cho][AA] ion pair was determined, revealing negligible variations in the binding energy between choline cation ([Cho]+) and different amino acid anions ([AA]−). To further characterize the intermolecular interactions between charged entities, the atoms in molecules theory (AIM) and independent gradient model based on Hirshfeld partition (IGMH) were employed, which highlight the presence of robust hydrogen bonds between [Cho]+ and [AA] −, underscoring their significance in shaping the structure of [Cho][AA] ion pairs. The MD studies exhibited the aggregation of anions in [Cho][AA] ILs, which was attributed to side chain interactions of anions. This aggregation effect becomes more pronounced with the elongation of anionic alkyl chain. Simultaneously, radial distribution functions (RDF) analysis revealed a higher probability of hydrogen bonds occurring between [Cho]+ and [AA]− as the alkyl chain in [AA]− lengthened. Additionally, the hydrogen bond lifetimes roughly display an increasing trend as the alkyl chain of the [AA]− increases, which may be due to the reduced diffusion coefficients of ions. Hopefully, our results would contribute fundamental knowledge regarding the impact of alkyl chain length of anions on the structure and interactions in the choline-based amino acid ionic liquids.