
Using the tightbinding approach, we investigate the energy spectrum of
square, triangular and hexagonal MoS$_2$ quantum dots (QDs) in the presence of
a perpendicular magnetic field. Novel edge states emerge in MoS$_2$ QDs, which
are distributed over the whole edge which we call ring states. The ring states
are robust in the presence of spinorbit coupling (SOC). The corresponding
energy levels of the ring states oscillate as function of the perpendicular
magnetic field which are related to AharonovBohm oscillations. Oscillations in
the magnetic field dependence of the energy levels and the peaks in the
magnetooptical spectrum emerge (disappear) as the ring states are formed
(collapsed). The period and the amplitude of the oscillation decreases with the
size of the MoS$_2$ QDs.

Efficient polarization of organic molecules is of extraordinary relevance
when performing nuclear magnetic resonance (NMR) and imaging. Commercially
available routes to dynamical nuclear polarization (DNP) work at extremely
lowtemperatures, thus bringing the molecules out of their ambient thermal
conditions and relying on the solidification of organic samples. In this work
we investigate polarization transfer from opticallypumped nitrogen vacancy
centers in diamond to external molecules at room temperature. This polarization
transfer is described by both an extensive analytical analysis and numerical
simulations based on spin bath bosonization and is supported by experimental
data in excellent agreement. These results set the route to hyperpolarization
of diffusive molecules in different scenarios and consequently, due to
increased signal, to highresolution NMR.

We systematically investigate the magnetic properties and local structure of
Ba2YIrO6 to demonstrate that Y and Ir lattice defects in the form of antiphase
boundary or clusters of antisite disorder affect the magnetism observed in this
$d^4$ compound. We compare the magnetic properties and atomic imaging of (1) a
slow cooled crystal, (2) a crystal quenched from 900\degree C after growth, and
(3) a crystal grown using a faster cooling rate than the slow cooled one.
Atomic imaging by scanning transmission electron microscopy (STEM) shows that
quenching from 900oC introduces antiphase boundary to the crystals, and a
faster cooling rate during crystal growth leads to clusters of Y and Ir
antisite disorder. STEM study suggests the antiphase boundary region is Irrich
with a composition of Ba2YIrO6. The magnetic measurements show that Ba2YIrO6
crystals with clusters of antisite defects have a larger effective moment and a
larger saturation moment than the slowcooled crystals. Quenched crystals with
Irrich antiphase boundary shows a slightly suppressed saturation moment than
the slow cooled crystals, and this seems to suggest that antiphase boundary is
detrimental to the moment formation. Our DFT calculations suggest magnetic
condensation is unlikely as the energy to be gained from superexchange is small
compared to the spinorbit gap. However, once Y is replaced by Ir in the
antisite disordered region, the picture of local nonmagnetic singlets breaks
down and magnetism can be induced. This is because of (a) enhanced interactions
due to increased overlap of orbitals between sites, and, (b) increased number
of orbitals mediating the interactions. Our work highlights the importance of
lattice defects in understanding the experimentally observed magnetism in
Ba2YIrO6 and other J=0 systems.

Coexistence of different geometric shapes at low energies presents a
universal structure phenomenon that occurs over the entire chart of nuclides.
Studies of the shape coexistence are important for understanding the
microscopic origin of collectivity and modifications of shell structure in
exotic nuclei far from stability. The aim of this work is to provide a
systematic analysis of characteristic signatures of coexisting nuclear shapes
in different mass regions, using a global selfconsistent theoretical method
based on universal energy density functionals and the quadrupole collective
model. The lowenergy excitation spectrum and quadrupole shape invariants of
the two lowest $0^{+}$ states of eveneven nuclei are obtained as solutions of
a fivedimensional collective Hamiltonian (5DCH) model, with parameters
determined by constrained selfconsistent meanfield calculations based on the
relativistic energy density functional PCPK1, and a finiterange pairing
interaction. The theoretical excitation energies of the states: $2^+_1$,
$4^+_1$, $0^+_2$, $2^+_2$, $2^+_3$, as well as the $B(E2; 0^+_1\to 2^+_1)$
values, are in very good agreement with the corresponding experimental values
for 621 eveneven nuclei. Quadrupole shape invariants have been implemented to
investigate shape coexistence, and the distribution of possible
shapecoexisting nuclei is consistent with results obtained in recent
theoretical studies and available data. The present analysis has shown that,
when based on a universal and consistent microscopic framework of nuclear
density functionals, shape invariants provide distinct indicators and reliable
predictions for the occurrence of lowenergy coexisting shapes. This method is
particularly useful for studies of shape coexistence in regions far from
stability where few data are available.

