Abstract: While spectroscopic data on small hydrocarbons in interstellar media in combination with crossed molecular beam (CMB) experiments have provided a wealth of information on astrochemically relevant species, much of the underlying mechanistic pathways of their formation remain elusive. Therefore, in this work, the chemical reaction mechanisms of C(³P_J) + C₆H₆and C⁺(²P) + C₆H₆systems using the quantum mechanical molecular dynamics (QMMD) technique at the PBE0-D3(BJ) level of theory is investigated, mimicking a CMB experiment. Both the dynamics of the reactions as well as the electronic structure for the purpose of the reaction network are evaluated. The method is validated for the first reaction by comparison to the available experimental data. The reaction scheme for the C(³P_J) + C₆H₆system covers the literature data,<italic>e.g.</italic>the major products are the 1,2-didehydrocycloheptatrienyl radical (C₇H₅) and benzocyclopropenyl radical (C₆H₅–CH), and it reveals the existence of less common pathways for the first time. The chemistry of the C⁺(²P_J) + C₆H₆system is found to be much richer, and we have found that this is because of more exothermic reactions in this system in comparison to those in the C(³P_J) + C₆H₆system. Moreover, using the QMMD simulation, a number of reaction paths have been revealed that produce three distinct classes of reaction products with different ring sizes. All in all, at all the collision energies and orientations, the major product is the heptagon molecular ion for the ionic system. It is also revealed that the collision orientation has a dominant effect on the reaction products in both systems, while the collision energy mostly affects the charged system. These simulations both prove the applicability of this approach to simulate crossed molecular beams, and provide fundamental information on reactions relevant for the interstellar medium.

Abstract: We study the structure and phase behavior of fullerene-cubane C60·C8H8 by Monte Carlo simulation. Using a simple potential model capturing the icosahedral and cubic symmetries of its molecular constituents, we reproduce the experimentally observed phase transition from a cubic to an orthorhombic crystal lattice and the accompanying rotational freezing of the C60 molecules. We elaborate a scheme to identify the low-temperature orientations of individual molecules and to detect a pattern of orientational ordering similar to the arrangement of C60 molecules in solid C60. Our configuration of orientations supports a doubled periodicity along one of the crystal axes.

Abstract: We report on a first-principles study of the structural and electronic properties of a Sigma3 (112) grain boundary model in CuInSe2. The study focuses on a coherent, stoichiometry preserving, cation–Se terminated grain boundary, addressing the properties of the grain boundary as such, as well as the effect
of well known defects in CuInSe2. We show that in spite of its apparent simplicity, such a grain boundary exhibits a very rich phenomenology, providing an explanation for several of the experimentally observed properties of grain boundaries in CuInSe2 thin films. In particular, we show that the combined effect of Cu vacancies and cation antisites can result in the observed Cu depletion with no In enrichment at the grain boundaries. Furthermore, Cu vacancies are unlikely to produce a hole barrier at the grain boundaries, but Na may indeed have such an effect. We find that Na-on-Cu defects will tend to form abundantly at
the grain boundaries, and can provide a mechanism for the carrier depletion and/or type inversion experimentally reported.

Abstract: We present a high-throughput computing scheme based on density functional theory (DFT) to generate a class of oxides and screen them with the aim of identifying those that might be electronically appropriate for transparent conducting oxide (TCO) applications. The screening criteria used are a minimum band gap to ensure sufficient transparency, a band edge alignment consistent with easy n- or p-type dopability, and a minimum thermodynamic phase stability to be experimentally synthesizable. Following this scheme we screened 23 binary and 1518 ternary bixbyite oxides in order to identify promising candidates, which can then be a subject of an in-depth study. The results for the known TCOs are in good agreement with the reported data in the literature. We suggest a list of several new potential TCOs, including both n- and p-type compounds.

Abstract: We have prepared tungsten oxide films decorated with gold particles on Si substrates by aerosol assisted chemical vapor deposition (AACVD) and characterized them using scanning electron microscopy (SEM), transmission electron microscopy (TEM) and X-ray photoelectron spectroscopy (XPS). SEM shows that the films are composed of needle-like structures and TEM shows that both the needles and the gold particles are crystalline. XPS indicates the presence of oxygen vacancies, i.e. the films are WO3−x, and hence the deposited material is composed of semiconducting nanostructures and that the interaction between the gold particles and the WO3 needles surface is weak. The synthesis of semiconducting tungsten oxide nanostructures decorated with metal particles represents an important step towards the development of sensing devices with optimal properties.

