Abstract: Motivated by recent studies demonstrating a high strength – high ductility potential of nano-scale metallic glass samples, the mechanical response of freestanding Zr65Ni35 film with sub-micron thickness has been investigated by combining advanced on-chip tensile testing and electron microscopy. Large deformation up to 15% is found for specimen thicknesses below 500 nm with variations depending on specimen size and frame compliance. The deformation is homogenous until fracture, with no evidence of shear banding. The yield stress is doubled when decreasing the specimen cross-section, reaching ~3 GPa for small cross-sections. The fracture strain variation is related to both the stability of the test device and to the specimen size. The study concludes on clear disconnect between the mechanisms controlling the onset of plasticity and the fracture process.

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: Most of the studied two-dimensional (2D) materials have been obtained by exfoliation of van der Waals crystals. Recently, there has been growing interest in fabricating synthetic 2D crystals which have no layered bulk analogues. These efforts have been focused mainly on the surface growth of molecules in high vacuum. Here, we report an approach to making 2D crystals of covalent solids by chemical conversion of van der Waals layers. As an example, we used 2D indium selenide (InSe) obtained by exfoliation and converted it by direct fluorination into indium fluoride (InF3), which has a nonlayered, rhombohedral structure and therefore cannot possibly be obtained by exfoliation. The conversion of InSe into InF3 is found to be feasible for thicknesses down to three layers of InSe, and the obtained stable InF3 layers are doped with selenium. We study this new 2D material by optical, electron transport, and Raman measurements and show that it is a semiconductor with a direct bandgap of 2.2 eV, exhibiting high optical transparency across the visible and infrared spectral ranges. We also demonstrate the scalability of our approach by chemical conversion of large-area, thin InSe laminates obtained by liquid exfoliation, into InF3 films. The concept of chemical conversion of cleavable thin van der Waals crystals into covalently bonded noncleavable ones opens exciting prospects for synthesizing a wide variety of novel atomically thin covalent crystals.

Abstract: Most farm adaptations are reactive actions that run the risk of locking farm systems into suboptimal long-term trajectories. This is especially the case with regard to water management as water scarcity will be aggravated by climate change. This paper looks into farm irrigation choices in combination with crop choices because a proper crop choice has the potential to reduce water requirements. It proposes an extended Ricardian model to capture multiple adaptation decisions explicitly. The new simultaneous irrigation-crop farm decision model uses spatially detailed farm-level data of over 18,000 European farms on irrigation and seven different crop choices. The analysis shows that larger farmers and farmers in less water-scarce regions that use irrigation are more sensitive to temperature increases than rain-fed agriculture. This might be explained by the fact that these farmers do not experience the real cost of water scarcity because of which they take less efficient decisions.

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: Moore’s Law has driven the electronics industry for the past half century. However, the doubling of transistors about every two years is beginning to break down, owing to fundamental limits that arise as they approach the atomic length. As a result, the search for new pathways for electronics has become crucial. Among potential candidates, the discovery of magnetic textures known as skyrmions has attracted considerable interest and attention in spintronic technology, which relies on both the electron charge and its spin. The unusual topological and particle-like behavior launched skyrmions into the spotlight of scientific research. Topological protected stability, nanoscale size, and low driving currents needed to move them make skyrmions promising candidates for future consumer nanoelectronics. Recent advances in the field have provided all of the basic functions needed for carrying and processing information. In this thesis, we procure to advance the current understanding of skyrmion physics, and explore their potential to replace conventional electronics technology. First, the fundamental properties and lifetimes of racetrack skyrmions at room temperature are investigated. We discover that skyrmions can easily collapse at the boundary in laterally finite systems, and propose ways to improve their stability for constrained geometries. Then, pinning of single skyrmions on atomic defects of distinct origins are studied. We reveal that the preferred pinning positions depend on the skyrmion size and type of defect being considered, and discuss applications where control of skyrmions by defects is of particular interest. Next, we explore other magnetic configurations that can compete with skyrmions when considering new materials, and describe a previously unseen mechanism for collapse of skyrmions into cycloidal spin backgrounds. Finally, switching and interactions between skyrmions with distinct topologies are reported. We find that skyrmions transition to higher or lower topologies by absorbing a unit spin texture. The interactions between skyrmions of different topological charges can be attractive or repulsive, leading to the formation of arranged clusters. We conclude with a numerical library for simulating magnetic skyrmions in various scenarios.

