Adiabatic density-functional perturbation theory

Electron ptychography reveals a ferroelectricity dominated by anion displacements

Sodium niobate, a lead-free ferroic material, hosts delicately balanced, competing order parameters, including ferroelectric states that can be stabilized by epitaxial strain. Here we show that the resulting macroscopic ferroelectricity exhibits an unconventional microscopic structure using multislice electron ptychography. This technique overcomes multiple scattering artefacts limiting conventional electron microscopy, enabling both lateral spatial resolution beyond the diffraction limit and recovery of three-dimensional structural information. These imaging capabilities allow us to separate the ferroelectric interior of the sample from the relaxed surface structure and identify the soft phonon mode and related structural distortions with picometre precision. Unlike conventional ferroelectric perovskites, we find that the polar distortion in this material involves minimal distortions of the cation sublattices and is instead dominated by anion displacements relative to the niobium sublattice. We establish limits on film thickness for interfacial octahedral rotation engineering and directly visualize a random octahedral rotation pattern, arising from the flat dispersion of the associated phonon mode.

Nitrile Groups as Build-In Molecular Sensors for Interfacial Effects at Electrocatalytically Active Carbon-Nitrogen Materials

Electrocatalytic reactions are influenced by various interfacial phenomena including nonspecific interaction forces. For many examples, their contributions to the catalytic cycle have yet to be identified. Noncovalent interactions between the electrode and the electrolyte can be described by the local electric field environment at the interface and are experimentally accessible based on the Vibrational Stark Effect. We herein present a carbon-based C2N-type electrocatalyst that is active for the hydrogen evolution reaction and that contains nitrile functions as Stark reporter groups. With this system, we expand the range of electrocatalytically active systems suitable for electrochemical Stark spectroscopy while taking a step away from pure model systems. The stretching mode ν(C≡N) was analyzed via experimental and calculated Raman spectroscopy, revealing a defect character of the inherent CN groups. The ν(C≡N) peak position was furthermore studied via in situ electrochemical Raman spectroscopy. At noncatalytic conditions, a linear dependence between an applied electric potential and ν(C≡N) peak shift is observed, resulting in a red-shift at a more negative potential. At catalytic conditions, deviations from the linearity occur, and a semipermanent blue-shift of the CN peak is observed after electrocatalysis, implying a restructuring of the electrochemical double layer and therefore a change in the local electric field environment due to the catalytic turnover and the associated interfacial processes.

Impact of Microstructural Units Configuration and Arrangement on Phase-Matching in Nonlinear Optical KGa5Se8 Crystals

Achieving phase-matching (PM) behavior, which is typically derived from sufficient birefringence (Δn), is crucial for maximizing the nonlinear optics (NLO) laser power output. Δn is largely attributed to the configuration and arrangement of microstructural units, and phase transitions will significantly alter the stacking method of structural units. In this study, we successfully synthesized two diamond-like isomers of α-KGa5Se8 (P21, 1) and β-KGa5Se8 (P1, 2), which consist of open-honeycomb-like anionic frameworks interspersed with embedded K+ cations. Importantly, our analysis of frequency doubling experiment data revealed a significant finding at a laser wavelength of 1910 nm: the α-phase lacks PM capability, whereas the β-phase exhibits PM behavior. Notably, through structural and computational analysis, the distinct Δn primarily stems from the different arrangements of tetrahedral units within the two phases. Specifically, the polarizability anisotropy of the tetrahedral units in β-KGa5Se8 is 3.7 times that of α-KGa5Se8. This uniform packing improves the Δn index from 0.015 (α-phase) to 0.023 (β-phase). Furthermore, β-KGa5Se8 not only demonstrates a robust PM second-harmonic generation (SHG) response (2.1 × AgGaS2 @1910 nm) but also possesses a wide band gap of 2.56 eV, surpassing the 2.33 eV threshold that corresponds to the energy of doubly frequency light using the 1064 nm laser, a characteristic often not exceeding in most selenides. This work underscores the importance of the configuration and arrangement of microstructural units for achieving PM in NLO materials.

Ultrafast dynamics of hot carriers: Theoretical approaches based on real-time propagation of carrier distributions

In recent years, computational approaches which couple density functional theory (DFT)-based description of the electron-phonon and phonon-phonon scattering rates with the Boltzmann transport equation have been shown to obtain the electron and thermal transport characteristics of many 3D and 2D semiconductors in excellent agreement with experimental measurements. At the same time, progress in the DFT-based description of the electron-phonon scattering has also allowed to describe the non-equilibrium relaxation dynamics of hot or photo-excited electrons in several materials, in very good agreement with time-resolved spectroscopy experiments. In the latter case, as the time-resolved spectroscopy techniques provide the possibility to monitor transient material characteristics evolving on the femtosecond and attosecond time scales, the time evolution of photo-excited, nonthermal carrier distributions has to be described. Similarly, reliable theoretical approaches are needed to describe the transient transport properties of devices involving high energy carriers. In this review, we aim to discuss recent progress in coupling the ab initio description of materials, especially that of the electron-phonon scattering, with the time-dependent approaches describing the time evolution of the out-of-equilibrium carrier distributions, in the context of time-resolved spectroscopy experiments as well as in the context of transport simulations. We point out the computational limitations common to all numerical approaches, which describe time propagation of strongly out-of-equilibrium carrier distributions in 3D materials, and discuss the methods used to overcome them.

Smallest [5,6]Fullerene as Building Blocks for 2D Networks with Superior Stability and Enhanced Photocatalytic Performance

The assembly of molecules to form covalent networks can create varied lattice structures with physical and chemical properties distinct from those of conventional atomic lattices. Using the smallest stable [5,6]fullerene units as building blocks, various 2D C24 networks can be formed with superior stability and strength compared to the recently synthesized monolayer polymeric C60. Monolayer C24 harnesses the properties of both carbon crystals and fullerene molecules, such as stable chemical bonds, suitable band gaps, and large surface area, facilitating photocatalytic water splitting. The electronic band gaps of C24 are comparable to those of TiO2, providing appropriate band edges with sufficient external potential for overall water splitting over the acidic and neutral pH range. Upon photoexcitation, strong solar absorption enabled by strongly bound bright excitons can generate carriers effectively, while the type-II band alignment between C24 and other 2D monolayers can separate electrons and holes in individual layers simultaneously. Additionally, the number of surface-active sites of C24 monolayers are three times more than that of their C60 counterparts in a much wider pH range, providing spontaneous reaction pathways for the hydrogen evolution reaction. Our work provides insights into materials design using tunable building blocks of fullerene units with tailored functions for energy generation, conversion, and storage.

Achieving Phase-Matching in Nonlinear Optical Materials CsM2In2S6 (M = Cd/In, Hg/In) by the Incorporation of Unprecedented Trigonal Planar MS3 Motifs

