Search

The response to compression of the clathrate type-II structure Ge(cF136) is investigated by means of \textit{ab initio} small-cell metadynamics at different temperatures and pressures. At lower pressure $p$=2.5 GPa the metastable metallic bct-5 phase competes against $\beta$-Sn Ge(tI4), which forms at higher pressures. On lowering temperature, the presence of amorphous intermediates obtained from Ge$_{136}$ is more pronounced and is instrumental to the formation of denser structural motifs, from which metallic bct-5 can form. Therein, anisotropic box fluctuations promote phase formation. The metadynamics runs are analysed in depth using a set of topological descriptors, including coordination sequence and ring statistics. Differences in the structural landscape history and amorphous intermediates are critical for the selective formation of particular metastable polymorphs, towards turning crystal structure predictions into actual materials.

TiO2 nanoparticles (NPs) are attracting materials for biomedical applications, provided that they are coated with polymers to improve solubility, dispersion and biocompatibility. Conformation, coverage density and solvent effects largely influence their functionality and stability. In this work, we use atomistic molecular dynamics simulations to study polyethylene glycol (PEG) grafting to highly curved TiO2 NPs (2-3 nm) in different solvents. We compare the coating polymer conformations on NPs with those on (101) flat surfaces. In water, the transition from mushroom to brush conformation starts only at high density ({\sigma} = 2.25 chains/nm2). In dichloromethane (DCM), at low-medium coverage ({\sigma} < 1.35 chains/nm2), several interactions between the PEG chains backbone and undercoordinated Ti atoms are established, whereas at {\sigma} = 2.25 chains/nm2 the conformation clearly becomes brush-like. Finally, we demonstrate that these spherical brushes, when immersed in water, but not in DCM, follow the Daoud-Cotton (DC) classical scaling model for the polymer volume fraction dependence with the distance from the center of star-shaped systems.

Metal oxide nanoparticles (NPs) are regarded as good candidates for many technological applications, where their functional environment is often an aqueous solution. The correct description of metal oxide electronic structure is still a challenge for local and semilocal density functionals, whereas hybrid functional methods provide an improved description, and local atomic function based codes such as CRYSTAL17 outperform plane wave codes when it comes to hybrid functional calculations. However, the computational cost of hybrids are still prohibitive for systems of real sizes, in a real environment. Therefore, we here present and critically assess the accuracy of our electrostatic embedding quantum mechanical/molecular mechanical (QM/MM) coupling between CRYSTAL17 and AMBER16, and demonstrate some of its capabilities via the case study of TiO2 NPs in water. First, we produced new Lennard-Jones (LJ) parameters that improve the accuracy of water-water interactions in the B3LYP/TIP3P coupling. Then, we applied our QM/MM coupling methodology to describe the interaction of a 1 nm wide multilayer of water surrounding a spherical TiO2 nanoparticle (NP). Optimizing the QM/MM water water parameters was found to have little to no effect on the local NP properties, which provide insights into the range of influence that can be attributed to the LJ term in the QM/MM coupling. The effect of adding additional water in an MM fashion on the geometry optimized nanoparticle structure is small, but more evident effects are seen in its electronic properties. We also show that there is good transferability of existing QM/MM LJ parameters for organic molecules water interactions to our QM/MM implementation, even though these parameters were obtained with a different QM code and QM/MM implementation, but with the same functional.

Computational modeling of titanium dioxide nanoparticles of realistic size is extremely relevant for the direct comparison with experiments but it is also a rather demanding task. We have recently worked on a multistep/scale procedure to obtain global optimized minimum structures for chemically stable spherical titania nanoparticles of increasing size, with diameter from 1.5 nm (~300 atoms) to 4.4 nm (~4000 atoms). We use first self-consistent-charge density functional tight-binding (SCC-DFTB) methodology to perform thermal annealing simulations to obtain globally optimized structures and then hybrid density functional theory (DFT) to refine them and to achieve high accuracy in the description of structural and electronic properties. This allows also to assess SCC-DFTB performance in comparison with DFT(B3LYP) results. As a further step, we investigate photoexcitation and photoemission processes involving electron/hole pair formation, separation, trapping and recombination in the nanosphere of medium size by hybrid DFT. Finally, we show how a recently defined new set of parameters for SCC-DFTB allows for a proper description of titania/water multilayers interface, which paves the way for modeling large realistic nanoparticles in aqueous environment.

