Abstract: Two-dimensional group-V elemental materials have attracted widespread attention due to their nonzero band gap while displaying high electron mobility. Using first-principles calculations, we propose a series of new elemental bilayers with group-V elements (Bi, Sb, As). Our study reveals the dynamical stability of four-, six-, and eight-atom ring structures, demonstrating their possible coexistence in such bilayer systems. The proposed structures for Sb and As are large-gap semiconductors that are potentially interesting for applications in future nanodevices. The Bi structures have nontrivial topological properties with a direct nontrivial band gap. The nontrivial gap is shown to arise from a band inversion at the Brillouin zone center due to the strong intrinsic spin-orbit coupling in Bi atoms. Moreover, we demonstrate the possibility of tuning the properties of these materials by enhancing the ratio of six-atom rings to four-and eight-atom rings, which results in wider nontrivial band gaps and lower formation energies.

Abstract: Two-dimensional materials (2DMs) continue to attract a lot of attention, particularly for their extreme flexibility and superior thermal properties. Molecular dynamics simulations are among the most powerful methods for computing these properties, but their reliability depends on the accuracy of interatomic interactions. While first principles approaches provide the most accurate description of interatomic forces, they are computationally expensive. In contrast, classical force fields are computationally efficient, but have limited accuracy in interatomic force description. Machine learning interatomic potentials, such as Gaussian Approximation Potentials, trained on density functional theory (DFT) calculations offer a compromise by providing both accurate estimation and computational efficiency. In this work, we present a systematic procedure to develop Gaussian approximation potentials for selected 2DMs, graphene, buckled silicene, and h-XN (X = B, Al, and Ga, as binary compounds) structures. We validate our approach through calculations that require various levels of accuracy in interatomic interactions. The calculated phonon dispersion curves and lattice thermal conductivity, obtained through harmonic and anharmonic force constants (including fourth order) are in excellent agreement with DFT results. HIPHIVE calculations, in which the generated GAP potentials were used to compute higher-order force constants instead of DFT, demonstrated the first-principles level accuracy of the potentials for interatomic force description. Molecular dynamics simulations based on phonon density of states calculations, which agree closely with DFT-based calculations, also show the success of the generated potentials in high-temperature simulations.

Abstract: Two-dimensional materials have been attracting increasing interests because of their outstanding properties for Lithium-ion battery applications. In particular, a material family called MXenes (Mn+1Cn, where n = 1, 2, 3) have been recently attracted immense interest in this respect due to their incomparable fast-charging properties and high capacity promises. In this article, we review the state-of-the-art computational progress on Li-ion battery applications of MXene materials in accordance with our systematical DFT calculations. Structural, mechanical, dynamical, and electrical properties of 20 distinct MXene (M: Sc, Ti, V, Cr, Nb, Mo, Hf, Ta, W, and Zr) have been discussed. The battery performances of these MXene monolayers are further investigated by Li-ion binding energies, open circuit voltage values, and Li migration energy barriers. The experimental and theoretical progress up to date demonstrates particularly the potential of non-terminated or pristine MXene materials in Li ion-storage applications. Stability analyses show most of the pristine MXenes should be achievable, however susceptible to the development progress on the experimental growth procedures. Among pristine MXenes, Ti2C, V2C, Sc2C, and Zr2C compounds excel with their high charge/discharge rate prospect due to their extremely low Li diffusion energy barriers. Considering also their higher predicted gravimetric capacities, Sc, Ti, V, and Zr containing MXenes are more promising for their utilization in energy storage applications.

Abstract: Two-dimensional systems of inverse patchy colloids modeled as disks with a central charge and having their surface decorated with oppositely pointlike charged patches are investigated using molecular dynamics simulations. The self-assembly of the patchy colloids leads to diverse ground state configurations ranging from crystalline arrangements of monomers to linear clusters, ramified linear clusters and to percolated configurations. Two structural phase diagrams are constructed: (1) as a function of the net charge and area fraction, and (2) as a function of the net charge and the range of the pair interaction potential. An interesting reentrant percolation transition is obtained as a function of the net charge of the colloids. We identify distinct mechanisms that lead to the percolation transition.

Abstract: Two-dimensional Transition Metal Dichalcogenides (TMDs) possessing extraordinary physical properties at reduced dimensionality have attracted interest due to their promise in electronic and optical device applications. However, TMD monolayers can show a broad range of different properties depending on their crystal phase; for example, H phases are usually semiconductors, while the T phases are metallic. Thus, controlling phase transitions has become critical for device applications. In this study, the energetically low-lying crystal structures of pristine and Janus TMDs are investigated by using ab initio Nudged Elastic Band and molecular dynamics simulations to provide a general explanation for their phase stability and transition properties. Across all materials investigated, the T phase is found to be the least stable and the H phase is the most stable except for WTe2, while the T' and T '' phases change places according to the TMD material. The transition energy barriers are found to be large enough to hint that even the higher energy phases are unlikely to undergo a phase transition to a more stable phase if they can be achieved except for the least stable T phase, which has zero barrier towards the T ' phase. Indeed, in molecular dynamics simulations the thermodynamically least stable T phase transformed into the T ' phase spontaneously while in general no other phase transition was observed up to 2100 K for the other three phases. Thus, the examined T ', T '' and H phases were shown to be mostly stable and do not readily transform into another phase. Furthermore, so-called mixed phase calculations considered in our study explain the experimentally observed lateral hybrid structures and point out that the coexistence of different phases is strongly stable against phase transitions. Indeed, stable complex structures such as metal-semiconductor-metal architectures, which have immense potential to be used in future device applications, are also possible based on our investigation.

