Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “first principles density functional theory”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Structure and Oxygen Evolution Activity of β-NiOOH: Where Are the Protons?

Ni oxides and oxyhydroxides (NiO x ) have been studied for a long time as cathode materials for alkaline batteries and electrocatalysts for the oxygen evolution reaction (OER). Yet, understanding of the connection between their atomic and electronic structures and electrochemical performance or stability is still incomplete. In this work, we use first-principles density functional theory (DFT) calculations to revisit the structure, electronic properties, and OER activity of β-NiOOH, the catalytically active phase of NiO x . Following extensive DFTbased screening, we identify a hitherto overlooked structure characterized by a uniform distribution of H atoms on the NiO 2 layers. All the Ni 3+ cations in this structure exhibit an identical t g 6 e g 1 electronic configuration with an occupied 3d z2 orbital. Comparison of the calculated bond lengths with extended X-ray absorption fine structure (EXAFS) data unequivocally supports this structure relative to all other low-energy configurations. Based on this structure, we uncover and detail defect-dominated OER mechanisms on the basal β-NiOOH (001) surface, with overpotentials as low as 0.39 V. The present results should provide a valuable contribution to ongoing efforts for understanding and developing enhanced transition-metal hydroxide catalysts for the OER.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Phase diagram of magnetic shape memory alloy Ni 50 Mn $50–x$ In $x$ , 0 < $x$ , 25 from first principles, via spin cluster expansion and phonon vibrational entropies

The metamagnetic shape memory Heusler alloy Ni 50 Mn $50–x$ In $x$ exhibits a rich phase diagram featuring competing magnetic states, coupled magnetic–structural phase transitions, and strong compositional sensitivity. Existing first-principles approaches struggletocapturetheintertwinedchemical, magnetic, andvibrationaleffectsinthesealloys, necessitating a more integrated modeling framework. We develop a spin cluster expansion (spin-CE) framework augmented by a quasi-harmonic phonon model to capture both configurational (chemical and magnetic) and vibrational contributions to the free energy of Ni 50 Mn $50–x$ In $x$ over the full compositional range 0 ≤ x ≤25. The spin-CE includes both chemical clusters and composition-dependent Ising spin interactions, with parameters fit to a first-principles density functional theory (DFT) dataset. Using this approach, we predict the complete magnetostructural phase diagram and transformation temperatures of Ni 50 Mn $50–x$ In $x$ across the composition space. We find that vibrational entropy alone is insufficient to reproduce the martensitic transformation in the magnetic shape memory alloy regime, highlighting the essential role of magnetism. Incorporating both magnetic and vibrational contributions allows us to reproduce all experimentally known phases, including the disappearance of the stable martensite phase at a critical In concentration and the Curie temperature of the austenite phase. The method also captures the transition with increasing In in martensite from antiferromagnetic to ferromagnetic order and predicts re-entrant ferromagnetism, though the latter occurs at higher In content than reported experimentally. We discuss possible sources of this discrepancy and highlight the broader applicability of the method to other magnetostructurally complex materials, where it may offer mechanistic insight and predictive design capabilities.

Cluster expansion↗
DOE OSTI2021

Navigating Transition-Metal Chemical Space: Artificial Intelligence for First-Principles Design

