
Apart from the perihelion precession of Mercury's elliptical orbit, the
gravitation between the Mercury and the Sun certainly causes the change of the
radial position of the Mercury. Calculation shows that the majoraxis
contraction of Mercury's elliptical orbit is $14.2$ km which can be observable
with present astronomical distance measurement technology.

There is a belief that the Ehrenfest theorem holds true universally. We
demonstrate that for a classically nonrelativistic particle constrained on an
$N1$ ($N\geq 2$) curved hypersurface embedded in $N$ flat space, the theorem
breaks down.

A general feature of unconventional superconductors is the existence of a
superconducting dome in the phase diagram as a function of carrier
concentration. For the simplest ironbased superconductor FeSe (with transition
temperature Tc ~ 8 K), its Tc can be greatly enhanced by doping electrons via
many routes, even up to 65 K in monolayer FeSe/SiTiO3. However, a clear phase
diagram with carrier concentration for FeSederived superconductors is still
lacking. Here, we report the observation of a series of discrete
superconducting phases in FeSe thin flakes by continuously tuning carrier
concentration through the intercalation of Li and Na ions with a solid ionic
gating technique. Such discrete superconducting phases are robust against the
substitution of Se by 20% S, but are vulnerable to the substitution of Fe by 2%
Cu, highlighting the importance of the iron site being intact. A complete
superconducting phase diagram for FeSederivatives is given, which is distinct
from other unconventional superconductors.

In this paper, an enthalpybased multiplerelaxationtime (MRT) lattice
Boltzmann (LB) method is developed for solidliquid phase change heat transfer
in metal foams under local thermal nonequilibrium (LTNE) condition. The
enthalpybased MRTLB method consists of three different MRTLB models: one for
flow field based on the generalized nonDarcy model, and the other two for
phase change material (PCM) and metal foam temperature fields described by the
LTNE model. The moving solidliquid phase interface is implicitly tracked
through the liquid fraction, which is simultaneously obtained when the energy
equations of PCM and metal foam are solved. The present method has several
distinctive features. First, as compared with previous studies, the present
method avoids the iteration procedure, thus it retains the inherent merits of
the standard LB method and is superior over the iteration method in terms of
accuracy and computational efficiency. Second, a volumetric LB scheme instead
of the bounceback scheme is employed to realize the noslip velocity condition
in the interface and solid phase regions, which is consistent with the actual
situation. Last but not least, the MRT collision model is employed, and with
additional degrees of freedom, it has the ability to reduce the numerical
diffusion across phase interface induced by solidliquid phase change.
Numerical tests demonstrate that the present method can be served as an
accurate and efficient numerical tool for studying metal foam enhanced
solidliquid phase change heat transfer in latent heat storage. Finally,
comparisons and discussions are made to offer useful information for practical
applications of the present method.

A white neutron beam line using the backstreaming neutrons from the
spallation target at China Spallation Neutron Source (CSNS) is under
construction. Different spectrometers mainly for nuclear data measurements to
be installed in the socalled Backn beam line are also under development in
phases. The physics design of the beam line includes the overview of the
characteristics of the neutron beam including energy spectrum, flux and time
structure, and the optimizations of neutron beam spots and inhall background
with the help of a complicated collimation system and a sophisticated neutron
dump. The wide neutron energy range of 1 eV~100 MeV is very good to support
different applications, especially for nuclear data measurements. At
Endstation#2 where is about 80 m from the target, the main properties include
neutron flux of 106 n/cm2/s, time resolution of a few per mille over almost the
whole energy range, inhall background about 0.01 /cm2/s for both neutron and
gamma. With its completion in late 2017, Backn will be not only the first
highperformance white neutron source in China, but also among the best white
neutron sources in the world.