The trigonal planar unit possesses significant hyperpolarizability and polarizability anisotropy, which makes it useful for optimizing nonlinear optical (NLO) materials, however, chalcogenide with this unit has seldom been reported. In this work, a novel approach is introduced by integrating the unprecedented trigonal planar MS3 (M = Cd/In, Hg/In) motifs into the nearly optically isotropic tetrahedral units, resulting in two novel chalcogenides CsM2In2S6 (M = Cd/In, 1; Hg/In, 2). Notably, structures 1 and 2 feature nearly planar triangular units at the center, encircled by three trimers, further interconnecting each other to create 3D frameworks. Importantly, phases 1 and 2 display phase-matching (PM) capabilities, primarily attributed to incorporating trigonal planar MS3 units that additionally enhance polarizability anisotropy. Furthermore, compounds 1 and 2 demonstrate moderate second-harmonic generation (SHG) signals (0.70 and 0.84 × AgGaS2@1.7 µm). This study pioneers an efficient strategy for the design of infrared NLO crystals with PM capabilities.

Achieving strong optical nonlinearity and wide bandgap of pnictides via ionic motif-driven directed assembly of covalent groups

Noncentrosymmetric (NCS) pnictides are indispensable for nonlinear optics, ferroelectrics, magnetic Weyl electronics, etc., areas, yet their structure design remains a substantial challenge. By using asymmetric ionic unit-driven covalent groups orienting and rigidity-flexibility coupling dual strategy, we successfully design and synthesize four NCS pnictides: [Sr4Br]2[MII3Si25P40] (MII = Mg, Cd) and [Ba3Br][MIIISi10P16] (MIII = Ga, In), which exhibit strong second harmonic generation effects (5.2 to 7.5 × AgGaS2), wide bandgaps (1.81 to 1.90 electron volts), and moderate birefringence (0.030 to 0.051). An unprecedented NCS structure-inducing mechanism analysis revealed that the (Sr4Br) and (Ba4Br) ionic units featuring the diamond-like electrostatic force field effectively break inversion symmetry and trigger uniform arrangement of the covalent tetrahedron groups. Furthermore, the nonlinear optical (NLO) properties and birefringence can be remarkably tuned by the secondary covalent building blocks (MII/IIIP4 tetrahedra) with distinct bond flexibility providing a broader space for regulating the key parameters. This work might expand chemical space for exploiting high-performance pnictide NLO materials.

Robust Computation and Analysis of Vibrational Spectra of Layered Framework Materials Including Host-Guest Interactions

Layered framework materials, a rapidly advancing class of porous materials, are composed of molecular components stitched together via covalent bonds and are usually synthesized through wet-chemical methods. Computational infrared (IR) and Raman spectra are among the most important characterization tools for this material class. Besides the a priori known spectra of the molecular building blocks and the solvent, they allow for in situ monitoring of the framework formation during synthesis. Therefore, they need to capture the additional peaks from host-guest interactions and the bands from emerging bonds between the molecular building blocks, verifying the successful synthesis of the desired material. In this work, we propose a robust computational framework based on ab initio molecular dynamics (AIMD), where we compute IR and Raman spectra from the time-correlation functions of dipole moments and polarizability tensors, respectively. As a case study, we apply our methodology to a covalent organic framework (COF) material, COF-1, and present its AIMD-computed IR and Raman spectra with and without 1,4-dioxane solvent molecules in its pores. To determine robust settings, we meticulously validate our model and explore how stacking disorder and different methods for computing dipole moments and polarizabilities affect IR and Raman intensities. Using our robust computational protocol, we achieve excellent agreement with experimental data. Furthermore, we illustrate how the computed spectra can be dissected into individual contributions from the solvent molecules, the molecular building blocks of COF-1, and the bonds connecting them.

Prediction of s^{±}-Wave Superconductivity Enhanced by Electronic Doping in Trilayer Nickelates La_{4}Ni_{3}O_{10} under Pressure

Motivated by the recently reported signatures of superconductivity in trilayer La_{4}Ni_{3}O_{10} under pressure, we comprehensively study this system using ab initio and random-phase approximation techniques. Without electronic interactions, the Ni d_{3z^{2}-r^{2}} orbitals show a bonding-antibonding and nonbonding splitting behavior via the O p_{z} orbitals inducing a "trimer" lattice in La_{4}Ni_{3}O_{10}, analogous to the dimers of La_{3}Ni_{2}O_{7}. The Fermi surface consists of three electron sheets with mixed e_{g} orbitals, and a hole and an electron pocket made up of the d_{3z^{2}-r^{2}} orbital, suggesting a Ni two-orbital minimum model. In addition, we find that superconducting pairing is induced in the s^{±}-wave channel due to partial nesting between the M=(π,π) centered pockets and portions of the Fermi surface centered at the Γ=(0,0) point. With changing electronic density n, the s^{±} instability remains leading and its pairing strength shows a domelike behavior with a maximum around n=4.2 (∼6.7% electron doping). The superconducting instability disappears at the same electronic density as that in the new 1313 stacking La_{3}Ni_{2}O_{7}, correlated with the vanishing of the hole pocket that arises from the trilayer sublattice, suggesting that the high-T_{c} superconductivity of La_{3}Ni_{2}O_{7} does not originate from a trilayer and monolayer structure. Furthermore, we confirm the experimentally proposed spin state in La_{4}Ni_{3}O_{10} with an in-plane (π, π) order and antiferromagnetic coupling between the top and bottom Ni layers, and spin zero in the middle layer.

Magnetic anisotropy and phononic properties of two-dimensional ferromagnetic Fe3GeS2 monolayer

In 2023, Fe3GeS2 monolayer with Curie temperature of 630 K is predicted, which is promising to be used in next-generation spintronic devices. However, its magnetic anisotropy and phononic properties are still unclear. In this paper, we implemented the first-principles calculations on Fe3GeS2 monolayer, and found its ferromagnetic ground state with robustness to the -1.5%-1.3% biaxial strain. Meanwhile, the out-of-plane magnetic anisotropy dominated by dipolar interaction is found in Fe3GeS2 monolayer. Finally, we studied the phononic properties to identify the dynamical stability of Fe3GeS2 monolayer and highlight the contribution from the anharmonic interaction of optical phonons to the thermal expansion coefficient. We also find two single-phonon modes can be used to design quantum mechanical resonators with a wide cool-temperature range. These results can provide a comprehensive understanding of the magnetism and phonon properties of two-dimensional (2D) Fe3GeS2, beneficial for the application of 2D Fe3GeS2 in spintronics.

Non-invasive detection of hazardous materials with a thermal-to-epithermal neutron station: a feasibility study towards practical application