We have performed magnetic susceptibility, heat capacity, muon spin
relaxation, and neutron scattering measurements on three members of the family
Ba3MRu2O9, where M = In, Y and Lu. These systems consist of mixedvalence Ru
dimers on a triangular lattice with antiferromagnetic interdimer exchange.
Although previous work has argued that charge order within the dimers or
intradimer double exchange plays an important role in determining the magnetic
properties, our results suggest that the dimers are better described as
molecular units due to significant orbital hybridization, resulting in one
spin1/2 moment distributed equally over the two Ru sites. These molecular
building blocks form a frustrated, quasitwodimensional triangular lattice.
Our zero and longitudinal field muSR results indicate that the molecular
moments develop a collective, static magnetic ground state, with oscillations
of the zero field muon spin polarization indicative of longrange magnetic
order in the Lu sample. The static magnetism is much more disordered in the Y
and In samples, but they do not appear to be conventional spin glasses.

The coherent highfidelity generation of nuclear spins in longlived singlet
states which may find application as quantum memory or sensor represents a
considerable experimental challenge. Here we propose a dissipative scheme that
achieves the preparation of pairs of nuclear spins in longlived singlet states
by a protocol that combines the interaction between the nuclei and a
periodically reset electron spin of an NV center with local rfcontrol of the
nuclear spins. The final state of this protocol is independent of the initial
preparation of the nuclei, is robust to external field fluctuations and can be
operated at room temperature. We show that a high fidelity singlet pair of a
13C dimer in a nuclear bath in diamond can be generated under realistic
experimental conditions.

We propose to use a dissipatively engineered nitrogen vacancy (NV) center as
a mediator of interaction between two nuclear spins that are protected from
decoherence and relaxation of the NV. Under ambient conditions this scheme
achieves highly selective highfidelity quantum gates between nuclear spins in
a quantum register even at large NVnuclear distances. Importantly, this method
allows for the use of nuclear spins as a sensor rather than a memory, while the
NV spin acts as an ancillary system for the initialization and read out of the
sensor. The immunity to the decoherence and relaxation of the NV center leads
to a tunable sharp frequency filter while allowing at the same time the
continuous collection of the signal to achieve simultaneously high spectral
selectivity and high signaltonoise ratio (SNR).

Over the past few decades, tremendous progress has been made in the
development of particlebased discrete simulation methods versus the
conventional continuumbased methods. In particular, the lattice Boltzmann (LB)
method has evolved from a theoretical novelty to a ubiquitous, versatile and
powerful computational methodology for both fundamental research and
engineering applications. It is a kineticbased mesoscopic approach that
bridges the microscales and macroscales, which offers distinctive advantages in
simulation fidelity and computational efficiency. Applications of the LB method
have been found in a wide range of disciplines including physics, chemistry,
materials, biomedicine and various branches of engineering. The present work
provides a comprehensive review of the LB method for thermofluids and energy
applications, focusing on multiphase flows, thermal flows and thermal
multiphase flows with phase change. The review first covers the theoretical
framework of the LB method, revealing the existing inconsistencies and defects
as well as common features of multiphase and thermal LB models. Recent
developments in improving the thermodynamic and hydrodynamic consistency,
reducing the spurious currents, enhancing the numerical stability, etc., are
highlighted. These efforts have put the LB method on a firmer theoretical
foundation with enhanced LB models that can achieve larger liquidgas density
ratio, higher Reynolds number and flexible surface tension. Examples of
applications are provided in fuel cells and batteries, droplet collision,
boiling heat transfer and evaporation, and energy storage. Finally, further
developments and future prospect of the LB method are outlined for thermofluids
and energy applications.

