Abstract: We analyzed pattern formation and identified various morphologies in a system of particles interacting through a non-monotonic potential with a competing range interaction characterized by a repulsive core (r < r(c)) and an attractive tail (r > r(c)), using molecular-dynamics simulations. Depending on parameters, the interaction potential models the inter-particle interaction in various physical systems ranging from atoms, molecules and colloids to vortices in low kappa type-II superconductors and in recently discovered 'type-1.5' superconductors. We constructed a 'morphology diagram' in the plane 'critical radius r(c)-density n' and proposed a new approach to characterizing the different types of patterns. Namely, we elaborated a set of quantitative criteria in order to identify the different pattern types, using the radial distribution function (RDF), the local density function and the occupation factor.

Abstract: We use kinetic Monte Carlo simulations to investigate current fluctuations in boundary driven generalized exclusion processes, in different dimensions. Simulation results are in full agreement with predictions based on the additivity principle and the macroscopic fluctuation theory. The current statistics are independent of the shape of the contacts with the reservoirs, provided they are macroscopic in size. In general, the current distribution depends on the spatial dimension. For the special cases of the symmetric simple exclusion process and the zero-range process, the current statistics are the same for all spatial dimensions.

Abstract: In superconducting films under an applied dc current, we analyze experimentally and theoretically the influence of engineered pinning on the vortex velocity at which the flux-flow dissipation undergoes an abrupt transition from low to high resistance. We argue, based on a nonuniform distribution of vortex velocity in the sample, that in strongly disordered systems the mean critical vortex velocity for flux-flow instability (i) has a nonmonotonic dependence on magnetic field and (ii) decreases as the pinning strength is increased. These findings challenge the generally accepted microscopic model of Larkin and Ovchinnikov (1979 J. Low. Temp. Phys. 34 409) and all subsequent refinements of this model which ignore the presence of pinning centers.

Abstract: Nanoparticle charging in a capacitively coupled radio frequency discharge in argon is studied using a particle in cell Monte Carlo collisions method. The plasma parameters and dust potential were calculated self-consistently for different unmovable dust profiles. A new method for definition of the dust floating potential is proposed, based on the information about electron and ion energy distribution functions, obtained during the kinetic simulations. This approach provides an accurate balance of the electron and ion currents on the dust particle surface and allows us to precisely calculate the dust floating potential. A comparison of the obtained floating potentials with the results of the traditional orbital motion limit (OML) theory shows that in the presence of the ion resonant charge exchange collisions, even when the OML approximation is valid, its results are correct only in the region of a weak electric field, where the ion drift velocity is much smaller than the thermal one. With increasing ion drift velocity, the absolute value of the calculated dust potential becomes significantly smaller than the theory predicts. This is explained by a non-Maxwellian shape of the ion energy distribution function for the case of fast ion drift.

Abstract: Using three-dimensional (3D) numerical discretization of the GinzburgLandau (GL) equations, we investigate the superconducting state of a sphere with a piercing hole in the presence of a magnetic field. In the case of samples with central perforation, in axially applied homogeneous magnetic field, we realized unconventional vortex states of broken symmetry due to complex, 3D competing interactions, which depend on the GL parameter ê. For certain sizes of the sample, non-hysteretic multi-vortex entry and exit is predicted with the non-existence of some vorticities as stable states. In a tilted magnetic field, we studied the gradual transformation of 3D flux patterns into 1D vortex chains, where vortices align along the perforation, and the evolvement of the multi-vortex entry as well. We analyze the flux-guiding ability of the hole in a tilted field, which leads to fractional flux response in magnetization M(H) curves.

Abstract: We theoretically investigate, through Ginzburg?Landau simulations, the possibility to induce an open circuit voltage in absence of applied current, by dragging superconducting vortices with a dynamic pinning array as for instance that created by a nearby sliding vortex lattice or moving laser spots. Different dynamic regimes, such as synchronous vortex motion or dynamic vortex chains consisting of laggard vortices, can be observed by varying the velocity of the sliding pinning potential and the applied magnetic field. Additionally, due to the edge barrier, significantly different induced voltage is found depending on whether the vortices are dragged along the superconducting strip or perpendicular to the lateral edges. The output voltage in the proposed mesoscopic superconducting dynamo can be tuned by varying size, density and directions of the sliding pinning potential.

Abstract: Dumbbell-like structures are recently found to be energetically favored in group IV two-dimensional (2D) materials, exhibiting rich physics and many interesting properties. In this paper, using first-principles calculations, we have investigated the oxidized form of the hexagonal honeycomb (ODB-h) and zigzag dumbbell silicene (ODB-z). We confirm that both oxidization processes are energetically favorable, and their phonon spectra further demonstrate the dynamic stability. Contrary to the pristine dumbbell silicene structures (PDB-h and PDB-z silicene), these oxidized products ODB-h and ODB-z silicene are both semimetals with Dirac cones at the Fermi level. The Dirac cones of ODB-h and ODB-z silicene are at the K point and between Y and Gamma points respectively, possessing high Fermi velocities of 3.1 x 10(5) m s(-1) (ODB-h) and 2.9-3.4 x 10(5) m s(-1) (ODB-z). The origin of the Dirac cones is further explained by tight-binding models. The semimetallic properties of ODB-h and ODB-z are sensitive to compression due to the self-absorption effect, but quite robust against the tensile strain. These outstanding properties make oxidized dumbbell silicene a promising material for quantum computing and high-speed electronic devices.

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: Motivated by the successful preparation of MoSi2N4 monolayers in the last year [Y.-L. Hong et al., Science, 2020, 369, 670-674], we investigate the structural, electronic and optical properties of the MoS2/MoSi2N4 heterostructure (HTS). The phonon dispersion and the binding energy calculations refer to the stability of the HTS. The heterostructure has an indirect bandgap of 1.26 (1.84) eV using PBE (HSE06) which is smaller than the corresponding value of MoSi2N4 and MoS2 monolayers. We find that the work function of the MoS2/MoSi2N4 HTS is smaller than the corresponding value of its individual monolayers. The heterostructure structure can enhance the absorption of light spectra not only in the ultraviolet region but also in the visible region as compared to MoSi2N4 and MoS2 monolayers. The refractive index behaviour of the HTS can be described as the cumulative effect which is well described in terms of a combination of the individual effects (the refractive index of MoSi2N4 and MoS2 monolayers).

