First-principles Density Functional Theory Research Articles

Electromechanical strain response of phosphorene nanotubes

Nanomaterials that undergo structural or other property changes upon application of external stimuli are called stimuli responsive materials and are particularly suited for drug delivery, biosensing or artificial muscle applications. Two-dimensional (2D) black phosphorus is an ideal material for such applications due to its remarkable electromechanical response. Given that one-dimensional (1D) black phosphorus nanotubes (PNTs) are calculated to be energetically stable, it is possible that they can undergo similar electromechanical responses to their 2D counterparts, allowing their potential application as nanochannel devices for drug delivery. Using first-principles density functional theory, we investigated the electromechanical response of different-sized PNTs upon charge injection. Upon hole injection, the diameter of the PNTs expands up to a maximum of 30.2% for a (0,15) PNT that is 0.24 nm in diameter. The PNTs become highly p-doped as the valence band maximum crosses the Fermi level and undergoes switching from a direct to indirect band gap. The mechanism behind the electromechanical response was determined through analysis of the structural deformations, charge density distribution and Bader partial charges. It was shown that injection of charge alters the Young’s Modulus of the PNTs, as hole injection weakens the structural integrity of the nanotube, allowing a greater electromechanical response, with PNT-15 showing the largest decrease in the Young’s Modulus of 15.34%. These findings show that 1D PNTs are promising materials for the development of nanoelectromechanical actuators which could be used for drug delivery, energy harvesting or similar applications.

Charge Redistribution in Mechanochemical Reactions for Solid Interfaces.

Mechanochemical strategies are widely used in various fields, ranging from friction and wear to mechanosynthesis, yet how the mechanical stress activates the chemical reactions at the electronic level is still open. We used first-principles density functional theory to study the rule of the stress-modified electronic states in transmitting mechanical energy to trigger chemical responses for different mechanochemical systems. The electron density redistribution among initial, transition, and final configurations is defined to correlate the energy evolution during reactions. We found that stress-induced changes in electron density redistribution are linearly related to activation energy and reaction energy, indicating the transition from mechanical work to chemical reactivity. The correlation coefficient is defined as the term "interface reactivity coefficient" to evaluate the susceptibility of chemical reactivity to mechanical action for material interfaces. The study may shed light on the electronic mechanism of the mechanochemical reactions behind the fundamental model as well as the mechanochemical phenomena.

Synthesis and characterization of low-dimensional N-heterocyclic carbene lattices.

The covalent interaction of N-heterocyclic carbenes (NHCs) with transition metal atoms gives rise to distinctive frontier molecular orbitals (FMOs). These emergent electronic states have spurred the widespread adoption of NHC ligands in chemical catalysis and functional materials. Although formation of carbene-metal complexes in self-assembled monolayers on surfaces has been explored, design and electronic structure characterization of extended low-dimensional NHC-metal lattices remains elusive. Here we demonstrate a modular approach to engineering one-dimensional (1D) metal-organic chains and two-dimensional (2D) Kagome lattices using the FMOs of NHC-Au-NHC junctions to create low-dimensional molecular networks exhibiting intrinsic metallicity. Scanning tunneling spectroscopy and first-principles density functional theory reveal the contribution of C-Au-C π-bonding states to dispersive bands that imbue 1D- and 2D-NHC lattices with exceptionally small work functions.

Seeing the Unseen: The Structural Influence of the Lone Pair Electrons in PbWO4.

Pair distribution function (PDF) analysis of the scheelite-type material PbWO4 reveals previously unidentified short-range structural distortions in the PbO8 polyhedra and WO4 tetrahedra not observed in the similarly structured CaWO4. These local distortions are a result of the structural influence of the Pb2+ 6s2 lone pair electrons. These are not evident from the Rietveld analysis of synchrotron X-ray or neutron powder diffraction data, nor do they strongly influence the X-ray PDF (XPDF). This illustrates the importance of neutron PDF (NPDF) in the study of such materials. First-principles density function theory (DFT) calculations show that the Pb2+ 6s2 electrons are hybridized with the O2- 2p electrons near the Fermi level. The presence of local-scale distortions has previously been neglected in studies of structure-functionality relationships in PbWO4 and other scheelite-structured photocatalytic materials, including BiVO4, and this observation opens new avenues for their optimization.