Abstract: Monte Carlo (MC) methods have a long-standing history as partners of molecular dynamics (MD) to simulate the evolution of materials at the atomic scale. Among these techniques, the uniform-acceptance force-bias Monte Carlo (UFMC) method [ G. Dereli Mol. Simul. 8 351 (1992)] has recently attracted attention [ M. Timonova et al. Phys. Rev. B 81 144107 (2010)] thanks to its apparent capacity of being able to simulate physical processes in a reduced number of iterations compared to classical MD methods. The origin of this efficiency remains, however, unclear. In this work we derive a UFMC method starting from basic thermodynamic principles, which leads to an intuitive and unambiguous formalism. The approach includes a statistically relevant time step per Monte Carlo iteration, showing a significant speed-up compared to MD simulations. This time-stamped force-bias Monte Carlo (tfMC) formalism is tested on both simple one-dimensional and three-dimensional systems. Both test-cases give excellent results in agreement with analytical solutions and literature reports. The inclusion of a time scale, the simplicity of the method, and the enhancement of the time step compared to classical MD methods make this method very appealing for studying the dynamics of many-particle systems.

Abstract: Monolayers of group VI transition metal dichalcogenides possess direct gaps in the visible spectrum with the exception of MoTe2, where its gap is suitably located in the infrared region but its stability is of particular interest, as tellurium compounds are acutely sensitive to oxygen exposure. Here, our environmental (time-dependent) measurements reveal two distinct effects on MoTe2 monolayers: For weakly luminescent monolayers, photoluminescence signal and optical contrast disappear, as if they are decomposed, but yet remain intact as evidenced by AFM and Raman measurements. In contrast, strongly luminescent monolayers retain their optical contrast for a prolonged amount of time, while their PL peak blue-shifts and PL intensity saturates to slightly lower values. Our X-ray photoelectron spectroscopy measurements and DFT calculations suggest that the presence of defects and functionalization of these defect sites with O-2 molecules strongly dictate their material properties and aging response by changing the excitonic dynamics due to deep or shallow states that are created within the optical band gap. Presented results not only shed light on environmental effects on fundamental material properties and excitonic dynamics of MoTe2 monolayers but also highlight striking material transformation for metastable 20 systems such as WTe2, silicone, and phosphorene.

Abstract: Monolayer chromium-trihalides, the archetypal two-dimensional (2D) magnetic materials, are readily suggested as a promising platform for high-frequency magnonics. Here we detail the spin-wave properties of monolayer CrBr₃and CrI₃, using spin-dynamics simulations parametrized from the first principles. We reveal that spin-wave dispersion can be tuned in a broad range of frequencies by strain, paving the way towards flexo-magnonic applications. We further show that ever-present halide vacancies in these monolayers host sufficiently strong Dzyaloshinskii-Moriya interaction to scatter spin-waves, which promotes design of spin-wave guides by defect engineering. Finally we discuss the spectra of spin-waves propagating across a moiré-periodic modulation of magnetic parameters in a van der Waals heterobilayer, and show that the nanoscale moiré periodicities in such samples are ideal for realization of a magnonic crystal in the terahertz frequency range. Recalling the additional tunability of magnetic 2D materials by electronic gating, our results situate these systems among the front-runners for prospective high-frequency magnonic applications.

