Abstract: Using first-principles calculations, we study the effect of embedding various atoms into the nanopore sites of both C6N6 and C6N8 monolayers. Our results indicate that the embedded atoms significantly affect the electronic and magnetic properties of C6N6 and C6N8 monolayers and lead to extraordinary and multifarious electronic properties, such as metallic, half-metallic, spin-glass semiconductor and dilute-magnetic semiconductor behaviour. Our results reveal that the H atom concentration dramatically affects the C6N6 monolayer. On increasing the H coverage, the impurity states also increase due to H atoms around the Fermi-level. C6N6 shows metallic character when the H atom concentration reaches 6.25%. Moreover, the effect of charge on the electronic properties of both Cr@C6N6 and C@C6N8 is also studied. Cr@C6N6 is a ferromagnetic metal with a magnetic moment of 2.40 mu(B), and when 0.2 electrons are added and removed, it remains a ferromagnetic metal with a magnetic moment of 2.57 and 2.77 mu(B), respectively. Interestingly, one can observe a semi-metal, in which the VBM and CBM in both spin channels touch each other near the Fermi-level. C@C6N8 is a semiconductor with a nontrivial band gap. When 0.2 electrons are removed, it remains metallic, and under excess electronic charge, it exhibits half-metallic behaviour.

Abstract: In this paper, the effect of BC3, C3N and C3N4BC(3) and MoS2/C(3)N4 heterostructures are direct semiconductors with band gaps of 0.4 and 1.74 eV, respectively, while MoS2/C3N is a metal. Furthermore, the influence of strain and electric field on the electronic structure of these van der Waals heterostructures is investigated. The MoS2/BC3 heterostructure, for strains larger than -4%, transforms it into a metal where the metallic character is maintained for strains larger than -6%. The band gap decreases with increasing strain to 0.35 eV (at +2%), while for strain (>+6%) a direct-indirect band gap transition is predicted to occur. For the MoS2/C3N heterostructure, the metallic character persists for all strains considered. On applying an electric field, the electronic properties of MoS2/C3N4 are modified and its band gap decreases as the electric field increases. Interestingly, the band gap reaches 30 meV at +0.8 V/angstrom, and with increase above +0.8 V/angstrom, a semiconductor-to-metal transition occurs. Furthermore, we investigated effects of semi- and full-hydrogenation of MoS2/C3N and we found that it leads to a metallic and semiconducting character, respectively.

Abstract: Recently, a two-dimensional (2D) MoSi 2N 4 (MSN) structure has been successfully synthesized [Hong et al., Science 369(6504), 670-674 (2020)]. Motivated by this result, we investigate the structural, electronic, and optical properties of MSN monolayer (MSN-1L) and bilayer (MSN-2L) under the applied electric field (E-field) and strain using density functional theory calculations. We find that the MSN-2L is a semiconductor with an indirect bandgap of 1.60 (1.80)eV using Perdew-Burke-Ernzerhof (HSE06). The bandgap of MSN-2L decreases as the E-field increases from 0.1 to 0.6V/angstrom and for larger E-field up to 1.0V/angstrom the bilayer becomes metallic. As the vertical strain increases, the bandgap decreases; more interestingly, a semiconductor to a metal phase transition is observed at a strain of 12 %. Furthermore, the optical response of the MSN-2L is in the ultraviolet (UV) region of the electromagnetic spectrum. The absorption edge exhibits a blue shift by applying an E-field or a vertical compressive strain. The obtained interesting properties suggest MSN-2L as a promising material in electro-mechanical and UV opto-mechanical devices.

Abstract: First‐principles calculations, within the framework of density functional theory, are used to investigate the structural, electronic, optical, and thermoelectric properties of monolayer PbTe. The effect of layer thickness, electric field, strain, and vacancy defects on the electronic and magnetic properties is systematically studied. The results show that the bandgap decreases as the layer thickness increases from monolayer to bulk. With application of an electric field on bilayer PbTe, the bandgap decreases from 70 meV (0.2 V Å⁻¹) to 50 meV (1 V Å⁻¹) when including spin–orbit coupling (SOC). Application of uniaxial strain induces a direct‐to‐indirect bandgap transition for strain greater than +6%. In addition, the bandgap decreases under compressive biaxial strain (with SOC). The effect of vacancy defects on the electronic properties of PbTe is also investigated. Such vacancy defects turn PbTe into a ferromagnetic metal (single vacancy Pb) with a magnetic moment of 1.3 μB, and into an indirect semiconductor with bandgap of 1.2 eV (single Te vacancy) and 1.5 eV (double Pb + Te vacancy). In addition, with change of the Te vacancy concentration, a bandgap of 0.38 eV (5.55%), 0.43 eV (8.33%), and 0.46 eV (11.11%) is predicted.