Abstract: We establish the topological conditions underlying the thermal stability of C30Si30 clusters. Two topologies have been considered: a segregated one, where Si and C atoms lie on neighboring and yet, separated parts of the cage, and a non-segregated one, where the number of SiC bonds is maximized. The segregated network is energetically favored against the non-segregated one, both structures being fully relaxed at T = 0 K. Conversely, the non-segregated structure is the only one stable at finite temperatures, regardless of the nature of the local states (d or p) included in the KleynmanBylander construction.

Abstract: We calculate bubble nucleation rates in a Lennard-Jones fluid through explicit molecular dynamics simulations. Our approach-based on a recent free energy method (dubbed reweighted Jarzynski sampling), transition state theory, and a simple recrossing correction-allows us to probe a fairly wide range of rates in several superheated and cavitation regimes in a consistent manner. Rate predictions from this approach bridge disparate independent literature studies on the same model system. As such, we find that rate predictions based on classical nucleation theory, direct brute force molecular dynamics simulations, and seeding are consistent with our approach and one another. Published rates derived from forward flux sampling simulations are, however, found to be outliers. This study serves two purposes: First, we validate the reliability of common modeling techniques and extrapolation approaches on a paradigmatic problem in materials science and chemical physics. Second, we further test our highly generic recipe for rate calculations, and establish its applicability to nucleation processes.

Abstract: Using the simplified linear combination of atomic orbitals (LCAO) method in combination with ab initio calculations, we construct a tight-binding (TB) model for two different crystal structures of monolayer gallium: a(100)- and b(010)-Gallenene. The analytical expression for the Hamiltonian and numerical results for the overlap matrix elements between different orbitals of the Ga atoms and for the Slater and Koster (SK) integrals are obtained. We find that the compaction of different structures affects significantly the formation of the orbitals. The results for a(100)-Gallenene can be very well explained with an orthogonal basis set, while for b(010)-Gallenene we have to assume a non-orthogonal basis set in order to construct the TB model. Moreover, the transmission properties of nanoribbons of both monolayers oriented along the AC and ZZ directions are also investigated and it is shown that both AC- and ZZ-b(010)-Gallenene nanoribbons exhibit semiconducting behavior with zero transmission while those of a(100)-Gallenene nanoribbons are metallic.

Abstract: Using first-principles many-body perturbation theory, we investigate the optical properties of 8-Pmmn borophene at two levels of approximations; the GW method considering only the electron-electron interaction and the GW in combination with the Bethe-Salpeter equation including electron-hole coupling. The band structure exhibits anisotropic Dirac cones with semimetallic character. The optical absorption spectra are obtained for different light polarizations and we predict strong optical absorbance anisotropy. The absorption peaks undergo a global redshift when the electron-hole interaction is taken into account due to the formation of bound excitons which have an anisotropic excitonic wave function.

Abstract: Using first-principles calculations we investigate the interaction of various molecules, including H-2, N-2, CO, CO2, H2O, H2S, NH3, CH4 with a C3N nanosheet. Due to the weaker interaction between H-2, N-2, CO, CO2, H2O, H2S, NH3, and CH4 molecules with C3N, the adsorption energy is small and does not yield any significant distortion of the C3N lattice and the molecules are physisorbed. Calculated charge transfer shows that these molecules act as weak donors. However, adsorption of O-2, NO, NO2 and SO2 molecules are chemisorbed, they receive electrons from C3N and act as a strong acceptor. They interact strongly through hybridizing its frontier orbitals with the p-orbital of C3N, modifying the electronic structure of C3N. Our theoretical studies indicate that C3N-based sensor has a high potential for O-2, NO, NO2 and SO2 molecules detection.