Conspectus The variability of chemical bonding in open-shell transition-metal complexes not only motivates their study as functional materials and catalysts but also challenges conventional computational modeling tools. Here, tailoring ligand chemistry can alter preferred spin or oxidation states as well as electronic structure properties and reactivity, creating vast regions of chemical space to explore when designing new materials atom by atom. Although first-principles density functional theory (DFT) remains the workhorse of computational chemistry in mechanism deduction and property prediction, it is of limited use here. DFT is both far too computationally costly for widespread exploration of transition-metal chemical space and also prone to inaccuracies that limit its predictive performance for localized d electrons in transition-metal complexes. These challenges starkly contrast with the well-trodden regions of small-organic-molecule chemical space, where the analytical forms of molecular mechanics force fields and semiempirical theories have for decades accelerated the discovery of new molecules, accurate DFT functional performance has been demonstrated, and gold-standard methods from correlated wavefunction theory can predict experimental results to chemical accuracy. The combined promise of transition-metal chemical space exploration and lack of established tools has mandated a distinct approach. In this Account, we outline the path we charted in exploration of transition-metal chemical space starting from the first machine learning (ML) models (i.e., artificial neural network and kernel ridge regression) and representations for the prediction of open-shell transition-metal complex properties. The distinct importance of the immediate coordination environment of the metal center as well as the lack of low-level methods to accurately predict structural properties in this coordination environment first motivated and then benefited from these ML models and representations. Once developed, the recipe for prediction of geometric, spin state, and redox potential properties was straightforwardly extended to a diverse range of other properties, including in catalysis, computational “feasibility”, and the gas separation properties of periodic metal–organic frameworks. Interpretation of selected features most important for model prediction revealed new ways to encapsulate design rules and confirmed that models were robustly mapping essential structure–property relationships. Encountering the special challenge of ensuring that good model performance could generalize to new discovery targets motivated investigation of how to best carry out model uncertainty quantification. Distance-based approaches, whether in model latent space or in carefully engineered feature space, provided intuitive measures of the domain of applicability. With all of these pieces together, ML can be harnessed as an engine to tackle the large-scale exploration of transition-metal chemical space needed to satisfy multiple objectives using efficient global optimization methods. In practical terms, bringing these artificial intelligence tools to bear on the problems of transition-metal chemical space exploration has resulted in ML-model assessments of large, multimillion compound spaces in minutes and validated new design leads in weeks instead of decades.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Vibrational Entropy Contribution to Mixing Free Energy of Ni-Rich LiNi 1- y Co y O 2

First principles density-functional-theory calculations are performed of the mixing free energy of the layered oxide LiNi 1-y Co y O 2 , a model Li-ion-battery cathode material, at 1000 K. Analysis of internal energy, LN exchange, configurational entropy and vibrational entropy are combined to obtain the net mixing free energy. The vibrational entropy requires a more indirect treatment than the other terms; analysis based on heat capacity is presented here, to show that the vibrational entropy is of comparable magnitude to the configurational entropy, and it tends to promote solid-solution behavior in LiNi 1-y Co y O 2 .

25 ENERGY STORAGE↗
DOE OSTI2025

Coverage-dependent structures and thermodynamic stability of intercalated Gd layers beneath buffer-layer graphene on SiC(0001)

Electronic properties of two-dimensional (2D) materials are strongly influenced by their atomic arrangements, making the theoretically-aided characterization of experimentally-synthesized 2D structures crucial. Using first-principles density functional theory, we analyze nearly 200 configurations of intercalated Gd layers beneath buffer-layer graphene on SiC(0001) over a Gd coverage range of 0.01 < θ < 1.2. By fully relaxing selectively-constructed configurations at each coverage within a large, low-strain supercell, we determine the coverage dependence of the chemical potential for intercalated Gd structures. Thermodynamically-preferred configurations below θ ≈ 0.8 form single-atom-thick monolayers, while 3D-like or multilayer structures emerge beyond θ ≈ 0.9. Most structures are amorphous-like, including the configuration at the chemical potential minimum around θ ≈ 0.4. In contrast, a strongly stretched Gd(0001)-like monolayer at θ = 1/3 and a nearly perfect Gd(0001) monolayer at θ = 1 are significantly less favorable with 0.16 eV and 0.82 eV higher chemical potentials above the minimum, respectively. Furthermore, the graphene layer decoupled by intercalated Gd near the chemical potential minimum is significantly flatter compared to its morphology above intercalated 3D structures at higher coverages and nearly isolated Gd atoms in the lowest coverage region. In conclusion, these findings align with our experimental results and underscore the need for further research on this unique intercalated system, which holds significant potential for diverse applications.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Theoretical study of the crystal and electronic properties of α - RuI 3