Abstract: Herein, using first-principles calculations the structural and electronic properties of the Ti(2)CO(2)MXene monolayer with and without oxygen vacancies are systematically investigated with different defect concentrations and patterns, including partial, linear, local, and hexagonal types. The Ti(2)CO(2)monolayer is found to be a semiconductor with a bandgap of 0.35 eV. The introduction of oxygen vacancies tends to increase the bandgap and leads to electronic phase transitions from nonmagnetic semiconductors to half-metals. Moreover, the semiconducting characteristic of O-vacancy Ti(2)CO(2)can be adjusted via electric fields, strain, and F-atom substitution. In particular, an electric field can be used to alter the nonmagnetic semiconductor of O-vacancy Ti(2)CO(2)into a magnetic one or into a half-metal, whereas the electronic phase transition from a semiconductor to metal can be achieved by applying strain and F-atom substitution. The results provide a useful guide for practical applications of O-vacancy Ti(2)CO(2)monolayers in nanoelectronic and spinstronic nanodevices.

Abstract: In this work, we extensively investigate the structural and electronic properties of van der Waals heterostructures (HTs) constructed by MoS${2}$/$BC₃$, MoS${2}$/$C₃N$, MoS${2}$/$C₃N₄$, MoS${2}$/$C₄N₃$ and those using Janus MoSSe instead of MoS$2$ by performing density functional theory calculations. The electronic band structure calculations and the corresponding partial density of states reveal that the significant changes are driven by quite strong layer-layer interaction between the constitutive layers. Our results show that although all monolayers are semiconductors as free-standing layers, the MoS${2}$/$C₃N$ and MoS${2}$/$C₄N₃$ bilayer HTs display metallic behavior as a consequence of transfer of charge carriers between two constituent layers. In addition, it is found that in MoSSe/$C₃N$ bilayer HT, the degree of metallicity is affected by the interface chalcogen atom type when Se atoms are facing to $C₃N$ layer, the overlap of the bands around the Fermi level is smaller. Moreover, the half-metallic magnetic $C₄N₃$ is shown to form magnetic half-metallic trilayer HT with MoS$2$ independent of the stacking sequence, i.e. whether it is sandwiched or two $C₄N₃$ layer encapsulate MoS$2$ layer. We further analyze the trilayer HTs in which MoS$2$ is encapsulated by two different monolayers and it is revealed that at least with one magnetic monolayer, it is possible to construct a magnetic trilayer. While the trilayer of $C₄N₃$/MoS${2}$/$BC₃$ and $C₄N₃$/MoS${2}$/$C₃N₄$ exhibit half-metallic characteristics, $C₄N₃$/MoS${_2}$/$C₃$N possesses a magnetic metallic ground state. Overall, our results reveal that holly structures of BCN crystals are suitable for heterostructure formation even over van der Waals type interaction which significantly changes electronic nature of the constituent layers.

Abstract: Latest synthesis of ZnSb monolayer, encouraged us to conduct density functional theory (DFT) simulations in order to study the structural, magnetic, electronic/optical and mechanical features of the sp2-hybridized honeycomb ZnSb monolayer (ML-ZnSb) and bilayer (BL-ZnSb). Our structural optimizations reveal that ML-ZnSb is an anisotropic hexagonal structure while BL-ZnSb is composed of shifted ZnSb layers which are covalently binded. ML-ZnSb is found to be a ferromagnetic metal, in contrast BL-ZnSb has a non-magnetic indirect band gap semiconducting ground state. For the in-plane polarization, first absorption peak of ML-ZnSb and BL-ZnSb confirm the absorbance of the light within the infrared domain wand visible range, respectively. Moreover, our results reveal that the layer-layer chemical bonding in BL-ZnSb significantly enhances the mechanical response of ML-ZnSb whose in-plane stiness is the smallest among all 2D materials (2DM). Notably, the strong in-plane anisotropy of ML-ZnSb in its stiness reduces in BL-ZnSb.