Abstract: Monitoring and modelling tools may improve irrigation strategies in precision agriculture. We used non-invasive soil moisture monitoring, a crop growth and a soil hydrological model to predict soil water content fluctuations and crop yield in a heterogeneous sandy grassland soil under supplementary irrigation. The sensitivity of the soil hydrological model to hydraulic parameters, water stress, crop yield and lower boundary conditions was assessed after integrating models. Free drainage and incremental constant head conditions were implemented in a lower boundary sensitivity analysis. A time-dependent sensitivity analysis of the hydraulic parameters showed that changes in soil water content are mainly affected by the soil saturated hydraulic conductivity K-s and the Mualem-van Genuchten retention curve shape parameters n and alpha. Results further showed that different parameter optimization strategies (two-, three-, four- or six-parameter optimizations) did not affect the calculated water stress and water content as significantly as does the bottom boundary. In this case, a two-parameter scenario, where K-s was optimized for each layer under the condition of a constant groundwater depth at 135-140 cm, performed best. A larger yield reduction, and a larger number and longer duration of stress conditions occurred in the free drainage condition as compared to constant boundary conditions. Numerical results showed that optimal irrigation scheduling using the aforementioned water stress calculations can save up to 12-22 % irrigation water as compared to the current irrigation regime. This resulted in a yield increase of 4.5-6.5 %, simulated by the crop growth model.

Abstract: Molybdenum carbide has various applications for which studying the material using classical molecular dynamics simulations would be valuable. Here, we develop an empirical potential within the Tersoff formalism using particle swarm optimization for the orthorhombic phase of Mo2C. The developed potential is shown to predict lattice constants, elastic properties, and equation of state results that are consistent with current and previously reported results from experiments and first principles calculations. We demonstrate the potential with simulations of indentation using multiple indenter sizes that load and unload in three different directions relative to the crystallographic lattice of orthorhombic Mo2C. Direction-dependent force-displacement trends are analyzed and explained in terms of the spatial distributions of stress and strain within the material during indentation. This study reveals the anisotropic elasticity of orthorhombic Mo2C and, more generally, provides researchers with a new empirical potential that can be used to explore the properties and behavior of the material going forward.

Abstract: Molecular dynamics simulations were used to investigate the size-dependent melting mechanism of nickel nanoclusters of various sizes. The melting process was monitored by the caloric curve, the overall cluster Lindemann index, and the atomic Lindemann index. Size-dependent melting temperatures were determined, and the correct linear dependence on inverse diameter was recovered. We found that the melting mechanism gradually changes from dynamic coexistence melting to surface melting with increasing cluster size. These findings are of importance in better understanding carbon nanotube growth by catalytic chemical vapor deposition as the phase state of the catalyst nanoparticle codetermines the growth mechanism.

Abstract: Molecular dynamics simulations using the Brenner potential have been performed to investigate reaction mechanisms of various hydrocarbon radicals with low kinetic energies on amorphous hydrogenated carbon (a-C:H) surfaces and to simulate thin a-C:H film growth. Experimental data from an expanding thermal plasma setup were used as input for the simulations. The hydrocarbon reaction mechanisms were studied both during growth of the films and on a set of surface sites specific for a-C:H surfaces. Thin film growth was studied using experimentally detected growth species. It is found that the reaction mechanisms and sticking coefficients are dependent on the specific surface sites, and the structural properties of the growth radicals. Furthermore, it is found that thin a-C:H films can be densified using an additional H-flux towards the substrate.

Abstract: Molecular dynamics simulations are used to study the effects of different cations on the permeation of charged polymers through flat capillaries with heights below 2 nm. Interestingly, we found that, despite being monovalent, Li+ , Na+ , and K+ cations have different effects on polymer permeation, which consequently affects their transmission speed throughout those capillaries. We attribute this phenomenon to the interplay of the cations' hydration free energies and the hydrodynamic drag in front of the polymer when it enters the capillary. Different alkali cations exhibit different surface versus bulk preferences in small clusters of water under the influence of an external electric field. This paper presents a tool to control the speed of charged polymers in confined spaces using cations.

Abstract: Molecular dynamics simulations are used to investigate the effect of an AFM tip when indenting graphene nano bubbles filled by a noble gas (i.e. He, Ne and Ar) up to the breaking point. The failure points resemble those of viral shells as described by the Foppl-von Karman (FvK) dimensionless number defined in the context of elasticity theory of thin shells. At room temperature, He gas inside the bubbles is found to be in the liquid state while Ne and Ar atoms are in the solid state although the pressure inside the nano bubble is below the melting pressure of the bulk. The trapped gases are under higher hydrostatic pressure at low temperatures than at room temperature.