The material α-RuCl 3 , with a two-dimensional Ru honeycomb sublattice, has attracted considerable attention because it may be a realization of the Kitaev quantum spin liquid. Recently, a new honeycomb material, α-RuCl 3 , was prepared under moderately high pressure, and it is stable under ambient conditions. However, different from α-RuCl 3 , α-RuI 3 was reported to be a paramagnetic metal without long-range magnetic order down to 0.35 K. Here, the structural and electronic properties of the quasi-two-dimensional α-RuI 3 are theoretically studied. First, based on first-principles density functional theory calculations, the ABC stacking honeycomb-layer R$\bar{3}$ (No. 148) structure is found to be the most likely stacking order for α-RuI 3 along the c axis. Furthermore, both R$\bar{3}$ and P$\bar{3}$1c are dynamically stable because no imaginary frequency modes were obtained in the phononic dispersion spectrum without Hubbard U. Moreover, the different physical behavior of α-RuI 3 compared to α-RuCl 3 can be understood naturally. The strong hybridization between Ru 4d and I 5p orbitals decreases the “effective” atomic Hubbard repulsion, leading the electrons of RuI3 to be less localized than in RuCl 3 . As a consequence, the effective electronic correlation is reduced from Cl to I, leading to the metallic nature of α-RuI 3 . Based on the DFT+U (Ueff=2 eV) plus spin-orbital coupling, we obtained a spin-orbit Mott insulating behavior for α-RuCl 3 and, with the same procedure, a metallic behavior for α-RuI 3 , in good agreement with experimental results. Furthermore, when introducing large (unrealistic) U eff =6 eV, the spin-orbit Mott gap opens in α-RuI 3 as well, supporting the physical picture we are proposing. Our results provide guidance to experimentalists and theorists working on two-dimensional transition metal tri-iodide layered materials.

2-dimensional systems↗
DOE OSTI2025

A Panoramic View of MXenes via an Atomic Coordination‐Based Design Strategy

Two‐dimensional (2D) transition metal carbides and nitrides, known as MXenes, possess unique physical and chemical properties, enabling diverse applications in fields ranging from energy storage to communication, catalysis, sensing, healthcare, and beyond. Despite extensive research and notable advancements, a fundamental understanding of MXenes’ phase diversity and its connection to their hierarchical precursors, including the intermediate MAX phases and the ancestral bulk phases, remains limited. Here, in this study, it is hypothesized that the atomic coordination environments adopted by transition metal and nonmetallic atoms in their three‐dimensional (3D) bulk precursors may persist in 2D MXenes to govern their phase diversity. Using high‐throughput modeling based on first‐principles density functional theory, a wide range of MXene phases is unveiled and comprehensively evaluate their relative stabilities across a large chemical space. The key to the approach lies in considering various atomic coordination environments drawn from four types of ancestral bulk phases. Through this comprehensive structural library of MXenes, general guiding principles are uncovered, such as a close alignment between the phase stability of MXenes and that of their 3D precursors. These findings introduce a new design strategy in which the atomic coordination environments in bulk phases can serve as reliable predictors for accessing the diverse structural landscape of MXenes.

MXenes↗
DOE OSTI2025

Nanometer Resolution Structure‐Emission Correlation of Individual Quantum Emitters via Enhanced Cathodoluminescence in Twisted Hexagonal Boron Nitride

Understanding the atomic structure of quantum emitters, often originating from point defects or impuritie, is essential for designing and optimizing materials for quantum technologies such as quantum computing, communication, and sensing. Despite the availability of atomic-resolution scanning transmission electron microscopy and nanoscale cathodoluminescence microscopy, experimentally determining the atomic structure of individual emitters is challenging due to the conflicting needs for thick samples to generate strong cathodoluminescence signals and thin samples for structural analysis. To overcome this challenge, significantly enhanced cathodoluminescence at twisted interfaces is leveraged to achieve sub-nanometer localization precision for the first time in mapping individual quantum emitters in carbon-implanted hexagonal boron nitride. This unprecedent spatial sensitivity, together with correlative electron energy loss spectroscopy quantitative scanning transmission electron microscopy imaging, and first principles density functional theory calculations, enables the identification of the atomic structure of the 440 nm blue emitter in hexagonal boron nitride as a substituted vertical carbon dimer. Building on the atomic structure insights, nanoscale spatially precise creation of blue emitters is demonstrated by electron beam irradiation of carbon-coated hexagonal boron nitride. This advancement in correlating atomic structures with optical properties lays the foundation for a deeper understanding and precise engineering of quantum emitters, significantly advancing the development of cutting-edge quantum information technologies.