The Jiangmen Underground Neutrino Observatory (JUNO) is proposed to determine
the neutrino mass hierarchy using an underground liquid scintillator detector.
It is located 53 km away from both Yangjiang and Taishan Nuclear Power Plants
in Guangdong, China. The experimental hall, spanning more than 50 meters, is
under a granite mountain of over 700 m overburden. Within six years of running,
the detection of reactor antineutrinos can resolve the neutrino mass hierarchy
at a confidence level of 34$\sigma$, and determine neutrino oscillation
parameters $\sin^2\theta_{12}$, $\Delta m^2_{21}$, and $\Delta m^2_{ee}$ to
an accuracy of better than 1%. The JUNO detector can be also used to study
terrestrial and extraterrestrial neutrinos and new physics beyond the Standard
Model. The central detector contains 20,000 tons liquid scintillator with an
acrylic sphere of 35 m in diameter. $\sim$17,000 508mm diameter PMTs with high
quantum efficiency provide $\sim$75% optical coverage. The current choice of
the liquid scintillator is: linear alkyl benzene (LAB) as the solvent, plus PPO
as the scintillation fluor and a wavelengthshifter (BisMSB). The number of
detected photoelectrons per MeV is larger than 1,100 and the energy resolution
is expected to be 3% at 1 MeV. The calibration system is designed to deploy
multiple sources to cover the entire energy range of reactor antineutrinos, and
to achieve a fullvolume position coverage inside the detector. The veto system
is used for muon detection, muon induced background study and reduction. It
consists of a Water Cherenkov detector and a Top Tracker system. The readout
system, the detector control system and the offline system insure efficient and
stable data acquisition and processing.

Here we propose and analyse in detail protocols that can achieve rapid
hyperpolarization of 13C nuclear spins in randomly oriented ensembles of
nanodiamonds at room temperature. Our protocols exploit a combination of
optical polarization of electron spins in nitrogenvacancy centers and the
transfer of this polarization to 13C nuclei by means of microwave control to
overcome the severe challenges that are posed by the random orientation of the
nanodiamonds and their nitrogenvacancy centers. Specifically, these random
orientations result in exceedingly large energy variations of the electron spin
levels that render the polarization and coherent control of the
nitrogenvacancy center electron spins as well as the control of their coherent
interaction with the surrounding 13C nuclear spins highly inefficient. We
address these challenges by a combination of an offresonant microwave double
resonance scheme in conjunction with a realization of the integrated solid
effect which, together with adiabatic rotations of external magnetic fields or
rotations of nanodiamonds, leads to a protocol that achieves high levels of
hyperpolarization of the entire nuclearspin bath in a randomly oriented
ensemble of nanodiamonds even at room temperature. This hyperpolarization
together with the long nuclear spin polarization lifetimes in nanodiamonds and
the relatively high density of 13C nuclei has the potential to result in a
major signal enhancement in 13C nuclear magnetic resonance imaging and suggests
functionalized and hyperpolarized nanodiamonds as a unique probe for molecular
imaging both in vitro and in vivo.

During its first six years of operation, the Fermi Large Area Telescope (LAT)
has detected >30 MeV gammaray emission from more than 40 solar flares, nearly
a factor of 10 more than those detected by EGRET. These include detections of
impulsive and sustained emissions, extending up to 20 hours in the case of the
2012 March 7 Xclass flares. We will present an overview of solar flare
detections with LAT, highlighting recent results and surprising features,
including the detection of >100 MeV emission associated with flares located
behind the limb. Such flares may shed new light on the relationship between the
sites of particle acceleration and gammaray emission.