Prediction of giant anomalous Nernst effect in Sm(Co,Ni)5

Sm-Co bulk alloys are well-known permanent magnets having large remanent magnetizations and coercive forces and are widely used in many industrial products. Recently, a large transverse thermoelectric conversion was observed for SmCo5 over a wide temperature range in the absence of magnetic fields. The large thermoelectric conductivity tensors (αxy) was also confirmed by the first-principles density functional theory (DFT) calculations. In this study, we predicted further enhancement of the αxy by including Ni in Co site of SmCo5. We showed that the αxy of Sm(Co1−xNix)5 increases with increasing the Ni ratio and takes the maximum value αxy = 11.3 A K−1 m−1 around x = 0.08 at 300 K, which is about 77% enhancement of αxy = 6.4 A K−1 m−1 in SmCo5. We clarified that the band proximity points near the nodal line of Sm(Co0.92Ni0.08)5 are the main contributing factor to the large Berry curvature, providing the steep slope of the energy dependence in the anomalous Hall conductivity around the Fermi energy.

Understanding the phase stability in multi-principal-component AlCuFeMn alloy

Abstract Method(s) that can reliably predict phase evolution across thermodynamic parameter space, especially in complex systems are of critical significance in academia as well as in the manufacturing industry. In the present work, phase stability in equimolar AlCuFeMn multi-principal-component alloy (MPCA) was predicted using complementary first-principles density functional theory (DFT) calculations, and ab-initio molecular dynamics (AIMD) simulations. Temperature evolution of completely disordered, partially ordered, and completely ordered phases was examined based on Gibbs free energy. Configurational, electronic, vibrational, and lattice mismatch entropies were considered to compute the Gibbs free energy of the competing phases. Additionally, elemental segregation was studied using ab-initio molecular dynamics (AIMD). The predicted results at 300K align well with room-temperature experimental observations using x-ray diffraction, scanning and transmission electron microscopy on a sample prepared using commercially available pure elements. The adopted method could help in predicting plausible phases in other MPCA systems with complex phase stability.

The elastic, electronic, and optical properties of TinO2n-1: Studying of GGA(LDA)+U

Based on the first-principles of density functional theory, the elastic, electronic and optical properties of three titanium suboxides TinO2n-1(n=5，6 and 7) are calculated and analyzed using generalized gradient approximation (GGA+U) and local density approximation (LDA+U) under the pseudopotential of ultra-soft. The results show that all three titanium suboxides are elastically anisotropic, but Ti5O9 has the smallest anisotropy. All the three materials exhibit certain ductility and have certern hardness (∼10 GPa). Electronic structure analysis indicates that Ti5O9, Ti6O11 and Ti7O13 have half-metallic properties, while Ti5O9 also has magnetic property, that is mainly derived from the 3d orbit electrons of Ti atoms. Optical property analysis shows that in the far ultraviolet region (6–12 eV), the reflection coefficient, absorption coefficient and extinction coefficient form resonance.

Strategically designing and fabricating nitrogen and sulfur Co-doped g-C3N4 for accelerating photocatalytic H2 evolution