2D material↗
DOE OSTI2021

Thickness and Spin Dependence of Raman Modes in Magnetic Layered Fe 3 GeTe 2

2D layered Fe 3 GeTe 2 has attracted increasing attention due to its high magnetic ordering temperature and novel physical properties. Lattice dynamics is a fundamental property of Fe 3 GeTe 2 , and its relationships with the number of layers and interlayer spin ordering have not yet been explored in depth. Here, by first-principles density functional theory calculations, the phonon vibrations and Raman intensities of Fe 3 GeTe 2 are systematically studied from the bulk to monolayer structures. Furthermore, the spin-phonon coupling effect is investigated by considering different interlayer magnetic orderings: ferromagnetic and antiferromagnetic. It is found that the frequencies of Raman modes in Fe 3 GeTe 2 exhibit considerable dependence on the layer number and spin ordering. In this work, the results not only reveal the notable spin-phonon interactions in Fe 3 GeTe 2 , but also demonstrate that Raman modes can be utilized for characterizing the sample thickness and interlayer spin ordering in this 2D magnet.

2D magnets↗
DOE OSTI2022

The Importance of Avoided Crossings in Understanding High Valley Degeneracy in Half‐Heusler Thermoelectric Semiconductors

Abstract Half‐Heusler (hH) compounds are promising candidates for inexpensive, low‐toxicity thermoelectric materials. It is well known that engineering electronic bands with high valley degeneracy is an effective approach for enhancing the performance of thermoelectric materials, and there are several routes for achieving high valley degeneracy in hH systems. For instance, there are multiple locations in the first Brillouin zone where the valence band maximum can be found (at the Γ‐, L‐, or W‐point), and there are two competing low‐lying conduction bands at the X‐point, where the conduction band minimum is located. By converging the multiple valence band and conduction band extrema, the valley degeneracy, and hence, performance of these materials can be improved. Here, group theoretical and tight‐binding approaches, in addition to first‐principles density functional theory calculations, are used to study the chemical origins of various band extrema in both the n‐type and p‐type compounds, with particular focus on ZrNiSn and NbFeSb. Specifically, the importance of avoided crossings is explained. The results of this work can be used to better understand and develop design strategies for engineering better performing hH thermoelectrics.

Brod, Madison K.↗
DOE OSTI2024

Electronic Properties of Ultra‐Wide Bandgap B x Al 1− x N Computed from First‐Principles Simulations

Abstract Ultra‐wide bandgap (UWBG) materials such as AlN and BN hold great promise for future power electronics due to their exceptional properties. They exhibit large bandgaps, high breakdown fields, high thermal conductivity, and high mechanical strengths. AlN and BN have been extensively researched, however, their alloys, B x Al 1− x N, are much less studied despite their ability to offer tunable properties by adjusting x . In this article, the electronic properties of 17 recently predicted ground states of B x Al 1− x N in the x = 0 − 1 range are predicted using first‐principles density functional theory and many‐body perturbation theory within GW approximation. All the B x Al 1− x N structures are found to be UWBG materials and have bandgaps that vary linearly from that of wurtzite‐phase ( w ) AlN (6.19 eV) to that of w ‐BN (7.47 eV). The bandstructures of B x Al 1− x N show that a direct‐to‐indirect bandgap crossover occurs near x = 0.25. Furthermore, it is found that B x Al 1− x N alloys have much larger dielectric constants than the constituent bulk materials (AlN = 9.3 ɛ 0 or BN = 7.3 ɛ 0 ), with values reaching as high as 12.1 ɛ 0 . These alloys are found to exhibit large dielectric breakdown fields in the range 9–35 MV cm −1 with a linear dependence on x . This work provides the much needed advancement in the understanding of the properties of B x Al 1− x N to aid their application in next‐generation devices.

Milne, Cody L.↗
DOE OSTI2022

Single Atoms Anchored in Hexagonal Boron Nitride for Propane Dehydrogenation from First Principles