The growing scale of the devastation that even a single terrorist attack can cause requires more effective methods for the detection of hazardous materials. In particular, there are no solutions for effectively monitoring threats at sea, both for the off-shore infrastructure and ports. Currently, state-of-the-art detection methods determine the density distribution and the shapes of tested subjects but only allow for a limited degree of substance identification. This work aims to present a feasibility study of the possible usage of several methods available on the thermal-to-epithermal neutron station, VESUVIO, at the ISIS neutron and muon spallation source, UK, for the detection of hazardous materials. To this end, we present the results of a series of experiments performed concurrently employing neutron transmission and Compton scattering using melamine, a commonly used explosive surrogate, in order to determine its signal characteristics and limits of detection and quantitation. The experiments are supported by first-principles modelling, providing detailed scrutiny of the material structure and the nuclear dynamics behind the neutron scattering observables.

Quantitative regulation of electron-phonon coupling

Electron-phonon (e-p) coupling plays a crucial role in various physical phenomena, and regulation of e-p coupling is vital for the exploration and design of high-performance materials. However, the current research on this topic lacks accurate quantification, hindering further understanding of the underlying physical processes and its applications. In this work, we demonstrate quantitative regulation of e-p coupling, by pressure engineering andin-situspectroscopy. We successfully observe both a distinct vibrational mode and a strong Stokes shift in layered CrBr3, which are clear signatures of e-p coupling. This allows us to achieve precise quantification of the Huang-Rhys factorSat the actual sample temperature, thus accurately determining the e-p coupling strength. We further reveal that pressure efficiently regulates the e-p coupling in CrBr3, evidenced by a remarkable 40% increase inSvalue. Our results offer an approach for quantifying and modulating e-p coupling, which can be leveraged for exploring and designing functional materials with targeted e-p coupling strengths.

Merits and Demerits of Machine Learning of Ferroelectric, Flexoelectric, and Electrolytic Properties of Ceramic Materials

In the present review, the merits and demerits of machine learning (ML) in materials science are discussed, compared with first principles calculations (PDE (partial differential equations) model) and physical or phenomenological ODE (ordinary differential equations) model calculations. ML is basically a fitting procedure of pre-existing (experimental) data as a function of various factors called descriptors. If excellent descriptors can be selected and the training data contain negligible error, the predictive power of a ML model is relatively high. However, it is currently very difficult for a ML model to predict experimental results beyond the parameter space of the training experimental data. For example, it is pointed out that all-dislocation-ceramics, which could be a new type of solid electrolyte filled with appropriate dislocations for high ionic conductivity without dendrite formation, could not be predicted by ML. The merits and demerits of first principles calculations and physical or phenomenological ODE model calculations are also discussed with some examples of the flexoelectric effect, dielectric constant, and ionic conductivity in solid electrolytes.

Structural phase transition, s±-wave pairing, and magnetic stripe order in bilayered superconductor La3Ni2O7 under pressure

Motivated by the recently discovered high-Tc superconductor La3Ni2O7, we comprehensively study this system using density functional theory and random phase approximation calculations. At low pressures, the Amam phase is stable, containing the Y2- mode distortion from the Fmmm phase, while the Fmmm phase is unstable. Because of small differences in enthalpy and a considerable Y2- mode amplitude, the two phases may coexist in the range between 10.6 and 14 GPa, beyond which the Fmmm phase dominates. In addition, the magnetic stripe-type spin order with wavevector (π, 0) was stable at the intermediate region. Pairing is induced in the s±-wave channel due to partial nesting between the M = (π, π) centered pockets and portions of the Fermi surface centered at the X = (π, 0) and Y = (0, π) points. This resembles results for iron-based superconductors but has a fundamental difference with iron pnictides and selenides. Moreover, our present efforts also suggest La3Ni2O7 is qualitatively different from infinite-layer nickelates and cuprate superconductors.

Negative Refraction of Weyl Phonons at Twin Quartz Interfaces

In Nature, α-quartz crystals frequently form contact twins, which are two adjacent crystals with the same chemical structure but different crystallographic orientation, sharing a common lattice plane. As α-quartz crystallizes in a chiral space group, such twinning can occur between enantiomorphs with the same handedness or with opposite handedness. Here, we use first-principles methods to investigate the effect of twinning and chirality on the bulk and surface phonon spectra, as well as on the topological properties of phonons in α-quartz. We demonstrate that, even though the dispersion appears identical for all twins along all high-symmetry lines and at all high-symmetry points in the Brillouin zone, the dispersions can be distinct at generic momenta for some twin structures. Furthermore, when the twinning occurs between different enantiomorphs, the charges of all Weyl nodal points flip, which leads to mirror symmetric isofrequency contours of the surface arcs on certain surfaces. We show that this allows negative refraction to occur at interfaces between certain twins of α-quartz.

Machine Learning K-Means Clustering in Interpolative Separable Density Fitting Algorithm: Advancing Accurate and Efficient Cubic-Scaling Density Functional Perturbation Theory Calculations within Plane Waves

Density functional perturbation theory (DFPT) is a crucial tool for accurately describing lattice dynamics. The adaptively compressed polarizability (ACP) method reduces the computational complexity of DFPT calculations from O(N4) to O(N3) by combining the interpolative separable density fitting (ISDF) algorithm. However, the conventional QR factorization with column pivoting (QRCP) algorithm, used for selecting the interpolation points in ISDF, not only incurs a high cubic-scaling computational cost, O(N3), but also leads to suboptimal convergence. This convergence issue is particularly pronounced when considering the complex interplay between the external potential and atomic displacement in ACP-based DFPT calculations. Here, we present a machine learning K-means clustering algorithm to select the interpolation points in ISDF, which offers a more efficient quadratic-scaling O(N2) alternative to the computationally intensive cubic-scaling O(N3) QRCP algorithm. We implement this efficient K-means-based ISDF algorithm to accelerate plane-wave DFPT calculations in KSSOLV, which is a MATLAB toolbox for performing Kohn-Sham density functional theory calculations within plane waves. We demonstrate that this K-means algorithm not only offers comparable accuracy to QRCP in ISDF but also achieves better convergence for ACP-based DFPT calculations. In particular, K-means can remarkably reduce the computational cost of selecting the interpolation points by nearly 2 orders of magnitude compared to QRCP in ISDF.

Basis-Set-Error-Free Random-Phase Approximation Correlation Energies for Atoms Based on the Sternheimer Equation

