Abstract: Anatase nanosheets with exposed {001} facets
have gained increasing interest for photocatalytic applications. To
fully understand the structure-to-activity relation, combined
experimental and computational methods have been exploited.
Anatase nanosheets were prepared under hydrothermal conditions
in the presence of fluorine ions. High resolution scanning
transmission electron microscopy was used to fully characterize
the synthesized material, confirming the TiO2 nanosheet
morphology. Moreover, the surface structure and composition
of a single nanosheet could be determined by annular bright-field
scanning transmission electron microscopy (ABF-STEM) and
STEM electron energy loss spectroscopy (STEM-EELS). The photocatalytic activity was tested for the decomposition of organic
dyes rhodamine 6G and methyl orange and compared to a reference TiO2 anatase sample. The anatase nanosheets with exposed
{001} facets revealed a significantly lower photocatalytic activity compared to the reference. In order to understand the
mechanism for the catalytic performance, and to investigate the role of the presence of F−, light-induced electron paramagnetic
resonance (EPR) experiments were performed. The EPR results are in agreement with TEM, proving the presence of Ti3+
species close to the surface of the sample and allowing the analysis of the photoinduced formation of paramagnetic species.
Further, ab initio calculations of the anisotropic effective mass of electrons and electron holes in anatase show a very high effective
mass of electrons in the [001] direction, having a negative impact on the mobility of electrons toward the {001} surface and thus
the photocatalysis. Finally, motivated by the experimental results that indicate the presence of fluorine atoms at the surface, we
performed ab initio calculations to determine the position of the band edges in anatase slabs with different terminations of the
{001} surface. The presence of fluorine atoms near the surface is shown to strongly shift down the band edges, which indicates
another reason why it can be expected that the prepared samples with a large amount of {001} surface, but with fluorine atoms
near the surface, show only a low photocatalytic activity.

Abstract: Group II-VI semiconductor honeycomb monolayers have a noncentrosymmetric crystal structure and therefore are expected to be important for nano piezoelectric device applications. This motivated us to perform first principles calculations based on density functional theory to unveil the piezoelectric properties (i.e., piezoelectric stress (e(11)) and piezoelectric strain (d(11)) coefficients) of these monolayer materials with chemical formula MX (where M = Be, Mg, Ca, Sr, Ba, Zr, Cd and X = S, Se, Te). We found that these two-dimensional materials have peculiar piezoelectric properties with d(11) coefficients 1 order of magnitude larger than those of commercially utilized bulk materials. A clear trend in their piezoelectric properties emerges, which

Abstract: Device performance of Si|SiOx core-shell based nanowires critically depends on the exact control over the oxide thickness. Low-temperature plasma oxidation is a highly promising alternative to thermal oxidation allowing for improved control over the oxidation process, in particular for ultrasmall Si nanowires. We here elucidate the room temperature plasma oxidation mechanisms of ultrasmall Si nanowires using hybrid molecular dynamics / force-bias Monte Carlo simulations. We demonstrate how the oxidation and concurrent water formation mechanisms are a function of the oxidizing plasma species and we demonstrate how the resulting core-shell oxide thickness can be controlled through these species. A new mechanism of water formation is discussed in detail. The results provide a detailed atomic level explanation of the oxidation process of highly curved Si surfaces. These results point out a route toward plasma-based formation of ultrathin core-shell Si|SiOx nanowires at room temperature.