Abstract: Controlled transport of water molecules through membranes and capillaries is important in areas as diverse as water purification and healthcare technologies(1-7). Previous attempts to control water permeation through membranes (mainly polymeric ones) have concentrated on modulating the structure of the membrane and the physicochemical properties of its surface by varying the pH, temperature or ionic strength(3,8). Electrical control over water transport is an attractive alternative; however, theory and simulations(9-14) have often yielded conflicting results, from freezing of water molecules to melting of ice(14-16) under an applied electric field. Here we report electrically controlled water permeation through micrometre-thick graphene oxide membranes(17-21). Such membranes have previously been shown to exhibit ultrafast permeation of water(17,22) and molecular sieving properties(18,21), with the potential for industrial-scale production. To achieve electrical control over water permeation, we create conductive filaments in the graphene oxide membranes via controllable electrical breakdown. The electric field that concentrates around these current-carrying filaments ionizes water molecules inside graphene capillaries within the graphene oxide membranes, which impedes water transport. We thus demonstrate precise control of water permeation, from ultrafast permeation to complete blocking. Our work opens up an avenue for developing smart membrane technologies for artificial biological systems, tissue engineering and filtration.

Abstract: Two-dimensional atomic crystals can radically change their properties in response to external influences, such as substrate orientation or strain, forming materials with novel electronic structure(1-5). An example is the creation of weakly dispersive, 'flat' bands in bilayer graphene for certain 'magic' angles of twist between the orientations of the two layers(6). The quenched kinetic energy in these flat bands promotes electron-electron interactions and facilitates the emergence of strongly correlated phases, such as superconductivity and correlated insulators. However, the very accurate fine-tuning required to obtain the magic angle in twisted-bilayer graphene poses challenges to fabrication and scalability. Here we present an alternative route to creating flat bands that does not involve fine-tuning. Using scanning tunnelling microscopy and spectroscopy, together with numerical simulations, we demonstrate that graphene monolayers placed on an atomically flat substrate can be forced to undergo a buckling transition(7-9), resulting in a periodically modulated pseudo-magnetic field(10-14), which in turn creates a 'post-graphene' material with flat electronic bands. When we introduce the Fermi level into these flat bands using electrostatic doping, we observe a pseudogap-like depletion in the density of states, which signals the emergence of a correlated state(15-17). This buckling of two-dimensional crystals offers a strategy for creating other superlattice systems and, in particular, for exploring interaction phenomena characteristic of flat bands. Buckled monolayer graphene superlattices are found to provide an alternative to twisted bilayer graphene for the study of flat bands and correlated states in a carbon-based material.

Abstract: Depending on the Ginzburg-Landau parameter kappa, superconductors can either be fully diamagnetic if kappa < 1/root 2 (type I superconductors) or allow magnetic flux to penetrate through Abrikosov vortices if kappa > 1/root 2 (type II superconductors; refs 1,2). At the Bogomolny critical point, kappa = kappa(c) = 1/root 2, a state that is infinitely degenerate with respect to vortex spatial configurations arises(3,4). Despite in-depth investigations of conventional type I and type II superconductors, a thorough understanding of the magnetic behaviour in the near-Bogomolny critical regime at kappa similar to kappa(c) remains lacking. Here we report that in confined systems the critical regime expands over a finite interval of kappa forming a critical superconducting state. We show that in this state, in a sample with dimensions comparable to the vortex core size, vortices merge into a multi-quanta droplet, which undergoes Rayleigh instability(5) on increasing kappa and decays by emitting single vortices. Superconducting vortices realize Nielsen-Olesen singular solutions of the Abelian Higgs model, which is pervasive in phenomena ranging from quantum electrodynamics to cosmology(6-9). Our study of the transient dynamics of Abrikosov-Nielsen-Olesen vortices in systems with boundaries promises access to non-trivial effects in quantum field theory by means of bench-top laboratory experiments.

Abstract: Graphene’s remarkable electronic properties have fuelled the vision of a graphene-based platform for lighter, faster and smarter electronics and computing applications. One of the challenges is to devise ways to tailor graphene’s electronic properties and to control its charge carriers. Here we show that a single-atom vacancy in graphene can stably host a local charge and that this charge can be gradually built up by applying voltage pulses with the tip of a scanning tunnelling microscope. The response of the conduction electrons in graphene to the local charge is monitored with scanning tunnelling and Landau level spectroscopy, and compared to numerical simulations. As the charge is increased, its interaction with the conduction electrons undergoes a transition into a supercritical regime where itinerant electrons are trapped in a sequence of quasi-bound states which resemble an artificial atom. The quasi-bound electron states are detected by a strong enhancement of the density of states within a disc centred on the vacancy site which is surrounded by halo of hole states. We further show that the quasi-bound states at the vacancy site are gate tunable and that the trapping mechanism can be turned on and off, providing a mechanism to control and guide electrons in graphene.