Abstract: The crystallization of nanomaterials is a primary source of solid-state, photonic structures. Thus, a detailed understanding of this process is of paramount importance for the successful application of photonic nanomaterials in emerging optoelectronic technologies. While colloidal crystallization has been thoroughly studied, for example, with advanced in situ electron microscopy methods, the noncolloidal crystallization (freezing) of nanoparticles (NPs) remains so far unexplored. To fill this gap, in this work, we present proof-of principle experiments decoding a crystallization of reconfigurable assemblies of NPs at a solid state. The chosen material corresponds to an excellent testing bed, as it enables both in situ and ex situ investigation using X-ray diffraction (XRD), transmission electron microscopy (TEM), high-angle annular dark-field scanning transmission electron microscopy (HAADF-STEM), atomic force microscopy (AFM), and optical spectroscopy in visible and ultraviolet range (UV−vis) techniques. In particular, ensemble measurements with small-angle XRD highlighted the dependence of the correlation length in the NPs assemblies on the number of heating/cooling cycles and the rate of cooling. Ex situ TEM imaging further supported these results by revealing a dependence of domain size and structure on the sample preparation route and by showing we can control the domain size over 2 orders of magnitude. The application of HAADF-STEM tomography, combined with in situ thermal control, provided three-dimensional single-particle level information on the positional order evolution within assemblies. This combination of real and reciprocal space provides insightful information on the anisotropic, reversibly reconfigurable assemblies of NPs. TEM measurements also highlighted the importance of interfaces in the polydomain structure of nanoparticle solids, allowing us to understand experimentally observed differences in UV−vis extinction spectra of the differently prepared crystallites. Overall, the obtained results show that the combination of in situ heating HAADF-STEM tomography with XRD and ex situ TEM techniques is a powerful approach to study nanoparticle freezing processes and to reveal the crucial impact of disorder in the solid-state aggregates of NPs on their plasmonic properties.

Abstract: Triple-modal imaging magnetic nanocapsules, encapsulating hydrophobic superparamagnetic iron oxide nanoparticles, are formulated and used to magnetically target solid tumours after intravenous administration in tumour-bearing mice. The engineered magnetic polymeric nanocapsules m-NCs are ~200 nm in size with negative Zeta potential and shown to be spherical in shape. The loading efficiency of superparamagnetic iron oxide nanoparticles in the m-NC was ~100%. Up to ~3- and ~2.2-fold increase in tumour uptake at 1 and 24 h was achieved, when a static magnetic field was applied to the tumour for 1 hour. m-NCs, with multiple imaging probes (e.g. indocyanine green, superparamagnetic iron oxide nanoparticles and indium-111), were capable of triple-modal imaging (fluorescence/magnetic resonance/nuclear imaging) in vivo. Using triple-modal imaging is to overcome the intrinsic limitations of single modality imaging and provides complementary information on the spatial distribution of the nanocarrier within the tumour. The significant findings of this study could open up new research perspectives in using novel magnetically-responsive nanomaterials in magnetic-drug targeting combined with multi-modal imaging.

Abstract: The evaluation of nucleation rates from molecular dynamics trajectories is hampered by the slow nucleation time scale and impact of finite size effects. Here, we show that accurate nucleation rates can be obtained in a very general fashion relying only on the free energy barrier, transition state theory, and a simple dynamical correction for diffusive recrossing. In this setup, the time scale problem is overcome by using enhanced sampling methods, in casu metadynamics, whereas the impact of finite size effects can be naturally circumvented by reconstructing the free energy surface from an appropriate ensemble. Approximations from classical nucleation theory are avoided. We demonstrate the accuracy of the approach by calculating macroscopic rates of droplet nucleation from argon vapor, spanning 16 orders of magnitude and in excellent agreement with literature results, all from simulations of very small (512 atom) systems.

Abstract: We present an approach to incorporate the effect of vibrational nonequilibrium in molecular dynamics (MD) simulations. A perturbed canonical ensemble, in which selected modes are excited to higher temperature while all others remain equilibrated at low temperature, is simulated by applying a specifically tailored bias potential. Our method can be readily applied to any (classical or quantum mechanical) MD setup at virtually no additional computational cost and allows the study of reactions of vibrationally excited molecules in nonequilibrium environments such as plasmas. In combination with enhanced sampling methods, the vibrational efficacy and mode selectivity of vibrationally stimulated reactions can then be quantified in terms of chemically relevant observables, such as reaction rates and apparent free energy barriers. We first validate our method for the prototypical hydrogen exchange reaction and then show how it can capture the effect of vibrational excitation on a symmetric SN2 reaction and radical addition on CO2.