Abstract: CsPbI3 is a promising perovskite material for photovoltaic applications in its photoactive perovskite or black phase. However, the material degrades to a photovoltaically inactive or yellow phase at room temperature. Various mitigation strategies are currently being developed to increase the lifetime of the black phase, many of which rely on inducing strains in the material that hinder the black-to-yellow phase transition. Physical insight into how these strategies exactly induce strain as well as knowledge of the spatial extent over which these strains impact the material is crucial to optimize these approaches but is still lacking. Herein, we combine machine learning potential-based molecular dynamics simulations with our in silico strain engineering approach to accurately quantify strained large-scale atomic structures on a nanosecond time scale. To this end, we first model the strain fields introduced by atomic substitutions as they form the most elementary strain sources. We demonstrate that the magnitude of the induced strain fields decays exponentially with the distance from the strain source, following a decay rate that is largely independent of the specific substitution. Second, we show that the total strain field induced by multiple strain sources can be predicted to an excellent approximation by summing the strain fields of each individual source. Finally, through a case study, we illustrate how this additive character allows us to explain how complex strain fields, induced by spatially extended strain sources, can be predicted by adequately combining the strain fields caused by local strain sources. Hence, the strain additivity proposed here can be adopted to further our insight into the complex strain behavior in perovskites and to design strain from the atomic level onward to enhance their sought-after phase stability.

Abstract: Morphological characterization of nanoparticle assemblies and hybrid nanomaterials is critical in determining their structure-property relationships as well as in the development of structures with desired properties. Electron tomography has become a widely utilized technique for the three-dimensional characterization of nanoparticle assemblies. However, the extraction of quantitative morphological parameters from the reconstructed volume can be a complex and labor-intensive task. In this study, we aim to overcome this challenge by automating the volumetric segmentation process applied to three-dimensional reconstructions of nanoparticle assemblies. The key to enabling automated characterization is to assess the performance of different volumetric segmentation methods in accurately extracting predefined quantitative descriptors for morphological characterization. In our methodology, we compare the quantitative descriptors obtained through manual segmentation with those obtained through automated segmentation methods, to evaluate their accuracy and effectiveness. To show generality, our study focuses on the characterization of assemblies of CdSe/CdS quantum dots, gold nanospheres and CdSe/CdS encapsulated in polymeric micelles, and silica-coated gold nanorods decorated with both CdSe/CdS or PbS quantum dots. We use two unsupervised segmentation algorithms: the watershed transform and the spherical Hough transform. Our results demonstrate that the choice of automated segmentation method is crucial for accurately extracting the predefined quantitative descriptors. Specifically, the spherical Hough transform exhibits superior performance in accurately extracting quantitative descriptors, such as particle size and interparticle distance, thereby allowing for an objective, efficient, and reliable volumetric segmentation of complex nanoparticle assemblies.

Abstract: The synthesis of gold nanorods requires the presence of symmetry-breaking and shape-directing additives, among which bromide ions and quaternary ammonium surfactants have been reported as essential. As a result, hexadecyltrimethylammonium bromide (CTAB) has been selected as the most efficient surfactant to direct anisotropic growth. One of the difficulties arising from this selection is the low solubility of CTAB in water at room temperature, and therefore the seeded growth of gold nanorods is usually performed at 25 degrees C or above, which has restricted so far the analysis of kinetic effects derived from lower temperatures. We report a systematic study of the synthesis of gold nanorods from pentatwinned seeds using hexadecyltrimethylammonium chloride (CTAC) as the principal surfactant and a low concentration of bromide as shape-directing agent. Under these conditions, the synthesis can be performed at temperatures as low as 8 degrees C, and the corresponding kinetic effects can be studied, resulting in temperature-controlled aspect ratio tunability.

Abstract: Direct experimental measurement of thermal expansion coefficient without substrate effects is a challenging task for two-dimensional (2D) materials, and its accurate estimation with large-scale ab initio molecular dynamics is computationally very expensive. Machine learning-based interatomic potentials trained with ab initio data have been successfully used in molecular dynamics simulations to decrease the computational cost without compromising the accuracy. In this study, we investigated using Gaussian approximation potentials to reproduce the density functional theory-level accuracy for graphene within both lattice dynamical and molecular dynamical methods, and to extend their applicability to larger length and time scales. Two such potentials are considered, GAP17 and GAP20. GAP17, which was trained with pristine graphene structures, is found to give closer results to density functional theory calculations at different scales. Further vibrational and structural analyses verify that the same conclusions can be deduced with density functional theory level in terms of the reasoning of the thermal expansion behavior, and the negative thermal expansion behavior is associated with long-range out-of-plane phonon vibrations. Thus, it is argued that the enabled larger system sizes by machine learning potentials may even enhance the accuracy compared to small-size-limited ab initio molecular dynamics.