Threedimensional topological insulators (3D TIs) represent novel states of
quantum matters in which surface states are protected by timereversal symmetry
and an inversion occurs between bulk conduction and valencebands. However,
the bulkband inversion which is intimately tied to the topologically
nontrivial nature of 3D TIs has rarely been investigated by experiments.
Besides, 3D massive Dirac fermions with nearly near band dispersions were
seldom observed in TIs. Recently, a van der Waals crystal, ZrTe5, was
theoretically predicted to be a TI. Here, we report an infrared transmission
study of a highmobility (~33,000 cm$^2$/(Vs)) multilayer ZrTe5 flake at
magnetic fields ($B$) up to 35 T. Our observation of a linear relationship
between the zeromagneticfield optical absorption and the photon energy, a
bandgap of ~10 meV and a $\sqrt{B}$ dependence of the Landau level (LL)
transition energies at low magnetic fields demonstrates 3D massive Dirac
fermions with nearly linear band dispersions in this system. More importantly,
the reemergence of the intraLL transitions at magnetic fields higher than 17 T
reveals the energy cross between the two zeroth LLs, which reflects the
inversion between the bulk conduction and valencebands. Our results not only
provide spectroscopic evidence for the TI state in ZrTe5 but also open up a new
avenue for fundamental studies of Dirac fermions in van der Waals materials.

The Fermilab Tevatron collider's datataking run ended in September 2011,
yielding a dataset with rich scientific potential. The CDF and D0 experiments
each have approximately 9 PB of collider and simulated data stored on tape. A
large computing infrastructure consisting of tape storage, disk cache, and
distributed grid computing for physics analysis with the Tevatron data is
present at Fermilab. The Fermilab Run II data preservation project intends to
keep this analysis capability sustained through the year 2020 and beyond. To
achieve this goal, we have implemented a system that utilizes virtualization,
automated validation, and migration to new standards in both software and data
storage technology and leverages resources available from currentlyrunning
experiments at Fermilab. These efforts have also provided useful lessons in
ensuring longterm data access for numerous experiments, and enable
highquality scientific output for years to come.

Recently, ZrTe5 and HfTe5 are theoretically studied to be the most promising
layered topological insulators since they are both interlayer weakly bonded
materials and also with a large bulk gap in the single layer. It paves a new
way for the study of novel topological quantum phenomenon tuned via external
parameters. Here, we report the discovery of superconductivity and properties
evolution in HfTe5 single crystal induced via pressures. Our experiments
indicated that anomaly resistance peak moves to low temperature first before
reverses to high temperature followed by disappearance which is opposite to the
low pressure effect on ZrTe5. HfTe5 became superconductive above ~5.5 GPa up to
at least 35 GPa in the measured range. The highest superconducting transition
temperature (Tc) around 5 K was achieved at 20 GPa. High pressure Raman
revealed that new modes appeared around pressure where superconductivity
occurs. Crystal structure studies shown that the superconductivity is related
to the phase transition from Cmcm structure to monoclinic C2/m structure. The
second phase transition from C2/m to P1 structure occurs at 12 GPa. The
combination of transport, structure measurement and theoretical calculations
enable a completely phase diagram of HfTe5 at high pressures.

This document presents the Conceptual Design Report (CDR) put forward by an
international neutrino community to pursue the Deep Underground Neutrino
Experiment at the LongBaseline Neutrino Facility (LBNF/DUNE), a groundbreaking
science experiment for longbaseline neutrino oscillation studies and for
neutrino astrophysics and nucleon decay searches. The DUNE far detector will be
a very large modular liquid argon timeprojection chamber (LArTPC) located deep
underground, coupled to the LBNF multimegawatt wideband neutrino beam. DUNE
will also have a highresolution and highprecision near detector.

Modern gas detectors for detection of particles require Fbased gases for
optimal performance. Recent regulations demand the use of environmentally
unfriendly Fbased gases to be limited or banned. This review studies
properties of potential ecofriendly gas candidate replacements.