Doping engineering is an effective strategy for graphitic carbon nitride (g-C3N4) to improve its photocatalytic hydrogen evolution reaction (HER) performance. In this work, a novel nitrogen and sulfur co-doped g-C3N4 (N, S-g-C3N4) is elaborately designed on the basis of theoretical predictions of first-principle density functional theory (DFT). The calculated Gibbs free energy of adsorbed hydrogen (ΔGH*) for N, S-g-C3N4 at the N-doping active sites is extremely close to zero (0.01 eV). Inspired by the theoretical predictions, the N, S-g-C3N4 is successfully fabricated through ammonia-rich pyrolysis synthesis strategy, in which ammonia is in-situ obtained by pyrolyzing melamine. Subsequent characterizations indicate that the N, S-g-C3N4 possesses high specific surface area, outstanding light utilization, good hydrophilicity, and efficient carrier transfer efficiency. Consequently, the N, S-g-C3N4 displays an extremely high H2 evolution rate of 8269.9 µmol g−1 h−1, achieves an apparent quantum efficiency (AQE) of 3.24%, and also possesses outsatnding durability. Theoretical calculations further demonstrate that N and S dopants can not only introduce doping energy level to reduce the band gap, but also induce charge redistribution to facilitate hydrogen adsorption, thus promoting the photocatalytic HER process. Moreover, femtosecond transient absorption (fs-TA) spectroscopy further corroborates the efficient photogenerated carrier transport of N, S-g-C3N4. This research highlights a promising and reliable strategy to achieve superior photocatalytic activity, and exhibits significant guidance for precise designing high-efficiency photocatalysts.

Structure Low Dimensionality and Lone-Pair Stereochemical Activity: the Key to Low Thermal Conductivity in the Pb-Sn-S System.

Recently, metal sulfides have begun to receive attention as potential cost-effective materials for thermoelectric applications beyond optoelectronic and photovoltaic devices. Herein, based on a comparative analysis of the structural and transport properties of 2D PbSnS2 and 1D PbSnS3, we demonstrate that the intrinsic effects that govern the low lattice thermal conductivity (κL) of these sulfides originate from the combination of the low dimensionality of their crystal structures with the stereochemical activity of the lone-pair electrons of cations. The presence of weak bonds in these materials, responsible for phonon scattering, results in inherently low κL of 1.0 W/m K in 1D PbSnS3 and 0.6 W/m K in 2D PbSnS2 at room temperature. However, the nature of the thermal transport is quite distinct. 1D PbSnS3 exhibits a higher thermal conductivity with a crystalline-like peak at low temperatures, while 2D PbSnS2 demonstrates glassy thermal conductivity in the entire temperature range investigated. First-principles density functional theory calculations reveal that the presence of antibonding states below the Fermi level, especially in PbSnS2, contributes to the very low κL. In addition, the calculated phonon dispersions exhibit very soft acoustic phonon branches that give rise to soft lattices and very low speeds of sounds.

Influence of benzothiophene acceptor moieties on the non-linear optical properties of pyreno-based chromophores: first-principles DFT framework.

Herein, a series of heterocyclic organic compounds (PYFD1-PYFD7) are designed with different acceptor moieties at the terminal position of a reference compound (PYFR) for nonlinear optical (NLO) active materials. The optoelectronic characteristics of the designed chromophores were investigated using density functional theory (DFT) calculations with the M06/6-311G(d,p) functional. Frontier molecular orbital (FMO) analysis revealed a significant decrease in the energy of the band gaps (2.340-2.602 eV) for the derivatives as compared to the PYFR reference compound (3.12 eV). An efficient transfer of charge from the highest occupied molecular orbital (HOMO) to the lowest unoccupied molecular orbital (LUMO) was seen, which was further corroborated by the density of states (DOS) and transition density matrix (TDM) heat maps. The results of the global reactivity parameters (GRPs) indicated that all derivatives exhibited greater softness (σ = 0.384-0.427 eV) and lower hardness (η = 0.394-1.302 eV) as compared to PYFR, indicating a higher level of polarizability in the derivatives. Moreover, all of the derivatives showed significant findings in terms of nonlinear optical (NLO) results as compared to the reference chromophore. PYFD2 showed the most effective NLO response (α = 1.861 × 10-22 and βtot = 2.376 × 10-28 esu), including a lowered band gap of 2.340 eV, the maximum softness value of 0.4273 eV, and the lowest hardness value of 1.170 eV as compared to other chromophores. The incorporation of different acceptors and thiophene as a π-spacer in this structural alteration significantly contributed to achieving remarkable NLO responses. Therefore, our findings may motivate experimentalists to synthesize these designed NLO active materials for the current advanced technological applications.

Selective SERS Detection of TATB Explosives Based on Hydroxy-Terminal Nanodiamond-Multilayer Graphene Substrate.