Abstract: The quantum electrostatic heterostructure (QEH) model allows for efficient computation of the dielectric screening properties of layered van der Waals (vdW)-bonded heterostructures in terms of the dielectric functions of the individual two-dimensional (2D) layers. Here, we extend the QEH model by including (1) contributions to the dielectric function from infrared active phonons in the 2D layers, (2) screening from homogeneous bulk substrates, and (3) intraband screening from free carriers in doped 2D semiconductor layers. We demonstrate the potential of the extended QEH model by calculating the dispersion of coupled phonons in multilayer stacks of hexagonal boron-nitride (hBN), the strong hybridization of plasmons and optical phonons in graphene/hBN heterostructures, the effect of substrate screening on the exciton series of monolayer MoS2, and the properties of hyperbolic plasmons in a doped phosphorene sheet. The new QEH code is distributed as a Python package with a simple command line interface and a comprehensive library of dielectric building blocks for the most common 2D materials, providing an efficient open platform for dielectric modeling of realistic vdW heterostructures.

Abstract: Plasma is gaining interest for CH4 conversion into higher hydrocarbons and H2. However, the performance in terms of conversion and selectivity toward different hydrocarbons is different for different plasma types, and the underlying mechanisms are not yet fully understood. Therefore, we study here these mechanisms in different plasma sources, by means of a chemical kinetics model. The model is first validated by comparing the calculated conversions and hydrocarbon/H2 selectivities with experimental results in these different plasma types and over a wide range of specific energy input (SEI) values. Our model predicts that vibrational−translational nonequilibrium is negligible in all CH4 plasmas investigated, and instead, thermal conversion is important. Higher gas temperatures also lead to a more selective production of unsaturated hydrocarbons (mainly C2H2) due to neutral dissociation of CH4 and subsequent dehydrogenation processes, while three-body recombination reactions into saturated hydrocarbons (mainly C2H6, but also higher hydrocarbons) are dominant in low temperature plasmas.

Abstract: The hyperdynamics method is a powerful tool to simulate slow processes at the atomic level. However, the construction of an optimal hyperdynamics potential is a task that is far from trivial. Here, we propose a generally applicable implementation of the hyperdynamics algorithm, borrowing two concepts from metadynamics. First, the use of a collective variable (CV) to represent the accelerated dynamics gives the method a very large flexibility and simplicity. Second, a metadynamics procedure can be used to construct a suitable history-dependent bias potential on-the-fly, effectively turning the algorithm into a self-learning accelerated molecular dynamics method. This collective variable-driven hyperdynamics (CVHD) method has a modular design: both the local system properties on which the bias is based, as well as the characteristics of the biasing method itself, can be chosen to match the needs of the considered system. As a result, system-specific details are abstracted from the biasing algorithm itself, making it extremely versatile and transparent. The method is tested on three model systems: diffusion on the Cu(001) surface and nickel-catalyzed methane decomposition, as examples of reactive processes with a bond-length-based CV, and the folding of a long polymer-like chain, using a set of dihedral angles as a CV. Boost factors up to 109, corresponding to a time scale of seconds, could be obtained while still accurately reproducing correct dynamics.