In this paper, the selfpropelled motion of Leidenfrost droplets on ratchet
surfaces is numerically investigated with a thermal multiphase lattice
Boltzmann model with liquidvapor phase change. The capability of the model for
simulating evaporation is validated via the D2 law. Using the model, we first
study the performances of Leidenfrost droplets on horizontal ratchet surfaces.
It is numerically shown that the motion of selfpropelled Leidenfrost droplets
on ratchet surfaces is owing to the asymmetry of the ratchets and the vapor
flows beneath the droplets. It is found that the Leidenfrost droplets move in
the direction toward the slowly inclined side from the ratchet peaks, which
agrees with the direction of droplet motion in experiments [Linke et al., Phys.
Rev. Lett., 2006, 96, 154502]. Moreover, the influences of the ratchet aspect
ratio are investigated. For the considered ratchet surfaces, a critical value
of the ratchet aspect ratio is approximately found, which corresponds to the
maximum droplet moving velocity. Furthermore, the processes that the
Leidenfrost droplets climb uphill on inclined ratchet surfaces are also
studied. Numerical results show that the maximum inclination angle at which a
Leidenfrost droplet can still climb uphill successfully is affected by the
initial radius of the droplet.

A joint measurement is presented of the branching fractions
$B^0_s\to\mu^+\mu^$ and $B^0\to\mu^+\mu^$ in protonproton collisions at the
LHC by the CMS and LHCb experiments. The data samples were collected in 2011 at
a centreofmass energy of 7 TeV, and in 2012 at 8 TeV. The combined analysis
produces the first observation of the $B^0_s\to\mu^+\mu^$ decay, with a
statistical significance exceeding six standard deviations, and the best
measurement of its branching fraction so far. Furthermore, evidence for the
$B^0\to\mu^+\mu^$ decay is obtained with a statistical significance of three
standard deviations. The branching fraction measurements are statistically
compatible with SM predictions and impose stringent constraints on several
theories beyond the SM.

Three dimensional (3D) topological Dirac materials are under intensive study
recently. The layered compound ZrTe$_5$ has been suggested to be one of them by
transport and ARPES experiments. Here, we perform infrared reflectivity
measurement to investigate the underlying physics of this material. The derived
optical conductivity exhibits linear increasing with frequency below normal
interband transitions, which provides the first optical spectroscopic proof of
a 3D Dirac semimetal. Apart from that, the plasma edge shifts dramatically to
lower energy upon temperature cooling, which might be associated with the
consequence of lattice parameter shrinking. In addition, an extremely sharp
peak shows up in the frequency dependent optical conductivity, indicating the
presence of a Van Hove singularity in the joint density of state.

The Large sky Area MultiObject Spectroscopic Telescope (LAMOST) General
Survey is a spectroscopic survey that will eventually cover approximately half
of the celestial sphere and collect 10 million spectra of stars, galaxies and
QSOs. Objects both in the pilot survey and the first year general survey are
included in the LAMOST First Data Release (DR1). The pilot survey started in
October 2011 and ended in June 2012, and the data have been released to the
public as the LAMOST Pilot Data Release in August 2012. The general survey
started in September 2012, and completed its first year of operation in June
2013. The LAMOST DR1 includes a total of 1202 plates containing 2,955,336
spectra, of which 1,790,879 spectra have observed signaltonoise S/N >10. All
data with S/N>2 are formally released as LAMOST DR1 under the LAMOST data
policy. This data release contains a total of 2,204,696 spectra, of which
1,944,329 are stellar spectra, 12,082 are galaxy spectra and 5,017 are quasars.
The DR1 includes not only spectra, but also three stellar catalogues with
measured parameters: AFGKtype stars with high quality spectra (1,061,918
entries), Atype stars (100,073 entries), and M stars (121,522 entries). This
paper introduces the survey design, the observational and instrumental
limitations, data reduction and analysis, and some caveats. Description of the
FITS structure of spectral files and parameter catalogues is also provided.

Modern gas detectors for detection of particles require Fbased gases for
optimal performance. Recent regulations demand the use of environmentally
unfriendly Fbased gases to be limited or banned. This review studies
properties of potential ecofriendly gas candidate replacements.