Abstract: Atomistic simulation methods for the quantification of free energies are in wide use. These methods operate by sampling the probability density of a system along a small set of suitable collective variables (CVs), which is, in turn, expressed in the form of a free energy surface (FES). This definition of the FES can capture the relative stability of metastable states but not that of the transition state because the barrier height is not invariant to the choice of CVs. Free energy barriers therefore cannot be consistently computed from the FES. Here, we present a simple approach to calculate the gauge correction necessary to eliminate this inconsistency. Using our procedure, the standard FES as well as its gauge-corrected counterpart can be obtained by reweighing the same simulated trajectory at little additional cost. We apply the method to a number of systems—a particle solvated in a Lennard-Jones fluid, a Diels–Alder reaction, and crystallization of liquid sodium—to demonstrate its ability to produce consistent free energy barriers that correctly capture the kinetics of chemical or physical transformations, and discuss the additional demands it puts on the chosen CVs. Because the FES can be converged at relatively short (sub-ns) time scales, a free energy-based description of reaction kinetics is a particularly attractive option to study chemical processes at more expensive quantum mechanical levels of theory.

Abstract: Uniform-acceptance force-bias Monte Carlo (fbMC) methods have been shown to be a powerful technique to access longer timescales in atomistic simulations allowing, for example, phase transitions and growth. Recently, a new fbMC method, the time-stamped force-bias Monte Carlo (tfMC) method, was derived with inclusion of an estimated effective timescale; this timescale, however, does not seem able to explain some of the successes the method. In this contribution, we therefore explicitly quantify the effective timescale tfMC is able to access for a variety of systems, namely a simple single-particle, one-dimensional model system, the Lennard-Jones liquid, an adatom on the Cu(100) surface, a silicon crystal with point defects and a highly defected graphene sheet, in order to gain new insights into the mechanisms by which tfMC operates. It is found that considerable boosts, up to three orders of magnitude compared to molecular dynamics, can be achieved for solid state systems by lowering of the apparent activation barrier of occurring processes, while not requiring any system-specific input or modifications of the method. We furthermore address the pitfalls of using the method as a replacement or complement of molecular dynamics simulations, its ability to explicitly describe correct dynamics and reaction mechanisms, and the association of timescales to MC simulations in general.

Abstract: Atomistic simulations can in principle provide an unbiased description of all mechanisms, intermediates, and products of complex chemical processes. However, due to the severe time scale limitation of conventional simulation techniques, unrealistically high simulation temperatures are usually applied, which are a poor approximation of most practically relevant low-temperature applications. In this work, we demonstrate the direct observation at the atomic scale of the pyrolysis and oxidation of n-dodecane at temperatures as low as 700 K through the use of a novel simulation technique, collective variable-driven hyperdynamics (CVHD). A simulated timescale of up to 39 seconds is reached. Product compositions and dominant mechanisms are found to be strongly temperature-dependent, and are consistent with experiments and kinetic models. These simulations provide a first atomic-level look at the full dynamics of the complicated fuel combustion process at industrially relevant temperatures and time scales, unattainable by conventional molecular dynamics simulations.

Abstract: Optimization of catalytic materials for a given application is greatly constrained by linear scaling relations. Recently, however, it has been demonstrated that it is possible to reversibly modulate the chemisorption of molecules on nanomaterials by charging (i.e., injection or removal of electrons) and hence reversibly and selectively modify catalytic activity beyond structure−activity correlations. The fundamental physical relation between the properties of the material, the charging process, and the chemisorption energy, however, remains unclear, and a systematic exploration and optimization of charge-switchable sorbent materials is not yet possible. Using hybrid DFT calculations of CO2 chemisorption on hexagonal boron nitride nanosheets with several types of defects and dopants, we here reveal the existence of fundamental correlations between the electron affinity of a material and charge-induced chemisorption, show how defect engineering can be used to modulate the strength and efficiency of the adsorption process, and demonstrate that excess electrons stabilize many topological defects. We then show how these insights could be exploited in the development of new electrocatalytic materials and the synthesis of doped nanomaterials. Moreover, we demonstrate that calculated chemical properties of charged materials are highly sensitive to the employed computational methodology because of the self-interaction error, which underlines the theoretical challenge posed by such systems.