In this paper, we aim to investigate the implementation of contact angles in
the pseudopotential lattice Boltzmann modeling of wetting at a large density
ratio. The pseudopotential lattice Boltzmann model [X. Shan and H. Chen, Phys.
Rev. E 49, 2941 (1994)] is a popular mesoscopic model for simulating multiphase
flows and interfacial dynamics. In this model, the contact angle is usually
realized by a fluidsolid interaction. Two widely used fluidsolid
interactions: the densitybased interaction and the pseudopotentialbased
interaction, as well as a modified pseudopotentialbased interaction formulated
in the present paper, are numerically investigated and compared in terms of the
achievable contact angles, the maximum and the minimum densities, and the
spurious currents. It is found that the pseudopotentialbased interaction works
well for simulating small static (liquid) contact angles, however, is unable to
reproduce static contact angles close to 180 degrees. Meanwhile, it is found
that the proposed modified pseudopotentialbased interaction performs better in
light of the maximum and the minimum densities and is overall more suitable for
simulating large contact angles as compared with the other two types of
fluidsolid interactions. Furthermore, the spurious currents are found to be
enlarged when the fluidsolid interaction force is introduced. Increasing the
kinematic viscosity ratio between the vapor and liquid phases is shown to be
capable of reducing the spurious currents caused by the fluidsolid
interactions.

The bioinertness of titanium and its alloys attracts their use as bone
implants. However a bioactive coating is usually necessary for improving the
bone bonding of such implants. In this study, electrophoretic deposition(EPD)
of hydroxyapatite (HA) powder on titanium plate was performed using butanol as
solvent under direct current (DC) and alternating current (AC) fields. The zeta
potential of the suspensions was measured to understand their stability and the
charge on the particles. Coating thickness was varied by adjusting the voltage
and time of deposition. Surface morphology and cross section thickness were
studied using scanning electron microscopy and image analysis software. Surface
crack density was calculated from the micrographs. The results showed that the
samples of similar thickness have higher grain density when coated using AC as
compared to DC EPD. This facile but novel test proves the capability of ACEPD
to attain denser and uniform HA coatings from nonaqueous medium.

In this paper we present a new multiscale simulation scheme for
nextgeneration electronic design automation for nanoelectronics. The scheme
features a combination of the firstprinciples quantum mechanical calculation,
semiclassical semiconductor device simulation, compact model generation and
circuit simulation. To demonstrate the feasibility of the proposed scheme, we
apply our newly developed quantum mechanics/electromagnetics method to simulate
the junctionless transistors. The simulation results are consistent with the
experimental measurements and provide new insights on the depletion effect of
the heterodoped gate on the drain current. Based on the calculated IV curves,
a compact model is then constructed for the junctionless transistors. The
validity of the compact model is further verified by the transient circuit
simulation of an inverter.

We describe the concepts and technical realization of the highresolution
softXray beamline ADRESS operating in the energy range from 300 to 1600 eV
and intended for Resonant Inelastic Xray Scattering (RIXS) and AngleResolved
Photoelectron Spectroscopy (ARPES). The photon source is an undulator of novel
fixedgap design where longitudinal movement of permanent magnetic arrays
controls not only the light polarization (including circular and 0180 deg
rotatable linear polarizations) but also the energy without changing the gap.
The beamline optics is based on the wellestablished scheme of plane grating
monochromator (PGM) operating in collimated light. The ultimate resolving power
E/dE is above 33000 at 1 keV photon energy. The choice of blazed vs lamellar
gratings and optimization of their profile parameters is described. Due to
glancing angles on the mirrors as well as optimized groove densities and
profiles of the gratings, high photon flux is achieved up to 1.0e13
photons/s/0.01%BW at 1 keV. Ellipsoidal refocusing optics used for the RIXS
endstation demagnifies the vertical spot size down to 4 um, which allows
slitless operation and thus maximal transmission of the highresolution RIXS
spectrometer delivering E/dE better than 11000 at 1 keV photon energy. Apart
from the beamline optics, we give an overview of the control system, describe
diagnostics and software tools, and discuss strategies used for the optical
alignment. An introduction to the concepts and instrumental realization of the
ARPES and RIXS endstations is given.

Direct numerical simulations of threedimensional (3D) homogeneous turbulence
under rapid rigid rotation are conducted to examine the predictions of resonant
wave theory for both small Rossby number and large Reynolds number. The
simulation results reveal that there is a clear inverse energy cascade to the
large scales, as predicted by 2D NavierStokes equations for resonant
interactions of slow modes. As the rotation rate increases, the
verticallyaveraged horizontal velocity field from 3D NavierStokes converges
to the velocity field from 2D NavierStokes, as measured by the energy in their
difference field. Likewise, the verticallyaveraged vertical velocity from 3D
NavierStokes converges to a solution of the 2D passive scalar equation. The
energy flux directly into small wave numbers in the $k_z=0$ plane from
nonresonant interactions decreases, while fastmode energy concentrates closer
to that plane. The simulations are consistent with an increasingly dominant
role of resonant triads for more rapid rotation.