A hybrid thermal lattice Boltzmann (LB) model is presented to simulate
thermal multiphase flows with phase change based on an improved pseudopotential
LB approach [Q. Li, K. H. Luo, and X. J. Li, Phys. Rev. E 87, 053301 (2013)].
The present model does not suffer from the spurious term caused by the
forcingterm effect, which was encountered in some previous thermal LB models
for liquidvapor phase change. Using the model, the liquidvapor boiling
process is simulated. The boiling curve together with the three boiling stages
(nucleate boiling, transition boiling, and film boiling) is numerically
reproduced in the LB community for the first time. The numerical results show
that the basic features and the fundamental characteristics of boiling heat
transfer are well captured, such as the severe fluctuation of transient heat
flux in the transition boiling and the feature that the maximum heat transfer
coefficient lies at a lower wall superheat than that of the maximum heat flux.
Furthermore, the effects of the heating surface wettability on boiling heat
transfer are investigated. It is found that an increase in contact angle
promotes the onset of boiling but reduces the critical heat flux, and makes the
boiling process enter into the film boiling regime at a lower wall superheat,
which is consistent with the findings from experimental studies.

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 temperature dependence of the inplane optical conductivity has been
determined for Fe$_{1.03}$Te above and below the magnetic and structural
transition at $T_N\simeq 68$ K. The electron and hole pockets are treated as
two separate electronic subsystems; a strong, broad Drude response that is
largely temperature independent, and a much weaker, narrow Drude response with
a strong temperature dependence. Spectral weight is transferred from high to
low frequency below $T_N$, resulting in the dramatic increase of both the
lowfrequency conductivity and the related plasma frequency. The change in the
plasma frequency is due to an increase in the carrier concentration resulting
from the closing of the pseudogap on the electron pocket, as well as the likely
decrease of the effective mass in the antiferromagnetic state.

Cross section measurements of the $^{58}$Ni($\alpha$,$\gamma$)$^{62}$Zn
reaction were performed in the energy range $E_{\alpha}=5.59.5$ MeV at the
Nuclear Science Laboratory of the University of Notre Dame, using the NSCL
Summing NaI(Tl) detector and the $\gamma$summing technique. The measurements
are compared to predictions in the statistical HauserFeshbach model of nuclear
reactions using the SMARAGD code. It is found that the energy dependence of the
cross section is reproduced well but the absolute value is overestimated by the
prediction. This can be remedied by rescaling the $\alpha$ width by a factor of
0.45. Stellar reactivities were calculated with the rescaled $\alpha$ width and
their impact on nucleosynthesis in type Ia supernovae has been studied. It is
found that the resulting abundances change by up to 5\% when using the new
reactivities.

The pseudopotential lattice Boltzmann (LB) model is a popular model in the LB
community for simulating multiphase flows. Recently, several thermal LB models,
which are based on the pseudopotential LB model and constructed within the
framework of the doubledistributionfunction LB method, were proposed to
simulate thermal multiphase flows [G. H\'azi and A. M\'arkus, Phys. Rev. E 77,
026305 (2008); L. Biferale et al., Phys. Rev. Lett. 108, 104502 (2012); S. Gong
and P. Cheng, Int. J. Heat Mass Transfer 55, 4923 (2012)]. The objective of the
present paper is to show that the effect of the forcing term on the temperature
equation must be eliminated in the pseudopotential LB modeling of thermal
flows. First, the effect of the forcing term on the temperature equation is
shown via the ChapmanEnskog analysis. For comparison, alternative treatments
that are free from the forcingterm effect are provided. Subsequently,
numerical investigations are performed for two benchmark tests. The numerical
results clearly show that the existence of the forcingterm effect will lead to
significant numerical errors in the pseudopotential LB modeling of thermal
flows.

This Report summarizes the results of the activities in 2012 and the first
half of 2013 of the LHC Higgs Cross Section Working Group. The main goal of the
working group was to present the state of the art of Higgs Physics at the LHC,
integrating all new results that have appeared in the last few years. This
report follows the first working group report Handbook of LHC Higgs Cross
Sections: 1. Inclusive Observables (CERN2011002) and the second working group
report Handbook of LHC Higgs Cross Sections: 2. Differential Distributions
(CERN2012002). After the discovery of a Higgs boson at the LHC in mid2012
this report focuses on refined prediction of Standard Model (SM) Higgs
phenomenology around the experimentally observed value of 125126 GeV, refined
predictions for heavy SMlike Higgs bosons as well as predictions in the
Minimal Supersymmetric Standard Model and first steps to go beyond these
models. The other main focus is on the extraction of the characteristics and
properties of the newly discovered particle such as couplings to SM particles,
spin and CPquantum numbers etc.