Abstract: Using first principles calculations, we study the structural and phononic properties of the recently synthesized Janus type single layers of molybdenum dichalcogenides. The Janus MoSSe single layer possesses 2H crystal structure with two different chalcogenide sides that lead to out-of-plane anisotropy. By virtue of the asymmetric structure of the ultra-thin Janus type crystal, we induced the out-of-plane anisotropy to show the distinctive vertical pressure effect on the vibrational properties of the Janus material. It is proposed that for the corresponding Raman active optical mode of the Janus structure, the phase modulation and the magnitude ratio of the strained atom and its first neighbor atom adjust the distinctive change in the eigen-frequencies and Raman activity. Moreover, a strong variation in the Raman activity of the Janus structure is obtained under bivertical and univertical strains. Not only eigen-frequency shifts but also Raman activities of the optical modes of the Janus structure exhibit distinguishable features. This study reveals that the vertical anisotropic feature of the Janus structure under Raman measurement allows us to distinguish which side of the Janus crystal interacts with the externals (substrate, functional adlayers, or dopants). Published by AIP Publishing.

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: Two-dimensional Transition Metal Dichalcogenides (TMDs) possessing extraordinary physical properties at reduced dimensionality have attracted interest due to their promise in electronic and optical device applications. However, TMD monolayers can show a broad range of different properties depending on their crystal phase; for example, H phases are usually semiconductors, while the T phases are metallic. Thus, controlling phase transitions has become critical for device applications. In this study, the energetically low-lying crystal structures of pristine and Janus TMDs are investigated by using ab initio Nudged Elastic Band and molecular dynamics simulations to provide a general explanation for their phase stability and transition properties. Across all materials investigated, the T phase is found to be the least stable and the H phase is the most stable except for WTe2, while the T' and T '' phases change places according to the TMD material. The transition energy barriers are found to be large enough to hint that even the higher energy phases are unlikely to undergo a phase transition to a more stable phase if they can be achieved except for the least stable T phase, which has zero barrier towards the T ' phase. Indeed, in molecular dynamics simulations the thermodynamically least stable T phase transformed into the T ' phase spontaneously while in general no other phase transition was observed up to 2100 K for the other three phases. Thus, the examined T ', T '' and H phases were shown to be mostly stable and do not readily transform into another phase. Furthermore, so-called mixed phase calculations considered in our study explain the experimentally observed lateral hybrid structures and point out that the coexistence of different phases is strongly stable against phase transitions. Indeed, stable complex structures such as metal-semiconductor-metal architectures, which have immense potential to be used in future device applications, are also possible based on our investigation.

Abstract: Two-dimensional (2D) materials have attracted growing interest due to their potential use in the next generation of nanoelectronic and optoelectronic applications. On the basis of first-principles calculations based on density functional theory, we first investigate the electronic and mechanical properties of single layer boron phosphide (h-BP). Our calculations show that h-BP is a mechanically stable 2D material with a direct band gap of 0.9 eV at the K-point, promising for both electronic and optoelectronic applications. We next investigate the electron transport properties of a p-n junction constructed from single layer boron phosphide (h-BP) using the non-equilibrium Green's function formalism. The n-and p-type doping of BP are achieved by substitutional doping of B with C and P with Si, respectively. C(Si) substitutional doping creates donor (acceptor) states close to the conduction (valence) band edge of BP, which are essential to construct an efficient p-n junction. By modifying the structure and doping concentration, it is possible to tune the electronic and transport properties of the p-n junction which exhibits not only diode characteristics with a large current rectification but also negative differential resistance (NDR). The degree of NDR can be easily tuned via device engineering.

Abstract: The thermodynamic aspects of various 2D materials are explored using Density Functional Theory (DFT). Various metal chalcogenides (MX2, M = metal, chalcogen X = S, Se, Te) are investigated with respect to their interaction and stability under different ambient conditions met in the integration process of a transistor device. Their interaction with high- dielectrics is also addressed, in order to assess their possible integration in Complementary Metal Oxide Semiconductor (CMOS) field effect transistors. 2D materials show promise for high performance nanoelectronic devices, but the presence of defects (vacancies, grain boundaries,...) can significantly impact their electronic properties. To assess the impact of defects, their enthalpies of formation and their signature levels in the density of states have been studied. We find, consistently with literature reports, that chalcogen vacancies are the most likely source of defects. It is shown that while pristine 2D materials are in general stable whenever set in contact with different ambient atmospheres, the presence of defective sites affects the electronic properties of the 2D materials to varying degrees. We observe that all the 2D materials studied in the present work show strong reactivity towards radical oxygen plasma treatments while reactivity towards other common gas phase chemical such as O-2 and H2O and groups present at the high- surface varies significantly between species. While energy band-gaps, effective masses and contact resistivities are key criteria in selection of 2D materials for scaled CMOS and tunneling based devices, the phase and ambient stabilities might also play a very important role in the development of reliable nanoelectronic applications.