Abstract: Lipid oxidation is associated with several inflammatory and neurodegenerative diseases, but many questions to unravel its effects on biomembranes are still open due to the complexity of the topic. For instance, recent studies indicated that phase-separated domains can have a significant effect on membrane function. It is reported that domain interfaces are “hot spots” for pore formation, but the underlying mechanisms and the effect of oxidation-induced phase separation on membranes remain elusive. Thus, to evaluate the permeability of the membrane coexisting of liquid-ordered (Lo) and liquid-disordered (Ld) domains, we performed atomistic molecular dynamics simulations. Specifically, we studied the membrane permeability of nonoxidized or oxidized homogeneous membranes (single-phase) and at the Lo/Ld domain interfaces of heterogeneous membranes, where the Ld domain is composed of either oxidized or nonoxidized lipids. Our simulation results reveal that the addition of only 1.5% of lipid aldehyde molecules at the Lo/Ld domain interfaces of heterogeneous membranes increases the membrane permeability, whereas their addition at homogeneous membranes does not have any effect. This study is of interest for a better understanding of cancer treatment methods based on oxidative stress (causing among others lipid oxidation), such as plasma medicine and photodynamic therapy.

Abstract: The chemical accessibility of the Celv oxidation state enables redox chemistry to be performed on the naturally coinagemetal -deficient phases CeM1-xSO (M = Cu, Ag). A metastable black compound with the PbFC1 structure type (space group P4/nmm: a = 3.8396(1) angstrom, c = 6.607(4) angstrom, V = 97.40(6) angstrom(3)) and a composition approaching CeSO is obtained by deintercalation of Ag from CeAg0.8SO. High-resolution transmission electron microscopy reveals the presence of large defect-free regions in CeSO, but stacking faults are also evident which can be incorporated into a quantitative model to account for the severe peak anisotropy evident in all the highresolution X-ray and neutron diffractograms of bulk CeSO samples; these suggest that a few percent of residual Ag remains. A strawcolored compound with the filled PbFCI (i.e., ZrSiCuAs- or HfCuSi2type) structure (space group P4/nmm: a = 3.98171(1) angstrom, c = 8.70913(5) angstrom, V = 138.075(1) angstrom 3) and a composition close to LiCeSO, but with small amounts of residual Ag, is obtained by direct reductive lithiation of CeAga8S0 or by insertion of Li into CeSO using chemical or electrochemical means. Computation of the band structure of pure, stoichiometric CeSO predicts it to be a Ce' compound with the 4f-states lying approximately 1 eV above the sulfide-dominated valence band maximum. Accordingly, the effective magnetic moment per Ce ion measured in the CeSO samples is much reduced from the value found for the Ce3+-containing LiCeSO, and the residual paramagnetism corresponds to the Ce3+ ions remaining due to the presence of residual Ag, which presumably reflects the difficulty of stabilizing Ce' in the presence of sulfide (S2-). Comparison of the behavior of CeCu0.8SO with that of CeCu0.8SO reveals much slower reaction kinetics associated with the Cu,_xS layers, and this enables intermediate CeCui LixSO phases to be isolated.

Abstract: Layered Li(M,Li)O2 (where M is a transition metal) ordered rock-salt-type structures are used in advanced metal-ion batteries as one of the best hosts for the reversible intercalation of Li ions. Besides the conventional redox reaction involving oxidation/reduction of the M cation upon Li extraction/insertion, creating oxygen-located holes because of the partial oxygen oxidation increases capacity while maintaining the oxidized oxygen species in the lattice through high covalency of the M–O bonding. Typical degradation mechanism of the Li(M,Li)O2 electrodes involves partially irreversible M cation migration toward the Li positions, resulting in gradual capacity/voltage fade. Here, using LiRhO2 as a model system (isostructural and isoelectronic to LiCoO2), for the first time, we demonstrate an intimate coupling between the oxygen redox and M cation migration. A formation of the oxidized oxygen species upon electrochemical Li extraction coincides with transformation of the layered Li1–xRhO2 structure into the γ-MnO2-type rutile–ramsdellite intergrowth LiyRh3O6 structure with rutile-like [1 × 1] channels along with bigger ramsdellite-like [2 × 1] tunnels through massive and concerted Rh migration toward the empty positions in the Li layers. The oxidized oxygen dimers with the O–O distances as short as 2.26 Å are stabilized in this structure via the local Rh–O configuration reminiscent to that in the μ-peroxo-μ-hydroxo Rh complexes. The LiyRh3O6 structure is remarkably stable upon electrochemical cycling illustrating that proper structural implementation of the oxidized oxygen species can open a pathway toward deliberate employment of the anion redox chemistry in high-capacity/high-voltage positive electrodes for metal-ion batteries. Upon chemical or electrochemical oxidation, layered LiRhO2 shows a unique structural transformation that involves both cation migration and oxidation of oxygen resulting in a stable tunnel-like rutile−ramsdellite intergrowth LiyRh3O6 structure. This structure demonstrates excellent performance with the steady and reversible capacity of ∼200 mAh/g. The stability of LiyRh3O6 is rooted in the accommodation of partially oxidized oxygen species through the formation of short O−O distances that are compatible with the connectivity of RhO6 octahedra.

