Abstract: Photocatalytic soot oxidation is studied on P25 TiO2 as an important model reaction for self-cleaning processes by means of electron paramagnetic resonance (EPR) and Fourier transform infrared (FTIR) spectroscopy. Contacting of carbon black with P25 leads on the one hand to a reduction of the local dioxygen concentration in the powder. On the other hand, the weakly adsorbed radicals on the carbon particles are likely to act as alternative traps for the photogenerated conduction-band electrons. We find furthermore that the presence of dioxygen and oxygen-related radicals is vital for the photocatalytic soot degradation. The complete oxidation of soot to CO2 is evidenced by in situ FTIR spectroscopy, no intermediate CO is detected during the photocatalytic process.

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: 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: Despite significant progress in single walled carbon nanotube (SWCNT) production by plasma enhanced chemical vapor deposition (PECVD), the growth mechanism in this method is not clearly understood. We employ reactive molecular dynamics simulations to investigate how plasma-based deposition allows growth at low temperature. We first investigate the SWCNT growth mechanism at low and high temperatures under conditions similar to thermal CVD and PECVD. We then show how ion bombardment during the nucleation stage increases the carbon solubility in the catalyst at low temperature. Finally, we demonstrate how moderate energy ions sputter amorphous carbon allowing for SWCNT growth at 500 K. (C) 2013 Elsevier B. V. All rights reserved.

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: Motivated by a recent experiment that reported the synthesis of a new 2D material nitrogenated holey graphene (C2N) [Mahmood et al., Nat. Commun., 2015, 6, 6486], the electronic, magnetic, and mechanical properties of nitrogenated (C2N), phosphorated (C2P) and arsenicated (C2As) monolayer holey graphene structures are investigated using first-principles calculations. Our total energy calculations indicate that, similar to the C2N monolayer, the formation of the other two holey structures are also energetically feasible. Calculated cohesive energies for each monolayer show a decreasing trend going from the C2N to C2As structure. Remarkably, all the holey monolayers considered are direct band gap semiconductors. Regarding the mechanical properties (in-plane stiffness and Poisson ratio), we find that C2N has the highest in-plane stiffness and the largest Poisson ratio among the three monolayers. In addition, our calculations reveal that for the C2N, C2P and C2As monolayers, creation of N and P defects changes the semiconducting behavior to a metallic ground state while the inclusion of double H impurities in all holey structures results in magnetic ground states. As an alternative to the experimentally synthesized C2N, C2P and C2As are mechanically stable and flexible semiconductors which are important for potential applications in optoelectronics.

Abstract: We have performed a first-principles study of the p- and n-type conductivity in CuIn1−xGaxSe2 due to native point defects, based on the HSE06 hybrid functional. Band alignment shows that the band gap becomes larger with x due to the increasing conduction band minimum, rendering it hard to establish n-type conductivity in CuGaSe2. From the defect formation energies, we find that In/GaCu is a shallow donor, while VCu, VIn/Ga and CuIn/Ga act as shallow acceptors. Using the total charge neutrality of ionized defects and intrinsic charge carriers to determine the Fermi level, we show that under In-rich growth conditions InCu causes strongly n-type conductivity in CuInSe2. Under increasingly In-poor growth conditions, the conductivity type in CuInSe2 alters to p-type and compensation of the acceptors by InCu reduces, as also observed in photoluminescence experiments. In CuGaSe2, the native acceptors pin the Fermi level far away from the conduction band minimum, thus inhibiting n-type conductivity. On the other hand, CuGaSe2 shows strong p-type conductivity under a wide range of Ga-poor growth conditions. Maximal p-type conductivity in CuIn1−xGaxSe2 is reached under In/Ga-poor growth conditions, in agreement with charge concentration measurements on samples with In/Ga-poor stoichiometry, and is primarily due to the dominant acceptor CuIn/Ga.

Abstract: A microscopic theory, which describes the orientational dynamics of C60 molecules in the face-centered cubic phase of C60-fullerite, is formulated or the case of a complex molecular structure. Interaction centers which comprise atoms, double bonds, and single bonds as molecular constituents contribute to the intermolecular potential. Orientation dependent physical properties are described in terms of symmetry-adapted rotator functions. It is found that a same set of rotator functions is sufficient even in the case of a complex molecular structure. Phase transition temperatures are discussed for various models of molecular structure.