We report single-atom catalysts embedded in N-doped graphene have attracted great interest recently, but the hexagonal boron nitride (h-BN) is much less explored as a support. Using first principles density function theory and molecular dynamics, here we investigate the stability of Pt, Au, and Ru single atoms anchored at B and N vacancies on h-BN. We find that Pt and Ru single atoms are much more stable than Au on h-BN. We further examine propane dehydrogenation on these single-atom catalysts and find that Pt 1 at the B vacancy in h-BN and Ru 1 at the N vacancy in h-BN show excellent activity for propane dehydrogenation, as evidenced by low energy barriers for both dehydrogenation steps. Our work suggests that Pt and Ru single atoms anchored at vacancy sites in h-BN could be promising for propane dehydrogenation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Understanding the Role of CdTe in Polycrystalline CdSe x Te 1– x /CdTe-Graded Bilayer Photovoltaic Devices

Grading of bandgap by alloying CdTe with selenium to form a CdSe x Te 1–x /CdTe-graded bilayer device has led to a device efficiency over 19%. A CdSe x Te 1–x absorber would increase the short-circuit current due to its lower bandgap but at the expense of open-circuit voltage. It has been demonstrated that adding a CdTe layer at the back of such a CdSe x Te 1–x film reduces the voltage deficit caused by the lower bandgap of absorber from selenium alloying while maintaining the higher short-circuit current. This leads to a photovoltaic device that draws advantage from both materials with an efficiency greater than either of them. In this report there is a detailed account using device data, ultraviolet photoelectron spectroscopy, electron microscopy, and first-principles density functional theory modeling is provided, which shows that CdTe acts as an electron reflector for CdSe x Te 1–x .

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Quaternary i-MAX Phases (Mo 2/3 RE 1/3 ) 2 AlC (RE: Dy, Tb, Er): Experimental Characterization and First-Principles Insights into their Fundamental Properties