The finite basis set errors for all-electron random-phase approximation (RPA) correlation energy calculations are analyzed for isolated atomic systems. We show that, within the resolution-of-identity (RI) RPA framework, the major source of the basis set errors is the incompleteness of the single-particle atomic orbitals used to expand the Kohn-Sham eigenstates, instead of the auxiliary basis set (ABS) to represent the density response function χ0 and the bare Coulomb operator v. By solving the Sternheimer equation for the first-order wave function on a dense radial grid, we are able to eliminate the major error─the incompleteness error of the single-particle atomic basis set─for atomic RPA calculations. The error stemming from a finite ABS can be readily rendered vanishingly small by increasing the size of the ABS, or by iteratively determining the eigenmodes of the χ0v operator. The variational property of the RI-RPA correlation energy can be further exploited to optimize the ABS in order to achieve fast convergence of the RI-RPA correlation energy. These numerical techniques enable us to obtain basis-set-error-free RPA correlation energies for atoms, and in this work, such energies for atoms from H to Kr are presented. The implications of the numerical techniques developed in the present work for addressing the basis set issue for molecules and solids are discussed.

Calculation of phonons in real-space density functional theory

We present an accurate and efficient formulation for the calculation of phonons in real-space Kohn-Sham density functional theory. Specifically, employing a local exchange-correlation functional, norm-conserving pseudopotential in the Kleinman-Bylander representation, and local form for the electrostatics, we derive expressions for the dynamical matrix and associated Sternheimer equation that are particularly amenable to the real-space finite-difference method, within the framework of density functional perturbation theory. In particular, the formulation is applicable to insulating as well as metallic systems of any dimensionality, enabling the efficient and accurate treatment of semi-infinite and bulk systems alike, for both orthogonal and nonorthogonal cells. We also develop an implementation of the proposed formulation within the high-order finite-difference method. Through representative examples, we verify the accuracy of the computed phonon dispersion curves and density of states, demonstrating excellent agreement with established plane-wave results.

Anisotropic Mechanical Properties of Orthorhombic SiP2 Monolayer: A First-Principles Study

In recent years, the two-dimensional (2D) orthorhombic SiP2 flake has been peeled off successfully by micromechanical exfoliation and it exhibits an excellent performance in photodetection. In this paper, we investigated the mechanical properties and the origin of its anisotropy in an orthorhombic SiP2 monolayer through first-principles calculations, which can provide a theoretical basis for utilizing and tailoring the physical properties of a 2D orthorhombic SiP2 in the future. We found that the Young's modulus is up to 113.36 N/m along the a direction, while the smallest value is only 17.46 N/m in the b direction. The in-plane anisotropic ratio is calculated as 6.49, while a similar anisotropic ratio (~6.55) can also be observed in Poisson's ratio. Meanwhile, the in-plane anisotropic ratio for the fracture stress of the orthorhombic SiP2 monolayer is up to 9.2. These in-plane anisotropic ratios are much larger than in black phosphorus, ReS2, and biphenylene. To explain the origin of strong in-plane anisotropy, the interatomic force constants were obtained using the finite-displacement method. It was found that the maximum of interatomic force constant along the a direction is 5.79 times of that in the b direction, which should be considered as the main origin of the in-plane anisotropy in the orthorhombic SiP2 monolayer. In addition, we also found some negative Poisson's ratios in certain specific orientations, allowing the orthorhombic SiP2 monolayer to be applied in next-generation nanomechanics and nanoelectronics.

Natural Optical Activity from Density-Functional Perturbation Theory

We present an accurate and computationally efficient first-principles methodology to calculate natural optical activity. Our approach is based on the long-wave density-functional perturbation theory and includes self-consistent field terms naturally in the formalism, which are found to be of crucial importance. The final result is expressed exclusively in terms of response functions to uniform field perturbations and avoids troublesome summations over empty states. Our strategy is validated by computing the natural optical activity tensor in representative chiral crystals (trigonal Se, α-HgS, and α-SiO_{2}) and molecules (C_{4}H_{4}O_{2}), finding excellent agreement with experiment and previous theoretical calculations.

Discovery of enhanced lattice dynamics in a single-layered hybrid perovskite

Layered hybrid perovskites exhibit emergent physical properties and exceptional functional performances, but the coexistence of lattice order and structural disorder severely hinders our understanding of these materials. One unsolved problem regards how the lattice dynamics are affected by the dimensional engineering of the inorganic frameworks and their interaction with the molecular moieties. Here, we address this question by using a combination of spontaneous Raman scattering, terahertz spectroscopy, and molecular dynamics simulations. This approach reveals the structural dynamics in and out of equilibrium and provides unexpected observables that differentiate single- and double-layered perovskites. While no distinct vibrational coherence is observed in double-layered perovskites, an off-resonant terahertz pulse can drive a long-lived coherent phonon mode in the single-layered system. This difference highlights the dramatic change in the lattice environment as the dimension is reduced, and the findings pave the way for ultrafast structural engineering and high-speed optical modulators based on layered perovskites.

Opening a Band Gap in Biphenylene Monolayer via Strain: A First-Principles Study

A biphenylene network is a novel 2D allotropy of carbon with periodic 4-6-8 rings, which was synthesized successfully in 2021. In recent years, although the mechanical properties and thermal transport received a lot of research attention, how to open the Dirac cone in the band structure of a biphenylene network is still a confused question. In this work, we utilized uniaxial and biaxial lattice strains to manipulate the electronic properties and phonon frequencies of biphenylene, and we found an indirect band gap under 10% biaxial strain through the first-principles calculations. This indirect band gap is caused by the competition between the band-edge state A and the Dirac cone for the conduction band minimum (CBM). Additionally, the lightest carrier's effective mass in biphenylene is 0.184 m₀ for electrons along x (Γ→X) direction, while the effective mass for holes shows a remarkable anisotropy, suggesting the holes in the tensile biphenylene monolayer are confined within a one-dimensional chain along x direction. For phonon dispersion, we discovered that the Raman-active Ag3 phonon mode shows a robust single phonon mode character under both compressive and tensile strain, but its frequency is sensitive to lattice strain, suggesting the lattice strain in biphenylene can be identified by Raman spectroscopy.

Vibrational properties of metastable orthorhombic Bi2Se3

Chalcogenide materials are being considered as some of the most promising systems for energy harvesting and energy conversion. Among them, the orthorhombic family of compounds X₂Y₃ (with X = Bi, Sb and Y = S, Se) has attracted special attention due to its interesting atomic structure and thermoelectric and optical properties. While Bi₂S₃ and Sb₂Se₃ have already been applied to solar cells, practical application of the new metastable Bi₂Se₃ is still a challenge due to the lack of data and knowledge on its properties. Here, the vibrational and structural properties of the orthorhombic metastable phase of Bi₂Se₃ are investigated by using Raman spectroscopy and ab initio calculations. We perform Raman spectroscopy measurements along with in situ thermal treatment on samples grown by electrochemical deposition. We show that by increasing the temperature an improved crystallization occurs in the orthorhombic structure, followed by recrystallization to the usual rhombohedral phase. Our results point out specific Raman modes of the orthorhombic phase. First principles computational results based on the density functional theory support the experimental data and describe three singlet Raman active vibrational modes, such as B(2)3g, B(2)2g and A(6)g.