Selective surface-enhanced Raman scattering (SERS) detection of target explosives with good reproducibility is very important for monitoring soldiers' health and ecological environment. Here, the specific charge transfer pathway was constructed between a stable nanodiamond-multilayer graphene (MGD) film substrate and the target explosives. Two-step wet chemical oxidation methods of H2O2 (30%) and HNO3 (65%) solutions were used to regulate the terminal structure of MGD films. The experimental results showed that the hydroxyl (-OH) functional groups are successfully modified on the surface of MGD thin films, and the MGD-OH substrates having good selectivity for 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) explosive in mixed solutions of the TATB, 2,2-dinitroethene-1,1-diamine, 2,4,6-trinitrotoluene, and 1,3,5-trinitroperhydro-1,3,5-triazine explosives compared with MGD substrates were demonstrated. Finally, first-principles density functional theory simulations revealed that the SERS enhancement of the MGD-OH substrate is mainly attributed to the transferred electrons between the -NO2 groups of TATB and the -OH groups of the MGD-OH substrate.

Investigation of structural, electronic, optical and thermoelectric properties of hallide double perovskite Rb2AlAgX6 (X = Cl, I): A first-principles DFT study

In order to meet the demands of the energy consumption, double perovskites (DP) are a dependable green energy source. Consequently, there is a great deal of promise for thermoelectric and optoelectronic device applications when studying these perovskite halides. Our work examined the thermoelectric and optoelectronic properties of Rb2AlAgCl6 and Rb2AlAgI6 halides using Density Functional Theory (DFT). The tolerance factor and energy of formation for the halides that were examined show that they are stable in the cubic phase, structurally and thermodynamically. We calculated the bandgaps for Rb2AlAgCl6 (Eg = 3.82 eV) and Rb2AlAgI6 (Eg = 1.321 eV) using mBJ potentials, and compared them to the reported values. The optical characteristics of the materials under investigation were revealed through the use of a complicated dielectric function. The computed optical findings show that the materials are suitable for optoelectronic applications, with maximal light absorption occurring in the ultraviolet (UV) and infrared (IR) ranges. Optical parameters obtained includes large values of absorption coefficients (∼105cm−1), low reflectivity (∼1–15 %), and high optical conductivity (∼1015 sec−1). Thermoelectric properties revealed enhanced power factor for Rb2AlAgI6 as compared to Rb2AlAgCl6. These results imply that these materials are appropriate for application in optical and thermoelectric devices.

Thermal expansion, lattice vibration, and isotope effect on hydrogen diffusion in BCC Fe, Cr, and W from first-principles calculations.

Fe, Cr, and W are important elements in the alloys of in-reactor materials and operate in high-temperature environments with thermal expansion. Their tritium-impeding abilities are crucial to the radiation safety of various nuclear reactors. In this study, first-principles density functional theory is combined with quasi-harmonic approximation to evaluate factors that can affect the interstitial formation energy and diffusion coefficient of hydrogen isotopes in body-centered cubic (BCC) Fe, Cr, and W, including thermal expansion, metal host lattice vibrations, phonon density-of-states (pDOS) coupling diffusing atoms, and isotope effects. Calculation results indicate that the interstitial formation energy decreases as lattice expansion increases, whereas the jump barriers remain almost constant. Thermal expansion, host lattice vibration, and pDOS coupling minimally affect the diffusion coefficients of hydrogen isotopes in Fe, Cr, and W. The diffusion coefficient ratios between hydrogen isotopes are higher than the inverse ratio of the square root of the isotope mass at low temperatures. However, they decrease to the inverse ratio of the square root of the isotope mass at temperatures exceeding 800K. This study comprehensively investigates factors that affect the diffusion coefficients of hydrogen isotopes in BCC Fe, Cr, and W, thus providing a firm theoretical foundation for predicting the diffusion coefficients of tritium at different temperatures using protium/deuterium diffusion coefficients.

Relationship between synthesis method-crystal structure-melting properties in cocrystals: the case of caffeine-citric acid.