Rare earth (RE)-based materials have unique electronic, magnetic, and optical properties, leading to the recent discovery of atomically layered solids with the chemical formula (M' 2/3 RE 1/3 ) 2 AlC, which have since garnered significant attention in the scientific community. This study aims to synthesize, characterize, and investigate the structural and thermal stability of the RE i-MAX phases. We prepared i-MAX phases using molybdenum (Mo) as M′ and RE elements as Dy, Tb, and Er, namely (Mo 2/3 Dy 1/3 ) 2 AlC, (Mo 2/3 Tb 1/3 ) 2 AlC, and (Mo 2/3 Er 1/3 ) 2 AlC. Structural characterization through x-ray diffraction (XRD) and Raman spectroscopy confirms the formation of the RE-based i-MAX phase, along with the presence of minor impurity phases in the alloys. Thermogravimetric analysis (TGA) conducted up to 1000°C under ambient conditions reveals that the i-MAX phases remain thermally stable up to approximately 450°C, beyond which oxidation leads to a noticeable weight gain in all samples. Differential scanning calorimetry (DSC) measurements during heating and cooling cycles show endothermic and exothermic peaks for (Mo 2/3 Dy 1/3 ) 2 AlC i-MAX in the 410–420°C range, indicating a temperature-induced minor atomic arrangement. In contrast, these peaks are absent in the Tb- and Er-based i-MAX phases. These findings offer valuable insights into the thermal behavior and stability of these i-MAX phases under thermal stress, contributing to a deeper understanding of their unique properties. Furthermore, first-principles density functional theory (DFT) calculations were performed to investigate the electronic and optical properties of the i-MAX phases. The results reveal their metallic nature, with pronounced contributions from Mo and RE elements near the Fermi level and within the conduction band.

Rare earth↗
DOE OSTI2020

First-principles study on the mechanical and thermodynamic properties of MoNbTaTiW

Refractory high-entropy alloys (RHEAs) are emerging as new materials for high temperature structural applications because of their stable mechanical and thermal properties at temperatures higher than 2273 K. In this study, the mechanical properties of MoNbTaTiW REDEA are examined by applying calculations based on first-principles density functional theory (DFT) and using a large unit cell with 100 randomized atoms. The phase calculation of MoNbTaTiW with CALPHAD method shows the existence of a stable body-centered cubic structure at a high temperature and a hexagonal closely packed phase at a low temperature. The predicted phase, shear modulus, Young’s modulus, Poisson’s ratio, and hardness values are consistent with available experimental results. The linear thermal expansion coefficient, vibrational entropy, and vibrational heat capacity of MoNbTaTiW RHEA are investigated in accordance with Debye-Grüneisen theory. These results may provide a basis for future research related to the application of RHEAs.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Effect of disorder on thermodynamic instability of binary Rare-earth – Nickel – Palladium compounds

In this work, we have investigated the thermodynamic stability of disordered rare-earth phases SmX 2 and Sm 10 X 21 (X=Ni, Pd) using machine-learning based analytical descriptor and first-principles density functional theory methods. The absence of Laves phase compounds in R-Pd binary systems is a longstanding problem of rare earth science: even though Ni and Pd belong to the same group of the periodic table and have similar electronic structure, the Pd compound crystallizes in a monoclinic (C2/m) phase with 10:21 stoichiometry, i.e., Sm 10 Pd 21 , while the Ni compound adopts a cubic Laves phase (MgCu 2 ) structure. To understand this contrasting phase stability, we performed thermodynamic convex hull analysis of Sm x Ni 1-x and Sm x Pd 1-x binary systems, which is experimentally validated using powder X-ray diffraction (PXRD) analyzes of polycrystalline Sm(Ni x Pd 1-x ) 2 samples with x=0, 0.5, and 1. A detailed electronic-structure (band-structure, charge density, and Fermi-surface) analysis of the differences between SmNi 2 /SmPd 2 and Sm 10 Ni 21 /Sm 10 Pd 21 compounds provides the quantum mechanical origin of the unfavorable mixing of Pd with Ni in cubic Laves phase. We show that the stability of Sm-Pd in 10:21 stoichiometry arises from improved intra-/inter-layer 5d-4d bonding compared to the 1:2 stoichiometry. Our work emphasizes the importance of ab-initio methods and computationally inexpensive analytical descriptors for the detailed analysis of thermodynamic and electronic properties of hard-to-prepare rare-earth compounds.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Doping transition-metal atoms in graphene for atomic-scale tailoring of electronic, magnetic, and quantum topological properties

Atomic-scale fabrication is an outstanding challenge and overarching goal for the nanoscience community. The practical implementation of moving and fixing atoms to a structure is non-trivial considering that one must spatially address the positioning of single atoms, provide a stabilizing scaffold to hold structures in place, and understand the details of their chemical bonding. Free-standing graphene offers a simplified platform for the development of atomic-scale fabrication and the focused electron beam in a scanning transmission electron microscope can be used to locally induce defects and sculpt the graphene. In this scenario, the graphene forms the stabilizing scaffold and the experimental question is whether a range of dopant atoms can be attached and incorporated into the lattice using a single technique and, from a theoretical perspective, we would like to know which dopants will create technologically interesting properties. Here in this paper, we demonstrate that the electron beam can be used to selectively and precisely insert a variety of transition metal atoms into graphene with highly localized control over the doping locations. We use first-principles density functional theory calculations with direct observation of the created structures to reveal the energetics of incorporating metal atoms into graphene and their magnetic, electronic, and quantum topological properties.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A deep learning interatomic potential developed for atomistic simulation of carbon materials

Interatomic potentials based on neural-network machine learning method have attracted considerable attention in recent years owing to their outstanding ability to balance the accuracy and efficiency in atomistic simulations. In this work, a neural-network potential (NNP) for carbon is generated to simulate the structural properties of various carbon structures. The potential is trained using a database consisting of crystalline and liquid structures obtained by the first-principles density functional theory (DFT) calculations. The developed potential accurately predicts the energies and forces in crystalline and liquid carbon structures, the energetic stability of defected graphene, and the structures of amorphous carbon as the function of density. As a result, the excellent accuracy and transferability of the NNP provide a promising tool for accurate atomistic simulations of various carbon materials with faster speed and much lower cost.

36 MATERIALS SCIENCE↗