Enhancing the Curie Temperature in Cr2Ge2Te6 via Charge Doping: A First-Principles Study

In this work, we explore the impacts of charge doping on the magnetism of a Cr₂Ge₂Te₆ monolayer using first-principles calculations. Our results reveal that doping with 0.3 electrons per unit cell can enhance the ferromagnetic exchange constant in a Cr₂Ge₂Te₆ monolayer from 6.874 meV to 10.202 meV, which is accompanied by an increase in the Curie temperature from ~85 K to ~123 K. The enhanced ratio of the Curie temperature is up to 44.96%, even higher than that caused by surface functionalization on monolayer Cr₂Ge₂Te₆, manifesting the effectiveness of charge doping by improving the magnetic stability of 2D magnets. This remarkable enhancement in the ferromagnetic exchange constant and Curie temperature can be attributed to the increase in the magnetic moment on the Te atom, enlarged Cr-Te-Cr bond angle, reduced Cr-Te distance, and the significant increase in super-exchange coupling between Cr and Te atoms. These results demonstrate that charge doping is a promising route to improve the magnetic stability of 2D magnets, which is beneficial to overcome the obstacles in the application of 2D magnets in spintronics.

[K2 PbX][Ga7 S12 ] (X = Cl, Br, I): The First Lead-Containing Cationic Moieties with Ultrahigh Second-Harmonic Generation and Band Gaps Exceeding the Criterion of 2.33 eV

In contrast to anionic group theory of nonlinear optical (NLO) materials that second-harmonic generation (SHG) responses mainly originate from anionic groups, structural regulation on the cationic groups of salt-inclusion chalcogenides (SICs) is performed to make them also contribute to the NLO effects. Herein, the stereochemically active lone-electron-pair Pb2+ cation is first introduced to the cationic groups of NLO SICs, and the resultant [K2 PbX][Ga7 S12 ] (X = Cl, Br, I) are isolated via solid-state method. The features of their three-dimensional structures comprise highly oriented [Ga7 S12 ]3- and [K2 PbX]3+ frameworks derived from AgGaS2 , which display the largest phase-matching SHG intensities (2.5-2.7 × AgGaS2 @1800 nm) among all SICs. Concurrently, three compounds manifest band gap values of 2.54, 2.49, and 2.41 eV (exceeding the criterion of 2.33 eV), which can avoid two-photon absorption under the fundamental laser of 1064 nm, along with the relatively low anisotropy of thermal expansion coefficients, leading to improved laser-induced damage thresholds (LIDTs) values of 2.3, 3.8, and 4.0 times that of AgGaS2 . In addition, the density of states and SHG coefficient calculations demonstrate that the Pb2+ cations narrow the band gaps and benefit SHG responses.

Vibrational circular dichroism spectra of natural products by means of the nuclear velocity perturbation theory

We present the application of the recently implemented nuclear velocity perturbation theory, using the combined Gaussian and plane waves approach in CP2K, to the vibrational circular dichroism (VCD) spectra of a set of natural products. Even though the calculations were carried out for isolated molecules in the gas-phase limit, neglecting inter-molecular interactions and anharmonic effects, the match between simulated and experimental spectra is reasonable. We also study the influence of different density functionals on the conformational search and the resulting VCD spectra via group coupling matrices (GCMs). The GCM analysis reveals that the VCD signal can in some cases arise from moieties which are close to each other and in other cases from moieties far from each other. Differences in spectra obtained using different exchange-correlation density functionals can be attributed to interaction terms between different moieties in the molecules changing their sign.

Insights into selected 2D piezo Rashba semiconductors for self-powered flexible piezo spintronics: material to contact properties

The new paradigm in electronics consists in realizing the seamless integration of many properties latent in nanomaterials, such as mechanical flexibility, strong spin-orbit coupling (Rashba spin splitting-RSS), and piezoelectricity. Taking cues from the pointers given on 1D ZnO nanowires (ACS Nano2018121811-20), the concept can be extended to multifunctional two-dimensional (2D) materials, which can serve as an ideal platform in next-generation electronics such as self-powered flexible piezo-spintronic device. However, a microscopically clear understanding reachable from the state-of-the-art density functional theory-based approaches is a prerequisite to advancing this research domain. Atomic-scale insights gained from meticulously performed scientific computations can firmly anchor the growth of this important research field, and that is of undeniable relevance from scientific and technological outlooks. This article reviews the scientific advance in understanding 2D materials hosting all the essential properties, i.e. flexibility, piezoelectricity, and RSS. Important 2D semiconducting monolayers that deserve a special mention, include monolayers of buckled MgX (X = S, Se, Te), CdTe, ZnTe, Janus structures of transition metal trichalcogenides, Janus tellurene and 2D perovskites. van Der Waals multilayers are also built to design multifunctional materials via modulation of the stacking sequence and interlayer coupling between the constituent layers. External electric field, strain engineering and charge doping are perturbations mainly used to tune the spintronic properties. Finally, the contact properties of these monolayers are also crucial for their actual implementation in electronic devices. The nature of the contacts, Schottky/Ohmic, needs to be carefully examined first as it controls the device's performance. In this regard, the rare occurrence of Ohmic contact in graphene/MgS van der Waals hetero bilayer has been presented in this review article.

A variational formulation of the Harris functional as a correction to approximate Kohn-Sham density functional theory

Accurate descriptions of intermolecular interactions are of great importance in simulations of molecular liquids. We present an electronic structure method that combines the accuracy of the Harris functional approach with the computational efficiency of approximately linear-scaling density functional theory (DFT). The non-variational nature of the Harris functional has been addressed by constructing a Lagrangian energy functional, which restores the variational condition by imposing stationarity with respect to the reference density. The associated linear response equations may be treated with linear-scaling efficiency in an atomic orbital based scheme. Key ingredients to describe the structural and dynamical properties of molecular systems are the forces acting on the atoms and the stress tensor. These first-order derivatives of the Harris Lagrangian have been derived and implemented in consistence with the energy correction. The proposed method allows for simulations with accuracies close to the Kohn-Sham DFT reference. Embedded in the CP2K program package, the method is designed to enable ab initio molecular dynamics simulations of molecular solutions for system sizes of several thousand atoms. Available subsystem DFT methods may be used to provide the reference density required for the energy correction at near linear-scaling efficiency. As an example of production applications, we applied the method to molecular dynamics simulations of the binary mixtures cyclohexane-methanol and toluene-methanol, performed within the isobaric-isothermal ensemble, to investigate the hydrogen bonding network in these non-ideal mixtures.