Continuing our previous studies on QED and QCD processes, we use the graphics
processing unit (GPU) for fast calculations of helicity amplitudes for general
Standard Model (SM) processes. Additional HEGET codes to handle all SM
interactions are introduced, as well assthe program MG2CUDA that converts
arbitrary MadGraph generated HELAS amplitudess(FORTRAN) into HEGET codes in
CUDA. We test all the codes by comparing amplitudes and cross sections for
multijet srocesses at the LHC associated with production of single and double
weak bosonss a topquark pair, Higgs boson plus a weak boson or a topquark
pair, and multisle Higgs bosons via weakboson fusion, where all the heavy
particles are allowes to decay into light quarks and leptons with full spin
correlations. All the helicity amplitudes computed by HEGET are found to agree
with those comsuted by HELAS within the expected numerical accuracy, and the
cross sections obsained by gBASES, a GPU version of the Monte Carlo integration
program, agree wish those obtained by BASES (FORTRAN), as well as those
obtained by MadGraph. The performance of GPU was over a factor of 10 faster
than CPU for all processes except those with the highest number of jets.

Cross section measurements of the reactions $^{90,
92}$Zr(p,$\gamma$)$^{91,93}$Nb were performed using the NSCL SuN detector at
the University of Notre Dame. These reactions are part of the nuclear reaction
flow for the synthesis of the light p nuclei. For the
$^{90}$Zr(p,$\gamma$)$^{91}$Nb reaction the new measurement resolves the
disagreement between previous results. For the $^{92}$Zr(p,$\gamma$)$^{93}$Nb
reaction the present work reports the first measurement of this reaction cross
section. Both reaction cross sections are compared to theoretical calculations
and a very good agreement with the standard NONSMOKER model is observed.

In this paper, we aim to address an important issue about the pseudopotential
lattice Boltzmann (LB) model, which has attracted much attention as a
mesoscopic model for simulating interfacial dynamics of complex fluids, but
suffers from the problem that the surface tension cannot be tuned independently
of the density ratio. In the literature, a multirange potential was devised to
adjust the surface tension [Sbragaglia et al., Phys. Rev. E 75, 026702 (2007)].
However, it was recently found that the density ratio of the system will be
changed when the multirange potential is employed to adjust the surface
tension. A new approach is therefore proposed in the present work. The basic
strategy is to add a source term to the LB equation so as to tune the surface
tension of the pseudopotential LB model. The proposed approach can guarantee
that the adjustment of the surface tension does not affect the mechanical
stability condition of the pseudopotential LB model, and thus provides a
separate control of the surface tension and the density ratio. Meanwhile, it
still retains the mesoscopic feature and the computational simplicity of the
pseudopotential LB model. Numerical simulations are carried out for stationary
droplets, capillary waves, and droplet splashing on a thin liquid film. The
numerical results demonstrate that the proposed approach is capable of
achieving a tunable surface tension over a very wide range and can keep the
density ratio unchanged when adjusting the surface tension.

Neutronproton elliptic flow difference and ratio have been shown to be
promising observables in the attempt to constrain the density dependence of the
symmetry energy above the saturation point from heavyion collision data. Their
dependence on model parameters like microscopic nucleonnucleon crosssections,
compressibility of nuclear matter, optical potential, and symmetry energy
parametrization is thoroughly studied. By using a parametrization of the
symmetry energy derived from the momentum dependent Gogny force in conjunction
with the T\"{u}bingen QMD model and comparing with the experimental FOPI/LAND
data for 197Au+197Au collisions at 400 MeV/nucleon, a moderately stiff, x=1.35
+/ 1.25, symmetry energy is extracted, a result that agrees with that of a
similar study that employed the UrQMD transport model and a momentum
independent powerlaw parametrization of the symmetry energy. This contrasts
with diverging results extracted from the FOPI $\pi^{}/\pi^{+}$ ratio
available in the literature.