Abstract: Synchrotron small angle X-ray scattering (SAXS) was used to characterize silicalite-1 zeolite crystallization from TEOS/TPAOH/water clear sol. SAXS patterns were recorded over a broad range of length scales, enabling the simultaneous monitoring of nanoparticles and crystals occurring at various stages of the synthesis. A simple two-population model accurately described the patterns. Nanoparticles were modeled by polydisperse coreshell spheres and crystals by monodisperse oblate ellipsoids. These models were consistent with TEM images. The SAXS results, in conjunction with in situ light scattering, showed that nucleation of crystals occurred in a short period of time. Crystals were uniform in size and shape and became increasingly anisotropic during growth. In the presence of nanoparticles, crystal growth was fast. During crystal growth, the number of nanoparticles decreased gradually but their size was constant. These observations suggested that the nanoparticles were growth units in an aggregative crystal growth mechanism. Crystals grown in the presence of nanoparticles developed a faceted habit and intergrowths. In the final stages of growth, nanoparticles were depleted. Concurrently, the crystal growth rate decreased significantly.

Abstract: The in-plane stiffness (C), Poisson's ratio (nu), Young's modulus and ultimate strength (sigma) along two different crystallographic orientations are calculated for the single layer crystals: MoS2, HfS2 and TiS3 in 1H, 1T and monoclinic phases. We find that MoS2 and HfS2 have isotropic in-plane stiffnesses of 124.24 N m(-1) and 79.86 N m(-1), respectively. While for TiS3 the in-plane stiffness is highly anisotropic due to its monoclinic structure, with C-x = 83.33 N m(-1) and C-y = 133.56 N m(-1) (x and y are parallel to its longer and shorter in-plane lattice vectors.). HfS2 which is in the 1T phase has the smallest anisotropy in its ultimate strength, whereas TiS3 in the monoclinic phase has the largest. Along the armchair direction MoS2 has the largest sigma of 23.48 GPa, whereas along y TiS3 has the largest sigma of 18.32 GPa. We have further analyzed the band gap response of these materials under uniaxial tensile strain, and find that they exhibit different behavior. Along both armchair and zigzag directions, the band gap of MoS2 (HfS2) decreases (increases) as strain increases, and the response is almost isotropic. For TiS3, the band gap decreases when strain is along x, while if strain is along y, the band gap increases first and then decreases beyond a threshold strain value. The different characteristics observed in these sulphides with different structures shed light on the relationship between the structure and properties, which is useful for applications in nanotechnology.

Abstract: When confined between graphene layers, water behaves differently from the bulk and exhibits unusual properties such as fast water flow and ordering into a crystal. The hydrogen-bonded network is affected by the limited space and by the characteristics of the confining walls. The presence of an extraordinary number of hydronium and hydroxide ions in narrow channels has the following effects: (i) they affect water permeation through the channel, (ii) they may interact with functional groups on the graphene oxide surface and on the edges, and (iii) they change the thermochemistry of water, which are fundamentally important to understand, especially when confined water is subjected to an external electric field. Here we study the physical properties of water when confined between two graphene sheets and containing hydronium and hydroxide. We found that: (i) there is a disruption in the solvation structure of the ions, which is also affected by the layered structure of confined water, (ii) hydronium and hydroxide occupy specific regions inside the nanochannel, with a prevalence of hydronium (hydroxide) ions at the edges (interior), and (iii) ions recombine more slowly in confined systems than in bulk water, with the recombination process depending on the channel height and commensurability between the size of the molecules and the nanochannel height – a decay of 20% (40%) in the number of ions in 8 ps is observed for a channel height of h = 7 angstrom (bulk water). Our work reveals distinctive properties of water confined in a nanocapillary in the presence of additional hydronium and hydroxide ions.

Abstract: The interface heat transfer of two layers induced by van der Waals (vdW) contacts is theoretically investigated, based on first-principles calculations at low temperatures. The results suggest that out-of-plane acoustic phonons with low frequencies dominate the interface thermal transport due to the vdW interaction. The interface thermal conductivity is proportional to the cubic of temperature at very low temperatures, but becomes linearly proportional to temperature as temperature increases. We show that manipulating the strain alters vdW coupling, leading to increased interfacial thermal conductivity at the interface. Our findings provide valuable insights into the interface heat transport in vdW heterostructures and support further design and optimization of electronic and optoelectronic nanodevices based on vdW contacts. The heat transfer induced by van der Waals contacts is dominated by ZA phonons. The interface thermal conductivity is proportional to the cubic of temperature, but becomes linearly proportional to temperature as temperature increases.