Reversible Transformation and Distribution Determination of Diverse Pt Single-Atom Species

In single-atom catalysts (SACs), the complexity of the support anchoring sites creates a vast diversity of single-atom species with varied coordination environments. To date, the quantitative distribution of these diverse single-atom species in a given SAC has remained elusive. Recently, CeO₂-supported metal SACs have been extensively studied by modulating their local environments via numerous synthetic strategies. However, owing to the absence of a quantitative description, unraveling the site-specific reactivity and regulating their transformation remain challenging. Here, we show that two distinct Pt/CeO₂ SACs can be reversibly generated by oxidative and nonoxidative dispersions, which contain varied Pt₁O_n-Ce^δ+ single-atom species despite similar Pt charge states and coordination numbers. By means of Raman spectroscopy and computational studies, we semiquantitatively reveal the distribution of diverse Pt₁O_n-Ce^δ+ species in each specific SACs. Remarkably, the minority species of Pt₁O₄-Ce³⁺-O_v accounting for only 14.2% affords the highest site-specific reactivity for low-temperature CO oxidation among the other abundant counterparts, i.e., Pt₁O₄-Ce⁴⁺ and Pt₁O₆-Ce⁴⁺. The second nearest oxygen vacancy (O_v) not only acts synergistically with the nearby active metal sites to lower the reaction barrier but also facilitates the dynamic transformation from six-coordinated to four-coordinated sites during cyclic nonoxidative and oxidative dispersions. This work elucidates the quantitative distribution and dynamic transformation of varied single-atom species in a given SAC, offering a more intrinsic descriptor and quantitative measure to depict the inhomogeneity of SACs.

Stability and Strength of Monolayer Polymeric C60

Two-dimensional fullerene networks have been synthesized in several forms, and it is unknown which monolayer form is stable under ambient conditions. Using first-principles calculations, I show that the believed stability of the quasi-tetragonal phases is challenged by mechanical, dynamic, or thermodynamic stability. For all temperatures, the quasi-hexagonal phase is thermodynamically the least stable. However, the relatively high dynamic and mechanical stabilities suggest that the quasi-hexagonal phase is intrinsically stronger than the other phases under various strains. The origin of the high stability and strength of the quasi-hexagonal phase can be attributed to the strong covalent C-C bonds that strongly hold the linked C₆₀ clusters together, enabling the closely packed hexagonal network. These results rationalize the experimental observations that so far only the quasi-hexagonal phase has been exfoliated experimentally as monolayers.

K2CdGe3S8: A New Infrared Nonlinear Optical Sulfide

A quaternary metal chalcogenide, namely K2CdGe3S8 (I), is obtained through a high-temperature solid-state approach. Compound I crystallizes with the non-centrosymmetric space group P212121. It features a 2D layer structure with [CdGe3S8] layers consisting of tetrahedral GeS4 and CdS4 units, and counter K+ embedded between the layers. The compound exhibits a powder second-harmonic generation (SHG) response of ~0.1 times that of KH2PO4 (KDP) with phase-matchable character at the laser wavelength of 1064 nm. Remarkably, it has a wide band gap (3.20 eV), which corresponds to a favorable high laser-induced damage threshold of 6.7 times that of AgGaS2. In addition, the calculated birefringence (Δn) is 0.039 at the wavelength of 1064 nm, which satisfies the Δn criteria for a promising infrared NLO material.

Theoretical study on ferroelectric nitrides with super-wurtzite structures for solar energy conversion applications

Polarized structured nitride semiconductors are attractive due to their unique and environment-friendly electronic properties. The stability, ferroelectricity and photocatalytic and photovoltaic properties of super-wurtzite Mg₂XN₃ (X = Bi, Mo, Nb, Sb, Ta, Tc and W) were determined based on first principles calculations in this study. The calculated results indicate that Mg₂XN₃ (X = Sb, Ta, Bi and Nb) are stable polar nitrides by phonon frequencies, elastic coefficients and ferroelectric analysis. Mg₂XN₃ (X = Sb, Ta and Nb) with large ferroelectric polarization strength could absorb ultraviolet light to promote photocatalytic water splitting for hydrogen production. Mg₂BiN₃ is a new excellent photovoltaic candidate due to its ideal energy band, high electron mobility, high absorption coefficient and large ferroelectric polarization strength.

Single-Crystal Diffraction, Raman Spectroscopy, and Density Functional Theory of DTO [N-(1,7-Dinitro-1,2,6,7-tetrahydro-[1,3,5]triazino[1,2-c][1,3,5]oxadiazin-8(4H)-ylidene)nitramide]

A combined experimental and modeling study of energetic compound N-(1,7-dinitro-1,2,6,7-tetrahydro-[1,3,5]triazino[1,2-c][1,3,5]oxadiazin-8(4H)-ylidene)nitramide [C₅H₆N₈O₇, (DTO)] has been performed. We report its crystal structure, solid-phase heat of formation, and its vibrational and electronic structure obtained by single-crystal X-ray diffractometry, Raman spectroscopy, and density functional theory (DFT). DTO exhibits two adjoining six-membered rings, a triazine ring (C₃N₃) and an oxadiazine ring (C₃N₂O) ring containing two nitro functional groups and one nitroamino group. DTO crystallizes with four molecules in its unit cell and presents a density of 1.862 kg/m³ at 298 K, in excellent agreement with both DFT calculations performed both at the molecular level using the B3LYP with the 6-311+G** basis set and the solid-state level using the hybrid functional HSE6 optimized with norm-conserving pseudopotentials. The calculated vibrational structure allows for the symmetry assignment of key Raman modes in terms of atomic movements, and the calculated frequency values are in good agreement with experiment. The solid-phase DFT calculations reveal that the N atoms of the triazine ring contribute mostly to the density of states at the Fermi level. In addition, we present and discuss the computed solid-phase heat of formation (215.9 kJ/mol) and molecular electrostatic potential surface of DTO and compare them to complementary materials.

First-principles study of two-dimensional C-silicyne nanosheet as a promising anode material for rechargeable Li-ion batteries

Li-ion batteries are one of the sustainable alternatives to meet the growing energy demands of an increasing population. However, finding a suitable negative electrode is key for improving battery performance. In the present work, first principles-based investigations are carried out to explore the capability of a planar 2D C-silicyne nanosheet - which is a Si analogue of α-graphyne having -CC- substitution - as an anode for improving the performance of Li-ion batteries. Thermally and dynamically stable C-silicyne sheets exhibit a metallic nature as inferred from the density of states studies. The average adsorption energies for sequential adsorption of the Li atom over the monolayer range from -1.35 to -0.46 eV, implying favourable interactions between the monolayer and the Li atom which indicate that during the lithiation process, clustering amongst the metal atoms is not preferred. The energy barrier for the migration of Li-ions is 0.21 eV, indicating an active charge/discharge process. A high storage capacity of 836.07 mA h g^-1 and a working potential of 0.60 V is obtained. A negligible amount of volume change of the C-silicyne monolayer after full lithiation is observed which implies good cyclability. All these outcomes imply that C-silicyne nanosheets are a potential anode material for next-generation LIBs.