Abstract: Typical Raman spectra of transition-metal dichalcogenides (TMDs) display two prominent peaks, E-2g and A(1g), that are well separated from each other. We find that these modes are degenerate in bulk WSe2 yielding one single Raman peak in contrast to other TMDs. As the dimensionality is lowered, the observed peak splits in two. In contrast, our ab initio calculations predict that the degeneracy is retained even for WSe2 monolayers. Interestingly, for minuscule biaxial strain, the degeneracy is preserved, but once the crystal symmetry is broken by a small uniaxial strain, the degeneracy is lifted. Our calculated phonon dispersion for uniaxially strained WSe2 shows a good match to the measured Raman spectrum, which suggests that uniaxial strain exists in WSe2 flakes, possibly induced during the sample preparation and/or as a result of the interaction between WSe2 and the substrate. Furthermore, we find that WSe2 undergoes an indirect-to-direct band-gap transition from bulk to monolayers, which is ubiquitous for semiconducting TMDs. These results not only allow us to understand the vibrational and electronic properties of WSe2, but also point to effects of the interaction between the monolayer TMDs and the substrate on the vibrational and electronic properties. DOI: 10.1103/PhysRevB.87.165409

Abstract: Ultranarrow superconductors are in the strong quantum confinement regime with formation of multiple coherent condensates associated with the many subbands of the electronic structure. Here, we analyze the multiband BCS-BEC crossover induced by the chemical potential tuned close to a subband bottom, in correspondence of a superconducting shape resonance. The evolution of the condensate fraction and of the pair correlation length in the ground state as functions of the chemical potential demonstrates the tunability of the BCS-BEC crossover for the condensate component of the selected subband. The extension of the crossover regime increases when the pairing strength and/or the characteristic energy of the interaction get larger. Our results indicate the coexistence of large and small Cooper pairs in the crossover regime, leading to the optimal parameter configuration for high transition temperature superconductivity.

Abstract: Ultrathin two-dimensional Janus-type platinum dichalcogenide crystals formed by two different atoms at opposite surfaces are investigated by performing state-of-the-art density functional theory calculations. First, it is shown that single-layer PtX2 structures (where X = S, Se, or Te) crystallize into the dynamically stable IT phase and are indirect band gap semiconductors. It is also found that the substitutional chalcogen doping in all PtX2 structures is favorable via replacement of surface atoms with a smaller chalcogen atom, and such a process leads to the formation of Janus-type platinum dichalcogenides (XPtY, where X and Y stand for S, Se, or Te) which are novel single-layer crystals. While all Janus structures are indirect band gap semiconductors as their binary analogues, their Raman spectra show distinctive features that stem from the broken out-of-plane symmetry. In addition, it is revealed that the construction of Janus crystals enhances the piezoelectric constants of PtX2 crystals significantly both in the in plane and in the out-of-plane directions. Moreover, it is shown that vertically stacked van der Waals heterostructures of binary and ternary (Janus) platinum dichalcogenides offer a wide range of electronic features by forming bilayer heterojunctions of type-I, type-II, and type-III, respectively. Our findings reveal that Janus-type ultrathin platinum dichalcogenide crystals are quite promising materials for optoelectronic device applications.

Abstract: Understanding the appearance of local electrical connectivity in liquid filled layered graphene oxide subjected to an external electric field is important to design electrically controlled smart permeable devices and also to gain insight into the physics behind electrical effects on confined water permeation. Motivated by recent experiments [K. G. Zhou et al. Nature (London) 559, 236 (2018)], we introduce a new model with random percolating paths for electrical connectivity in micron thick water filled layered graphene oxide, which mimics parallel resistors connected across the top and bottom electrodes. We find that a strong nonuniform radial electric field of the order similar to 10-50 mV/nm can be induced between layers depending on the current flow through the formed conducting paths. The maxima of the induced fields are not necessarily close to the electrodes and may be localized in the middle region of the layered material. The emergence of electrical connectivity and the associated electrical effects have a strong influence on the surrounding fluid in terms of ionization and wetting which subsequently determines the permeation properties.