Abstract: Atomic arrangements in the mixed-valence oxide U3O7 are refined from high-resolution neutron scattering data. The crystallographic model describes a long-range structural order in a U60O140 primitive cell (space group P42/n) containing distorted cuboctahedral oxygen clusters. By combining experimental data and electronic structure calculations accounting for spin–orbit interactions, we provide robust evidence of an interplay between charge localization and the magnetic moments carried by the uranium atoms. The calculations predict U3O7 to be a semiconducting solid with a band gap of close to 0.32 eV, and a more pronounced charge-transfer insulator behavior as compared to the well-known Mott insulator UO2. Most uranium ions (56 out of 60) occur in 9-fold and 10-fold coordinated environments, surrounding the oxygen clusters, and have a tetravalent (24 out of 60) or pentavalent (32 out of 60) state. The remaining uranium ions (4 out of 60) are not contiguous to the oxygen cuboctahedra and have a very compact, 8-fold coordinated environment with two short (2 × 1.93(3) Å) “oxo-type” bonds. The higher Hirshfeld charge and the diamagnetic character point to a hexavalent state for these four uranium ions. Hence, the valence state distribution corresponds to 24/60 × U(IV) + 32/60 U(V) + 4/60 U(VI). The tetravalent and pentavalent uranium ions are predicted to carry noncollinear magnetic moments (with amplitudes of 1.6 and 0.8 μB, respectively), resulting in canted ferromagnetic order in characteristic layers within the overall fluorite-related structure.

Abstract: The mineral Zlatogorite, CuNiSb2, was synthesized in the laboratory for the first time by annealing elements at ambient pressure (CuNiSb2-AP). Rietveld refinement of synchrotron powder X-ray diffraction data indicates that CuNiSb2-AP crystallizes in the NiAs-derived structure (P (3) over bar m1, #164) with Cu and Ni ordering. The structure consists of alternate NiSb6 and CuSb6 octahedral layers via face-sharing. The formation of such structure instead of metal disordered NiAs-type structure (P6(3)/mmc, #194) is validated by the lower energy of the ordered phase by first-principle calculations. Interatomic crystal orbital Hamilton population, electron localization function, and charge density analysis reveal strong Ni-Sb, Cu-Sb, and Cu-Ni bonding and long weak Sb-Sb interactions in CuNiSb2-AP. The magnetic measurement indicates that CuNiSb2-AP is Pauli paramagnetic. First-principle calculations and experimental electrical resistivity measurements reveal that CuNiSb2-AP is a metal. The low Seebeck coefficient and large thermal conductivity suggest that CuNiSb2 is not a potential thermoelectric material. Single crystals were grown by chemical vapor transport. The high pressure sample (CuNiSb2-8 GPa) was prepared by pressing CuNiSb2-AP at 700 degrees C and 8 GPa. However, the structures of single crystal and CuNiSb2-8 GPa are best fit with a disordered metal structure in the P (3) over bar m1 space group, corroborated by transmission electron microscopy.