Comparison of Length, Velocity, and Symmetric Gauges for the Calculation of Absorption and Electric Circular Dichroism Spectra with Real-Time Time-Dependent Density Functional Theory

A velocity and symmetric gauge implementation for real-time time-dependent density functional theory (RT-TDDFT) in the CP2K package using a Gaussian and plane wave approach is presented, including the explicit gauge-transformed contributions due to the nonlocal part of pseudopotentials. Absorption spectra of gas-phase α-pinene are calculated in length and velocity gauges in the long-wavelength approximation for the application of a δ pulse in linear and full order. The velocity gauge implementation is also applied to a solvated uracil molecule to showcase its use within periodic boundary conditions (PBC). For the calculation of the expectation value of the electric dipole moment in PBC, both the velocity representation and the modern theory of polarization give equivalent absorption spectra if a distributed reference point is used for the nonlocal term of the velocity operator. The discussion of linear response theory takes place in a unified framework in terms of linear response functions in propagator notation, distinguishing the parts of the linear response functions associated with perturbation and response. To further investigate gauge dependence, electric circular dichroism (ECD) spectra of α-pinene were calculated either as magnetic response to an electric field perturbation, in length or velocity gauge, or as electric response to a magnetic field perturbation in the symmetric gauge. Both approaches, electric and magnetic perturbations, have been found to yield equivalent ECD spectra.

Molecular Dipole Moment Learning via Rotationally Equivariant Gaussian Process Regression with Derivatives in Molecular-orbital-based Machine Learning

This study extends the accurate and transferable molecular-orbital-based machine learning (MOB-ML) approach to modeling the contribution of electron correlation to dipole moments at the cost of Hartree--Fock computations. A molecular-orbital-based (MOB) pairwise decomposition of the correlation part of the dipole moment is applied, and these pair dipole moments could be further regressed as a universal function of molecular orbitals (MOs).The dipole MOB features consist of the energy MOB features and their responses to electric fields. An interpretable and rotationally equivariant Gaussian process regression (GPR) with derivatives algorithm is introduced to learn the dipole moment more efficiently. The proposed problem setup, feature design, and ML algorithm are shown to provide highly-accurate models for both dipole moment and energies on water and fourteen small molecules. To demonstrate the ability of MOB-ML to function as generalized density-matrix functionals for molecular dipole moments and energies of organic molecules, we further apply the proposed MOB-ML approach to train and test the molecules from the QM9 dataset. The application of local scalable GPR with Gaussian mixture model unsupervised clustering (GMM/GPR) scales up MOB-ML to a large-data regime while retaining the prediction accuracy. In addition, compared with literature results, MOB-ML provides the best test MAEs of 4.21 mDebye and 0.045 kcal/mol for dipole moment and energy models, respectively, when training on 110000 QM9 molecules. The excellent transferability of the resulting QM9 models is also illustrated by the accurate predictions for four different series of peptides.

Using Feature-Assisted Machine Learning Algorithms to Boost Polarity in Lead-Free Multicomponent Niobate Alloys for High-Performance Ferroelectrics

To expand the unchartered materials space of lead-free ferroelectrics for smart digital technologies, tuning their compositional complexity via multicomponent alloying allows access to enhanced polar properties. The role of isovalent A-site in binary potassium niobate alloys, (K,A)NbO3 using first-principles calculations is investigated. Specifically, various alloy compositions of (K,A)NbO3 are considered and their mixing thermodynamics and associated polar properties are examined. To establish structure-property design rules for high-performance ferroelectrics, the sure independence screening sparsifying operator (SISSO) method is employed to extract key features to explain the A-site driven polarization in (K,A)NbO3 . Using a new metric of agreement via feature-assisted regression and classification, the SISSO model is further extended to predict A-site driven polarization in multicomponent systems as a function of alloy composition, reducing the prediction errors to less than 1%. With the machine learning model outlined in this work, a polarity-composition map is established to aid the development of new multicomponent lead-free polar oxides which can offer up to 25% boosting in A-site driven polarization and achieving more than 150% of the total polarization in pristine KNbO3 . This study offers a design-based rational route to develop lead-free multicomponent ferroelectric oxides for niche information technologies.

Direct Observation of Coherent Longitudinal and Shear Acoustic Phonons in TaAs Using Ultrafast X-Ray Diffraction

Using femtosecond time-resolved x-ray diffraction, we investigated optically excited coherent acoustic phonons in the Weyl semimetal TaAs. The low symmetry of the (112) surface probed in our experiment enables the simultaneous excitation of longitudinal and shear acoustic modes, whose dispersion closely matches our simulations. We observed an asymmetry in the spectral line shape of the longitudinal mode that is notably absent from the shear mode, suggesting a time-dependent frequency chirp that is likely driven by photoinduced carrier diffusion. We argue on the basis of symmetry that these acoustic deformations can transiently alter the electronic structure near the Weyl points and support this with model calculations. Our study underscores the benefit of using off-axis crystal orientations when optically exciting acoustic deformations in topological semimetals, allowing one to transiently change their crystal and electronic structures.

Implementation of Nuclear Velocity Perturbation and Magnetic Field Perturbation Theory in CP2K and Their Application to Vibrational Circular Dichroism

We present the implementation of nuclear velocity perturbation theory (NVPT), using a pioneering combination of atom-centered (velocity-dependent) Gaussian basis functions and plane waves in the CP2K package. The atomic polar tensors (APTs) and atomic axial tensors (AATs) are evaluated in the velocity representation using efficient density functional perturbation theory. The presence of nonlocal pseudopotentials, the representation of potentials on numerical integration grids, and effects arising from the basis functions being centered on the atoms have been considered in the implementation. The Magnetic Field Perturbation Theory (MFPT) using gauge-including atomic orbitals is implemented in the same code and compared to the NVPT. Our implementation is the first to compare both approaches (MFPT and NVPT) in the same code. The implementation has been verified via sum rules and by investigating the gauge origin dependence of the AATs for a set of small molecules, oxirane, and fluoro-oxirane. We also present vibrational circular dichroism spectra that are related to the APTs and AATs, applying both theories.

A Concise Review on Recent Developments of Machine Learning for the Prediction of Vibrational Spectra