The influence of the crystal synthesis method on the crystallographic structure of caffeine-citric acid cocrystals was analyzed thanks to the synthesis of a new polymorphic form of the cocrystal. In order to compare the new form to the already known forms, the crystal structure of the new cocrystal (C8H10N4O2·C6H8O7) was solved by powder X-ray diffraction thanks to synchrotron experiments. The structure determination was performed using `GALLOP', a recently developed hybrid approach based on a local optimization with a particle swarm optimizer, particularly powerful when applied to the structure resolution of materials of pharmaceutical interest, compared to classical Monte-Carlo simulated annealing. The final structure was obtained through Rietveld refinement, and first-principles density functional theory (DFT) calculations were used to locate the H atoms. The symmetry is triclinic with the space group P\overline{1} and contains one molecule of caffeine and one molecule of citric acid per asymmetric unit. The crystallographic structure of this cocrystal involves different hydrogen-bond associations compared to the already known structures. The analysis of these hydrogen bonds indicates that the cocrystal obtained here is less stable than the cocrystals already identified in the literature. This analysis is confirmed by the determination of the melting point of this cocrystal, which is lower than that of the previously known cocrystals.

BC cone-shaped anodes for lithium-ion batteries

The efficiency of lithium-ion batteries (LIBs) depends upon anode materials possessing high capacity. In present study, we investigate boron carbide cone anode (BCC) for LIBs. The first-principles density functional theory (DFT) method is used to design anode materials based on BCC for application in LIBs. We found that the BCC shows negative Li adsorption energies. Our results also revealed the fast diffusions of Li ion on the BCC with the small energy barrier of 0.34 eV. The storage capacity of BCC system rank among the highest capacity of anodes for LIBs, with value of 785 mAh/g. The low average open-circuit voltage (OCV) value of BCC system indicates that this anode can handle a large operating voltage in LIBs.

Linear photogalvanic effect in two-dimensional PSi photodetector by vacancy- and substitution-doping

Using particle swarm optimization methodology for crystal structure prediction and first-principles density functional theory, a graphene-like honeycombed PSi monolayer has been found by our group with its thermodynamical, dynamical and mechanical stability having been listed. To extend the application of the PSi monolayer, the armchair and zigzag photodetector devices based on the photogalvanic effect have been built. When the linearly polarized light illuminates on the armchair or zigzag photodetector, due to the C1 symmetry of PSi monolayer the PSi photodetector can directly produce high photocurrents. In addition, the produced photocurrents show the relations cos(2θ) and sin(2θ) on the polarization angle θ for the armchair and zigzag photodetector, respectively. Moreover, the value of photon energy can effectively influence the distribution of generated photocurrents, namely the relations on the angle θ. Especially, the vacancy- and substitution-doping in the PSi photodetector can enhance the generated photocurrents. These results demonstrated great potential applications of the PSi monolayer on the low energy-consumption optoelectronics devices.

Nanocrystal-based thermoelectric SnTe-NaSbSe2 alloys with strengthened band convergence and reduced thermal conductivity

Ternary I-V-VI2 colloidal NaSbSe2 nanocrystals are herein used to improve the performance of lead-free SnTe thermoelectric materials. We showcase a versatile bottom-up engineering approach to produce nanocrystal-based SnTe-NaSbSe2 alloys from the rapid hot press of colloidal nanocrystal building blocks. The incorporation of NaSbSe2 nanocrystals significantly enhances the Seebeck coefficient of SnTe. The band convergence and simultaneous increasing band gap of SnTe-NaSbSe2 alloys are certified by the first-principles density functional theory calculations. Besides, defect engineering generated by the incorporation of NaSbSe2 nanocrystals such as Sn vacancies, substitution point defects, dense dislocations, and strains generated by the NaSbSe2 nanoparticles incorporation result in a dramatic reduction of the lattice thermal conductivity below the amorphous limit of pure SnTe, down to 0.38 W m−1 K−1. As a consequence, power factors enhance up to 1.77 mW m−1 K−2, which is ∼193 % higher than that of the pristine SnTe, and thermoelectric figures of merit up to 1.15 at 823 K for (SnTe)0.85(NaSbSe2)0.15 are achieved.

Chlorinated graphene and graphene nanoribbons: A density functional theory study

Functional groups based on halides, such as fluorine (F), chlorine (Cl), bromine (Br), and iodine (I), are crucial for understanding the chemical reactivity of graphitic nanomaterials. Except for I, halogens exhibit electronegativity greater than carbon (C); therefore, charge transfer from carbon to halogen is expected. First-principles density functional theory calculations were performed to determine the role of different Cl-functional groups (methyl-trichloride, ethyl-trichloride, chloride, acyl-chloride, vinyl-chloride, acetyl hypochlorite, chloramines, sulfonyl chloride, and more) on the electronic properties of graphene and graphene nanoribbons (GNRs). GNRs with zigzag edges (ZGNRs) and armchair edges (AGNRs) were studied. We analyzed the optimized structures, band structure, density of states, cohesive energy, and band gap. Our results revealed that the based-Cl functional groups can provide an alternative route to activate the borders and surfaces of sp2 carbon materials. Methyl-trichloride and acyl-chloride can induce magnetism and metallicity. Chloride and acyl-chloride are the most energetically stable functional groups attached to the edges. Surprisingly, methyl-trichloride or acyl-chloride functionalizing the surface of the AGNRs showed a direct (indirect) band gap for states with spin-up (spin-down). The results of aromatic (chlorobenzene- and dichlorobenzene-like structures) functionalization considering F, Cl, Br, and I are also shown. Finally, –F2 and –ClF functionalization cases are discussed.

First-principles study of helium solution and diffusion in tungsten borides

The behavior of helium (He) in plasma facing materials (PFMs) is crucial for the performance of PFMs under neutron irradiation in a fusion environment. First-principles density functional theory calculations have been performed to model the solution and diffusion of He in six stoichiometric ground-state tungsten borides (WxBy), including W2B, WB, WB2, W2B5, WB3 and WB4. The solution energy of He interstitial in W borides linearly decreases with increasing B content, except for the deviation in WB and WB3. The mechanical and electronic contributions to He solubility in WxBy have been discussed. The calculated diffusion energetics indicates that He transport in W2B5, WB3 and WB4 is much more difficult than in W2B, WB and WB2. Subsequently, we investigate the influence of vacancies on the He behavior. Except for WB3, He is more strongly bound in W vacancies rather than B vacancies and the He binding energy to a W vacancy generally decreases with increasing B content in WxBy. These findings can be used to predict the potential influence of He on the properties of different W borides in the fusion environment.

Structural, Energetic, and Electronic Properties of H-Interstitial in C-Monovacancy: A First-Principles Density Functional Theory

We discover a unique structural-modified diamond that exhibits similar symmetry and band gap energy to that of the pure diamond. We study a complex carbon-vacancy-hydrogen in the diamond using the density-functional-theory method. The defective models are created by adding H-interstitial (Hi, where i = 1, 2, 3 and 4) in the 3D diamond C-monovacancy. The result shows that carbon-vacancy-hydrogen defects significantly decreased the symmetry from Td to C2V. Likewise, the volumetric size of the systems is widening up to 48.70 %, while the optimized band gap energies are narrowing. Additional states appeared in the C-monovacancy, H1-V, H2-V, and H3-V systems which further improved electron mobility. The Hi compensates the C-monovacancy which further serves as a deep donor. Interestingly, H4-V exhibits similar symmetry and band gap energy to that of the pure diamond, but its volumetric size is 48.70 % wider. HIGHLIGHTS Study of complex carbon-vacancy-hydrogen defects through density-functional-theory calculations Discover a unique structural-modified diamond that is 48.70 % wider than the pure diamond, but exhibits similar symmetry and band gap energy. Complex carbon-vacancy-hydrogen defects significantly lower the symmetries, widening the volumetric sizes, and narrowing the band gap energies. The H-interstitial acts as a deep donor. GRAPHICAL ABSTRACT