TiO2 nanoparticles (NPs) are nowadays considered fundamental building blocks for many technological applications. Morphology is found to play a key role with spherical NPs presenting higher binding properties and chemical activity. From the experimental point of view, the characterization of these nano-objects is extremely complex, opening a large room for computational investigations. In this work, TiO2 spherical NPs of different size (from 300 to 4000 atoms) have been studied with a two-scale computational approach. Global optimization to obtain stable and equilibrated NSs was performed with a self-consistent charge density functional tight-binding (SCC-DFTB) simulated annealing process, causing a considerable atomic rearrangement within the nanospheres. Those SCC-DFTB relaxed structures have been then optimized at DFT(B3LYP) level of theory. We present a systematic and comparative SCC-DFTB vs DFT(B3LYP) study of the structural properties, with particular emphasis on the surface-to-bulk sites ratio, coordination distribution of surface sites and surface energy. From the electronic point of view, we compare HOMO-LUMO and Kohn-Sham gaps, total and projected density of states. Overall, the comparisons between DFTB and hybrid DFT show that DFTB provides a rather accurate geometrical and electronic description of these nanospheres of realistic size (up to a diameter of 4.4 nm) at an extremely reduced computational cost. This opens for new challenges in simulations of very large systems and more extended molecular dynamics.

A detailed understanding of the relation between microscopic structure and phonon propagation at the nan oscale is essential to design materials with desired phononic and thermal properties.Here we uncover a new mechanism of phonon interaction in surface oxidized membranes, i.e., native oxide layers interact with phonons in ultra-thin silicon membranes through local resonances. The local resonances reduce the low frequency phonon group velocities and shorten their mean free path. This effect opens up a new strategy for ultralow thermal conductivity design as it complements the scattering mechanism which scatters higher frequency modes effectively. The combination of native oxide layer and alloying with germanium in concentration as small as 5% reduces the thermal conductivity of silicon membranes to 100 time lower than the bulk. In addition, the resonance mechanism produced by native oxide surface layers is particularly effective for thermal condutivity reduction even at very low temperatures, at which only low frequency modes are populated., Comment: 6 pages, 5 figures, Accepted for publication in Physical Review B

We use ab initio density functional calculations to determine the interaction of a graphene monolayer with the Si(111) surface. We found that graphene forms strong bonds to the bare substrate and accommodates the 12% lattice mismatch by forming a wavy structure consisting of free-standing conductive ridges that are connected by ribbon-shaped regions of graphene, which bond covalently to the substrate. We perform quantum transport calculations for different geometries to study changes in the transport properties of graphene introduced by the wavy structure and bonding to the Si substrate. Our results suggest that wavy graphene combines high mobility along the ridges with efficient carrier injection into Si in the contact regions.

By means of ab initio metadynamics runs we explored the lower-pressure region of the phase diagram of germanium. A monoclinic germanium phase with four-membered rings, less dense than diamond and compressible into \beta-tin phase (tI4) was found. A metallic bct-5 phase, mechanically stable down to room conditions appeared between diamond and tI4. mC16 is a narrow-gap semiconductor, while bct-5 is metallic and potentially still superconducting in the very low pressure range. This finding may help resolving outstanding experimental issues., Comment: 6 figures

Four sp3 carbon allotropes with six, eight, and 16 atoms per primitive cell have been derived using a combination of metadynamics simulations and topological scan. A chiral orthorhombic phase oC16 (C2221) was found to be harder than monoclinic M-carbon and shows excellent stability in the high-pressure range. A second orthorhombic phase of Cmmm symmetry, by \sim 0.028 eV/atom energetically lower than W-carbon, can be formed from graphite at \sim 9 GPa. In general, the mechanical response under pressure was found to depend on the structure topology, which reflects the way rings are formed from an initial graphene layer stacking., Comment: 5 pages, 5 figures