We show that Kolmogorov multipliers in turbulence cannot be statistically
independent of others at adjacent scales (or even a finite range apart) by
numerical simulation of a shell model and by theory. As the simplest
generalization of independent distributions, we suppose that the steadystate
statistics of multipliers in the shell model are given by a
translationinvariant Gibbs measure with a shortrange potential, when
expressed in terms of suitable ``spin'' variables: realvalued spins that are
logarithms of multipliers and XYspins defined by local dynamical phases.
Numerical evidence is presented in favor of the hypothesis for the shell model,
in particular novel scaling laws and derivative relations predicted by the
existence of a thermodynamic limit. The Gibbs measure appears to be in a
hightemperature, uniquephase regime with ``paramagnetic'' spin order.

Experimental results are reported for the bulk motion induced in a bed of
granular matter contained in a cylindrical pan with a flat bottom subjected to
simultaneous vertical and horizontal vibrations. The motion in space of the
moving pan is quantified. A number of distinct bulk dynamical modes are
observed in which the particle bed adopts different shapes and motions. At the
lowest pan excitation frequency $\omega$, the bed forms a ``heap,'' and rotates
about the cylinder axis. As $\omega$ is increased, a more complex ``toroidal''
mode appears in which the bed takes the shape of a torus; in this mode,
circulation occurs both about the cylinder axis, and also radially, with
particles moving from the outer edge of the pan to the centre on the top
surface of the bed, and back to the outer edge along the pan bottom. At the
highest $\omega$, surface modulations (``surface waves'' and ``sectors'') of
the toroidal mode occur. The origin of this family of behavior in terms of the
pan motion is discussed.

The structure of the Al_{70}Pd_{21}Mn_{9} surface has been investigated using
high resolution scanning tunnelling microscopy (STM).
From two large fivefold terraces on the surface in a short decorated
Fibonacci sequence, atomically resolved surface images have been obtained. One
of these terraces carries a rare local configuration in a form of a ring.
The location of the corresponding sequence of terminations in the bulk model
M of icosahedral iAlPdMn based on the threedimensional tiling T*(2F) of an
Fphase has been estimated using this ring configuration and the requirement
from the LEED work of Gierer et al. that the average atomic density of the
terminations is 0.136 atoms per A^2.
A termination contains two atomic plane layers separated by a vertical
distance of 0.48 A. The position of the bulk terminations is fixed within the
layers of Bergman polytopes in the model M: they are 4.08 A in the direction of
the bulk from a surface of the most dense Bergman layers.
From the coding windows of the top planes in terminations in M we conclude
that a Penrose (P1) tiling is possible on almost all fivefold terraces. The
shortest edge of the tiling P1, is either 4.8 A or 7.8 A.
The experimentally derived tiling of the surface with the ring configuration
has an edgelength of 8.0 + 0.3 A and hence matches the minimal edgelength
expected from the model.

The present paper reports highaccuracy crosssection data for the 2H(n,nnp)
reaction in the neutronproton (np) and neutronneutron (nn)
finalstateinteraction (FSI) regions at an incident mean neutron energy of
13.0 MeV. These data were analyzed with rigorous threenucleon calculations to
determine the 1S0 np and nn scattering lengths, a_np and a_nn. Our results are
a_nn = 18.7 +/ 0.6 fm and a_np = 23.5 +/ 0.8 fm. Since our value for a_np
obtained from neutrondeuteron (nd) breakup agrees with that from free np
scattering, we conclude that our investigation of the nn FSI done
simultaneously and under identical conditions gives the correct value for a_nn.
Our value for a_nn is in agreement with that obtained in piondeuteron capture
measurements but disagrees with values obtained from earlier nd breakup
studies.