Machine learning has become more and more popular in computational chemistry, as well as in the important field of spectroscopy. In this concise review, we walk the reader through a short summary of machine learning algorithms and a comprehensive discussion on the connection between machine learning methods and vibrational spectroscopy, particularly for the case of infrared and Raman spectroscopy. We also briefly discuss state-of-the-art molecular representations which serve as meaningful inputs for machine learning to predict vibrational spectra. In addition, this review provides an overview of the transferability and best practices of machine learning in the prediction of vibrational spectra as well as possible future research directions.

High-Field Magnetoelectric and Spin-Phonon Coupling in Multiferroic (NH4)2[FeCl5·(H2O)]

We combine high field polarization, magneto-infrared spectroscopy, and lattice dynamics calculations with prior magnetization to explore the properties of (NH₄)₂[FeCl₅·(H₂O)]─a type II molecular multiferroic in which the mixing between charge, structure, and magnetism is controlled by intermolecular hydrogen and halogen bonds. Electric polarization is sensitive to the series of field-induced spin reorientations, increasing linearly with the field and reaching a maximum before collapsing to zero across the quasi-collinear to collinear-sinusoidal reorientation due to the restoration of inversion symmetry. Magnetoelectric coupling is on the order of 1.2 ps/m for the P∥c, H∥c configuration between 5 and 25 T at 1.5 K. In this range, the coupling takes place via an orbital hybridization mechanism. Other forms of mixing are active in (NH₄)₂[FeCl₅·(H₂O)] as well. Magneto-infrared spectroscopy reveals that all of the vibrational modes below 600 cm^-1 are sensitive to the field-induced transition to the fully saturated magnetic state at 30 T. We analyze these local lattice distortions and use frequency shifts to extract spin-phonon coupling constants for the Fe-O stretch, Fe-OH₂ rock, and NH₄⁺ libration. Inspection also reveals subtle symmetry breaking of the ammonium counterions across the ferroelectric transition. The coexistence of such varied mixing processes in a platform with intermolecular hydrogen- and halogen-bonding opens the door to greater understanding of multiferroics and magnetoelectrics governed by through-space interactions.

Li-free Ternary Sulphide Cs5Ga9S16 with Excellent Nonlinear Optical Performance Similar to Classic LiGaS2

The element with the strongest electropositivity, Cs, was introduced into an anionic framework [Ga9S16]5− constructed from GaS4 tetrahedrons, which are functional nonlinear optical motifs, to acquire the ternary noncentrosymmetric sulphide...

Spectroscopic Signatures of Hydrogen-Bonding Motifs in Protonic Ionic Liquid Systems: Insights from Diethylammonium Nitrate in the Solid State

Diethylammonium nitrate, [N0 0 2 2][NO3], and its perdeuterated analogue, [N D D 2 2] [NO3], were structurally characterized and studied by infrared, Raman, and inelastic neutron scattering (INS) spectroscopy. Using these experimental data along with state-of-the-art computational materials modeling, we report unambiguous spectroscopic signatures of hydrogen-bonding interactions between the two counterions. An exhaustive assignment of the spectral features observed with each technique has been provided, and a number of distinct modes related to NH···O dynamics have been identified. We put a particular emphasis on a detailed interpretation of the high-resolution, broadband INS experiments. In particular, the INS data highlight the importance of conformational degrees of freedom within the alkyl chains, a ubiquitous feature of ionic liquid (IL) systems. These findings also enable an in-depth physicochemical understanding of protonic IL systems, a first and necessary step to the tailoring of hydrogen-bonding networks in this important class of materials.

Developing a simple box-behnken experimental design on the removal of doxorubicin anticancer drug using Fe3O4/graphene nanoribbons adsorbent

This paper aims to develop a Box-Behnken experimental design system to optimize the removal process of doxorubicin anticancer drugs. For this goal, Fe₃O₄/graphene nanoribbons was selected as adsorbent and removal of doxorubicin anticancer drug optimized using Box-Behnken experimental design with a selection of four effective factors. A three-level, four-factor Box-Behnken experimental design was used to assess the relationship between removal percentage as a dependent variable with adsorption weight (0.0015-0.01 mg), pH (3-9), temperature (15-45 °C) and time (1-15 min) as independent variables. Optimized condition by Behnken experimental design (pH = 7.36; time = 15 min; adsorbent weight = 0.01 mg and temperature = 29.26 °C) improved removal of doxorubicin anticancer drug about 99.2% in aqueous solution. The dynamic behavior, adsorption properties and mechanism of doxorubicin molecule on Fe₃O₄/graphene nanoribbon were investigated based on ab initio molecular dynamics (AIMD) simulations and density functional theory calculations with dispersion corrections. A closer inspection of the adsorption configurations and binding energies revealed that π-π interactions were the driving force when the doxorubicin molecule adsorbed on Fe₃O₄/graphene nanoribbon. The observed negative adsorption energy signifies a favourable and exothermic adsorption process of the various adsorbate-substrate systems. Besides, AIMD and phonon dispersion calculations confirm the dynamic stability of Fe₃O₄/graphene nanoribbon.

Magnetic order-dependent phonon properties in 2D magnet CrI3

We carried out a systematic theoretical study on how spin affects the phononic properties of CrI3 monolayers. We find that the frequencies of two infrared-active (IR) modes are significantly influenced by the magnetic configuration. Thus an IR spectrum may be applied to identify the magnetic order by utilizing the spin-lattice correlation. The thermal expansion coefficients are 2.21, 3.35 and -5.58 × 10-6 K-1 for ferromagnetic (FM), antiferromagnetic (AFM) and paramagnetic (PM) phases at 30 K, because of the competition between the modes with negative and positive Grüneisen constants. Furthermore, the lattice thermal conductivity is also sensitive to the magnetic phase, which is attributed to the spin-dependent lattice anharmonicity. Our results provide fundamental insights into the spin-lattice coupling and clarify the potential of a spintronic monolayer as a thermal switching device for active heat flow control.

Dynamics & Spectroscopy with Neutrons-Recent Developments & Emerging Opportunities

This work provides an up-to-date overview of recent developments in neutron spectroscopic techniques and associated computational tools to interrogate the structural properties and dynamical behavior of complex and disordered materials, with a focus on those of a soft and polymeric nature. These have and continue to pave the way for new scientific opportunities simply thought unthinkable not so long ago, and have particularly benefited from advances in high-resolution, broadband techniques spanning energy transfers from the meV to the eV. Topical areas include the identification and robust assignment of low-energy modes underpinning functionality in soft solids and supramolecular frameworks, or the quantification in the laboratory of hitherto unexplored nuclear quantum effects dictating thermodynamic properties. In addition to novel classes of materials, we also discuss recent discoveries around water and its phase diagram, which continue to surprise us. All throughout, emphasis is placed on linking these ongoing and exciting experimental and computational developments to specific scientific questions in the context of the discovery of new materials for sustainable technologies.