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 ditelluride (MoTe2) nanosheets have displayed intriguing physicochemical properties and opto-electric characteristics as a result of their tunable and
small band gap (Eg ∼ 1 eV), facilitating concurrent electron and hole transport. Despite the numerous efforts devoted to the development of p-type MoTe2 field-effect transistors (FETs), the presence of tellurium (Te) point vacancies has caused serious reliability issues. Here, we overcome this major
limitation by treating the MoTe2 surface with thiolated molecules to heal Te vacancies. Comprehensive materials and electrical characterizations provided unambiguous evidence for the efficient chemisorption of butanethiol. Our thiol-treated MoTe2 FET exhibited a 10-fold increase in hole current and a positive threshold voltage shift of 25 V, indicative of efficient hole carrier doping. We demonstrated that our powerful molecular engineering strategy can be extended to the controlled formation of van der Waals heterostructures by developing an n-SnS2/thiol-MoTe2 junction FET (thiol-JFET). Notably, the thiol-JFET exhibited a significant negative photoresponse with a responsivity of 50 A W−1 and a fast response time of 80 ms based on band-to-band tunneling. More interestingly, the
thiol-JFET displayed a gate tunable trimodal photodetection comprising two photoactive modes (positive and negative photoresponse) and one photoinactive mode. These findings underscore the potential of molecular engineering approaches in
enhancing the performance and functionality of MoTe2-based nanodevices as key components in advanced 2D-based optoelectronics.

Abstract: Molecular-dynamic simulations were performed in order to investigate the melting processes of isotropically confined binary systems. We considered two species of particles, which differ by their amount of electric charge. A Lindemann type of criterion was used to determine the angular melting temperature. We demonstrate that the magic-to-normal cluster transition can evolve in two distinct ways, that is, through a structural phase transition of the first order or via a smooth transition where an increase of the shells' width leads to a continuous decreasing mechanical stability of the system. Moreover, for large systems, we demonstrate that the internal cluster exerts a minor effect on the mechanical stability of the external shell. Furthermore, we show that highly symmetric configurations, such as those found for multiple ring structures, have large mechanical stability, i.e., high angular melting temperature.

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 improved TersoffBrenner potential parameters were performed to investigate Cl+ etching of a {2×1} reconstructed Si(100) surface. Steady-state Si etching accompanying the Cl coverage of the surface is observed. Furthermore, a steady-state chlorinated reaction layer is formed. The thickness of this reaction layer is found to increase with increasing energy. The stoichiometry of SiClx species in the reaction layer is found to be SiCl:SiCl2:SiCl3 = 1.0:0.14:0.008 at 50 eV. These results are in excellent agreement with available experimental data. While elemental Si products are created by physical sputtering, most SiClx (0<x<4) etch products are produced by chemical-enhanced physical sputtering.

Abstract: Molecular dynamics simulations are used to investigate the structural properties of a two-dimensional ensemble of magnetic rods, which are modeled as aligned single dipolar beads. The obtained self-assembled configurations can be characterized as (1) clusters, (2) percolated, and (3) ordered structures, and their structural properties are investigated in detail. By increasing the aspect ratio of the magnetic rods, we show that the percolation transition is suppressed due to the reduced mobility of the rods in two dimensions. Such a behavior is opposite to the one observed in three dimensions. A magnetic bulk phase is found with local ferromagnetic order and an unusual nonmonotonic behavior of the nematic order is observed.

Abstract: Molecular dynamics simulations are performed on capped and uncapped PbSe nanocrystals, employing newly developed classical interaction potentials. Here, we show that two uncapped nanocrystals fuse efficiently via direct surface attachment, even if they are initially misaligned. In sharp contrast to the general belief, interparticle dipole interactions do not play a significant role in this oriented attachment process. Furthermore, it is shown that presumably polar, capped PbSe{111} facets are never fully Pb- or Se-terminated.

Abstract: Molecular dynamics simulations are carried out for describing growth of Pd and PdO nanoclusters using the ReaxFF force field. The resulting nanocluster structures are successfully compared to those of nanoclusters experimentally grown in a gas aggregation source. The PdO structure is quasi-crystalline as revealed by high resolution transmission microscope analysis for experimental PdO nanoclusters. The role of the nanocluster temperature in the molecular dynamics simulated growth is highlighted.

Abstract: Molecular dynamics (MD) simulations were performed to provide atomic scale insight in the initial interaction between hydroxyl radicals (OH) and peptide systems in solution. These OH radicals are representative reactive oxygen species produced by cold atmospheric plasmas. The use of plasma for biomedical applications is gaining increasing interest, but the fundamental mechanisms behind the plasma modifications still remain largely elusive. This study helps to gain more insight in the underlying mechanisms of plasma medicine but is also more generally applicable to peptide oxidation, of interest for other applications. Combining both reactive and nonreactive MD simulations, we are able to elucidate the reactivity of the amino acids inside the peptide systems and their effect on their structure up to 1 μs. Additionally, experiments were performed, treating the simulated peptides with a plasma jet. The computational results presented here correlate well with the obtained experimental data and highlight the importance of the chemical environment for the reactivity of the individual amino acids, so that specific amino acids are attacked in higher numbers than expected. Furthermore, the long time scale simulations suggest that a single oxidation has an effect on the 3D conformation due to an increase in hydrophilicity and intra- and intermolecular interactions.

Abstract: Moire patterns in the pseudo-magnetic field and in the strain profile of graphene (GE) when put on top of a hexagonal lattice substrate are predicted from elasticity theory. The van der Waals interaction between GE and the substrate induces out-of-plane deformations in graphene which results in a strain field, and consequently in a pseudo-magnetic field. When the misorientation angle is about 0.5 degrees, a three-fold symmetric strain field is realized that results in a pseudo-magnetic field very similar to the one proposed by F. Guinea, M. I. Katsnelson, and A. K. Geim [Nature Phys. 6, 30 (2010)]. Our results show that the periodicity and length of the pseudo-magnetic field can be tuned in GE by changing the misorientation angle and substrate adhesion parameters and a considerable energy gap (23 meV) can be obtained due to out-of-plane deformation of graphene which is in the range of recent experimental measurements (20-30 meV). (C) 2014 AIP Publishing LLC.

Abstract: MOF-5 nanoparticles were prepared by mixing a solution of [Zn4O(C6H5COO)(6)] with a solution of benzene-1,4-dicarboxylic acid in DMF at ambient conditions. The former species mimics as a secondary building unit (SBU), and the latter acts as linker. Mixing of the two solutions induced the formation of MOF-5 nanoparticles in dilute suspension. The applied conditions were identified as suitable for a closer investigation of the particle formation process by combined light and small angle neutron scattering (SANS). Scattering analysis revealed a significant impact of the molar ratio of the two components in the reaction mixture. Excessive use of the building unit slowed down the process. A similar effect was observed upon addition of 4n-decylbenzoic acid, which is supposed to act as a modulator. The formation mechanism leads to initial intermediates, which turn into cubelike nanoparticles with a diameter of about 60-80 nm. This initial stage is followed by an extended formation period, where nucleation proceeds over hours, leading to an increasing number of nanoparticles with the same final size of 60-80 nm.