Abstract: Understanding the coupling between spin and phonons is critical for controlling the lattice thermal conductivity (kappa ( l )) in magnetic materials, as we demonstrate here for CrX3 (X = Br and I) monolayers. We show that these compounds exhibit large spin-phonon coupling (SPC), dominated by out-of-plane vibrations of Cr atoms, resulting in significantly different phonon dispersions in ferromagnetic (FM) and paramagnetic (PM) phases. Lattice thermal conductivity calculations provide additional evidence for strong SPC, where particularly large kappa ( l ) is found for the FM phase. Most strikingly, PM and FM phases exhibit radically different behavior with tensile strain, where kappa ( l ) increases with strain for the PM phase, and strongly decreases for the FM phase-as we explain through analysis of phonon lifetimes and scattering rates. Taken all together, we uncover the high significance of SPC on the phonon transport in CrX3 monolayers, a result extendable to other 2D magnetic materials, that will be useful in further design of thermal spin devices.

Abstract: Using a combination of the phenomenological GinzburgLandau theory and micromagnetic simulations, we study properties of a superconducting film with an array of soft magnetic dots on top. An external in-plane magnetic field gradually drives the magnets from an out-of-plane or magnetic vortex state to an in-plane single-domain state, which changes spatially the distribution of the superconducting condensate. If induced by the magnets, the vortexantivortex molecules exhibit rich transitions as a function of the applied in-plane field. At the same time, we show how the magnetic dots act as very effective dynamic pinning centers for vortices in an applied perpendicular magnetic field.

Abstract: Using a particle based model, we investigate the skyrmion dynamical behavior in a channel where the upper wall contains divots of one depth and the lower wall contains divots of a different depth. Under an applied driving force, skyrmions in the channels move with a finite skyrmion Hall angle that deflects them toward the upper wall for -x direction driving and the lower wall for +x direction driving. When the upper divots have zero height, the skyrmions are deflected against the flat upper wall for -x direction driving and the skyrmion velocity depends linearly on the drive. For +x direction driving, the skyrmions are pushed against the lower divots and become trapped, giving reduced velocities and a nonlinear velocity-force response. When there are shallow divots on the upper wall and deep divots on the lower wall, skyrmions get trapped for both driving directions; however, due to the divot depth difference, skyrmions move more easily under -x direction driving, and become strongly trapped for +x direction driving. The preferred -x direction motion produces what we call a Magnus diode effect since it vanishes in the limit of zero Magnus force, unlike the diode effects observed for asymmetric sawtooth potentials. We show that the transport curves can exhibit a series of jumps or dips, negative differential conductivity, and reentrant pinning due to collective trapping events. We also discuss how our results relate to recent continuum modeling on a similar skyrmion diode system.

Abstract: Using a particle-based model, we examine the collective dynamics of skyrmions interacting with a funnel potential under dc driving as the skyrmion density and relative strength of the Magnus and damping terms are varied. For driving in the easy direction, we find that increasing the skyrmion density reduces the average skyrmion velocity due to jamming of skyrmions near the funnel opening, while the Magnus force causes skyrmions to accumulate on one side of the funnel array. For driving in the hard direction, there is a critical skyrmion density below which the skyrmions become trapped. Above this critical value, a clogging effect appears with multiple depinning and repinning states where the skyrmions can rearrange into different clogged configurations, while at higher drives, the velocity-force curves become continuous. When skyrmions pile up near the funnel opening, the effective size of the opening is reduced and the passage of other skyrmions is blocked by the repulsive skyrmion-skyrmion interactions. We observe a strong diode effect in which the critical depinning force is higher and the velocity response is smaller for hard direction driving. As the ratio of Magnus force to dissipative term is varied, the skyrmion velocity varies in a non-linear and non-monotonic way due to the pile up of skyrmions on one side of the funnels. At high Magnus forces, the clogging effect for hard direction driving is diminished.

Abstract: Using advanced three-dimensional simulations, we show that an Abrikosov vortex, trapped inside a cavity perpendicular to an artificial Josephson junction, can serve as a very efficient source for generation of Josephson vortex-antivortex pairs in the presence of the applied electric current. In such a case, the nucleation rate of the pairs can be tuned in a broad range by an out-of-plane ac magnetic field in a broad range of frequencies. This parametrically amplified vortex-antivortex nucleation can be considered as a macroscopic analog of the dynamic Casimir effect, where fluxon pairs mimic the photons and the ac magnetic field plays the role of the oscillating mirrors. The emerging vortex pairs in our system can be detected by the pronounced features in the measured voltage characteristics, or through the emitted electromagnetic radiation, and exhibit resonant dynamics with respect to the frequency of the applied magnetic field. Reported tunability of the Josephson oscillations can be useful for developing high-frequency emission devices.

Abstract: Using an exact-diagonalization approach we show that one- and two-electron InAs quantum dots exhibit an avoided crossing in the energy spectra that is induced by the spin-orbit coupling in the presence of an in-plane external magnetic field. The width of the avoided crossings depends strongly on the orientation of the magnetic field, which reveals the intrinsic anisotropy of the spin-orbit-coupling interactions. We find that for specific orientations of the magnetic field avoided crossings vanish. A value of this orientation can be used to extract the ratio of the strength of Rashba and Dresselhaus interactions. The spin-orbit anisotropy effects for various geometries and orientations of the confinement potential are discussed. Our analysis explains the physics behind the recent measurements performed on a gated self-assembled quantum dot [ S. Takahashi et al. Phys. Rev. Lett. 104 246801 (2010)].