
Both multistreaming (random motion) and bulk motion cause the FingerofGod
(FoG) effect in redshift space distortion (RSD). We apply a direct measurement
of the multistreaming effect in RSD from simulations, proving that it induces
an additional, nonnegligible FoG damping to the redshift space density power
spectrum. We show that, including the multistreaming effect, the RSD modelling
is significantly improved. We also provide a theoretical explanation based on
halo model for the measured effect, including a fitting formula with one to two
free parameters. The improved understanding of FoG helps break the
$f\sigma_8\sigma_v$ degeneracy in RSD cosmology, and has the potential of
significantly improving cosmological constraints.

The superconductor to insulator transition(SIT) has been intensively studied
in thinfilm superconductors for a few decades by tuning their
thickness,disorder or external magnetic field.It is now established that, when
the normalstate conductance is high, superconductivity systematically gives
way to an anomalous quantum metallic state(AQM).In highly resistive systems,
the AQM intervenes the SI transition over a broad range of the tuning parameter
at zero temperature, while in clean systems it crosses over to Fermi liquid
metallic behavior at the upper critical field. The details of this
superconducting to quantum metal transition (SQMT) at T=0 invokes much interest
in the nature of this exotic ground state. However, the SIT/SQMT were not yet
investigated in a 2D superconductor with coexisting and fluctuating quantum
orders. Here, we report the observation of an AQM state as T approaches 0 in
iongel gated 1TTiSe2 in the 2D limit, driven by both magnetic field and
carrier doping. The gate tunability allows us to establish a threeparameter
phase diagram in terms of magnetic field, temperature and carrier density.
Scaling of the resistance at low temperatures indicates a fieldinduced
crossover from a Bose quantum metal to vortex quantum creeping with increasing
field. As it is well known that superconducting order coexists in this crystal
with a nearcommensurate charge density wave(CDW), we discuss the interplay
between superconducting and CDW fluctuations(discommensurations) and their
potential role in enabling the AQM phase. From our findings, gatetunable
1TTiSe2 emerges as a privileged platform to scrutinize, in a controlled way,
the details of the SQMT, the role of coexisting fluctuating orders and,
ultimately, obtain a deeper understanding of the fate of superconductivity in
strictly twodimensional crystals near zero temperature.

A long standing issue in peculiar velocity cosmology is whether the
halo/galaxy velocity bias $b_v=1$ at large scale. The resolution of this
important issue must resort to high precision cosmological simulations.
However, this is hampered by another long standing "sampling artifact" problem
in volume weighted velocity measurement. We circumvent this problem with a
hybrid approach. We first measure statistics free of sampling artifact, then
link them to volume weighted statistics in theory, and finally solve for the
velocity bias. $b_v$ determined by our method is not only free of sampling
artifact, but also free of cosmic variance. We apply this method to a
$\Lambda$CDM Nbody simulation of $3072^3$ simulation particles and $1200 {\rm
Mpc}/h$ box size. For the first time we determine the halo velocity bias of
various mass and redshift to $0.1\%$$1\%$ accuracy. Major findings are as
follows. (1) $b_v\neq 1$ at $k>0.1h/{\rm Mpc}$. The deviation from unity
($b_v1$) increases with $k$. Depending on halo mass and redshift, it may
reach $\mathcal{O}(0.01)$ at $k=0.2h/{\rm Mpc}$ and $\mathcal{O}(0.05)$ at
$k\sim 0.3h/{\rm Mpc}$. The discovered $b_v\neq 1$ has statistically
significant impact on structure growth rate measurement by spectroscopic
redshift surveys, including DESI, Euclid and SKA. (2) Both the sign and the
amplitude of $b_v1$ depend on mass and redshift. These results disagree with
the peak model prediction on protohalos in that $b_v$ has much weaker
deviation from unity, varies with redshift, and can be bigger than unity. (3)
Most of the mass and redshift dependences can be compressed into a single
dependence on the halo density bias. Based on this finding, we provide an
approximate twoparameter fitting formula.

We study the ergodic control problem for a class of jump diffusions in
$\mathbb{R}^d$, which are controlled through the drift with bounded controls.
The Levy measure is finite, but has no particular structure; it can be
anisotropic and singular. Moreover, there is no blanket ergodicity assumption
for the controlled process. Unstable behavior is `discouraged' by the running
cost which satisfies a mild coercive hypothesis (i.e., is nearmonotone). We
first study the problem in its weak formulation as an optimization problem on
the space of infinitesimal ergodic occupation measures, and derive the
HamiltonJacobiBellman equation under minimal assumptions on the parameters,
including verification of optimality results, using only analytical arguments.
We also examine the regularity of invariant measures. Then, we address the jump
diffusion model, and obtain a complete characterization of optimality.

Graph embedding has been proven to be efficient and effective in facilitating
graph analysis. In this paper, we present a novel spectral framework called
NOnBacktracking Embedding (NOBE), which offers a new perspective that
organizes graph data at a deep level by tracking the flow traversing on the
edges with backtracking prohibited. Further, by analyzing the nonbacktracking
process, a technique called graph approximation is devised, which provides a
channel to transform the spectral decomposition on an edgetoedge matrix to
that on a nodetonode matrix. Theoretical guarantees are provided by bounding
the difference between the corresponding eigenvalues of the original graph and
its graph approximation. Extensive experiments conducted on various realworld
networks demonstrate the efficacy of our methods on both macroscopic and
microscopic levels, including clustering and structural hole spanner detection.

SnSe is a promising thermoelectric material with recordbreaking figure of
merit, \textit{i.e., ZT}. As a semiconductor, optimal electrical dosage is the
key challenge to maximize \textit{ZT} in SnSe. However, to date a comprehensive
understanding of the electronic structure and most critically, the selfhole
doping mechanism in SnSe is still absent. Here, we report the highly
anisotropic electronic structure of SnSe investigated by both angleresolved
photoemission spectroscopy and quantum transport, in which a unique
"\textit{puddingmold}" shaped valence band with quasilinear energy dispersion
is revealed. We prove that the electrical doping in SnSe is extrinsically
controlled by the formation of SnSe$_{2}$ microdomains induced by local phase
segregation. Using different growth methods and conditions, we have achieved
wide tuning of hole doping in SnSe, ranging from intrinsic semiconducting
behaviour to typical metal with carrier density of $1.23\times 10^{18}$
cm$^{3}$ at room temperature. The resulting multivalley transport in $p$SnSe
is characterized by nonsaturating weak localization along the armchair axis,
due to strong intervalley scattering enhanced by inplane ferroelectric dipole
field of the puckering lattice. Strikingly, quantum oscillations of
magnetoresistance reveal threedimensional electronic structure with unusual
interlayer coupling strength in $p$SnSe, which is correlated to the
interweaving of SnSe individual layers by unique point dislocation defects. Our
results suggest that defect engineering may provide versatile routes in
improving the thermoelectric performance of the SnSe family.

Massive spectroscopic redshift surveys open a promising window to accurately
measure peculiar velocity at cosmological distances through redshift space
distortion (RSD). In paper I of this series of work we proposed to decompose
peculiar velocity into three eigenmodes (v_\delta, v_S and v_B) in order to
facilitate the RSD modeling and peculiar velocity reconstruction. In the
current paper we measure the dark matter RSD related statistics of the velocity
eigenmodes through a set of Nbody simulations, including the velocity power
spectra, correlation functions, onepoint probability distribution functions,
cumulants and the damping functions describing the Finger of God effect. (1)
The power spectrum measurement shows that these velocity components have
distinctly different spatial distribution and redshift evolution. In
particular, we measure the window function \tilde{W}(k,z), which describes the
impact of nonlinear evolution on the v_\deltadensity relation. We confirm that
it can induce a significant systematic error of O(10%) in RSD cosmology. We
demonstrate that \tilde{W} can be accurately described by a simple fitting
formula with one or two free parameters. (2) The correlation function
measurement shows that the correlation length is O(100), O(10) and O(1) Mpc for
v_\delta, v_S and v_B respectively. These correlation lengths determine where
we can treat the velocity fields as spatially uncorrelated. (3) The velocity
PDFs and cumulants quantify nonGaussianities of the velocity fields. We
confirm speculation in paper I that v_\delta is largely Gaussian, nevertheless
with nonnegligible nonGaussianity, v_B is significantly nonGaussian. We also
measure the damping functions. Despite the observed nonGaussianities, the
damping functions and hence the FOG effect are all well approximated as
Gaussian ones at scales of interest.

In this paper we develop a new framework that captures the common landscape
underlying the common nonconvex lowrank matrix problems including matrix
sensing, matrix completion and robust PCA. In particular, we show for all above
problems (including asymmetric cases): 1) all local minima are also globally
optimal; 2) no highorder saddle points exists. These results explain why
simple algorithms such as stochastic gradient descent have global converge, and
efficiently optimize these nonconvex objective functions in practice. Our
framework connects and simplifies the existing analyses on optimization
landscapes for matrix sensing and symmetric matrix completion. The framework
naturally leads to new results for asymmetric matrix completion and robust PCA.

Chiral anomaly induced negative magnetoresistance (NMR) has been widely used
as a critical transport evidence on the existence of Weyl fermions in
topological semimetals. In this mini review, we discuss the general observation
of the NMR phenomena in noncentrosymmetric NbP and NbAs. We show that NMR can
be contributed by intrinsic chiral anomaly of Weyl fermions and/or extrinsic
effects, such as superimposition of Hall signals, fielddependent inhomogeneous
current flow in the bulk, i.e. current jetting, and weak localization (WL) of
coexistent trivial carriers. Such WL controlled NMR is heavily dependent on
sample quality, and is characterized by pronounced crossover from positive to
negative MR growth at elevated temperatures, as a result of the competition
between the phase coherence time and the spinorbital scattering constant of
the bulk trivial pockets. Thus, the correlation of NMR and chiral anomaly needs
to be scrutinized, without the support of other complimentary techniques. Due
to the lifting of spin degeneracy, the spin orientations of Weyl fermions are
either parallel or antiparallel to the momentum, a unique physical property
known as helicity. The conservation of helicity provides strong protection for
the transport of Weyl fermions, which can only be effectively scattered by
magnetic impurities. Chemical doping of magnetic and nonmagnetic impurities
are thus more convincing in probing the existence of Weyl fermions than the NMR
method.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV groundbased dark
energy experiment that will study baryon acoustic oscillations (BAO) and the
growth of structure through redshiftspace distortions with a widearea galaxy
and quasar redshift survey. To trace the underlying dark matter distribution,
spectroscopic targets will be selected in four classes from imaging data. We
will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to
even higher redshift, DESI will target bright [O II] emission line galaxies up
to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying
dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the
Ly$\alpha$ forest absorption features in their spectra, which will be used to
trace the distribution of neutral hydrogen. When moonlight prevents efficient
observations of the faint targets of the baseline survey, DESI will conduct a
magnitudelimited Bright Galaxy Survey comprising approximately 10 million
galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy
and quasar redshifts will be obtained to measure the BAO feature and determine
the matter power spectrum, including redshift space distortions.

DESI (Dark Energy Spectropic Instrument) is a Stage IV groundbased dark
energy experiment that will study baryon acoustic oscillations and the growth
of structure through redshiftspace distortions with a widearea galaxy and
quasar redshift survey. The DESI instrument is a roboticallyactuated,
fiberfed spectrograph capable of taking up to 5,000 simultaneous spectra over
a wavelength range from 360 nm to 980 nm. The fibers feed ten threearm
spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500,
depending on wavelength. The DESI instrument will be used to conduct a
fiveyear survey designed to cover 14,000 deg$^2$. This powerful instrument
will be installed at prime focus on the 4m Mayall telescope in Kitt Peak,
Arizona, along with a new optical corrector, which will provide a threedegree
diameter field of view. The DESI collaboration will also deliver a
spectroscopic pipeline and data management system to reduce and archive all
data for eventual public use.

The mapping of dark matter clustering from real space to redshift space
introduces the anisotropic property to the measured density power spectrum in
redshift space, known as the redshift space distortion effect. The mapping
formula is intrinsically nonlinear, which is complicated by the higher order
polynomials due to indefinite cross correlations between the density and
velocity fields, and the FingerofGod effect due to the randomness of the
peculiar velocity field. Whilst the full higher order polynomials remain
unknown, the other systematics can be controlled consistently within the same
order truncation in the expansion of the mapping formula, as shown in this
paper. The systematic due to the unknown nonlinear density and velocity fields
is removed by separately measuring all terms in the expansion directly using
simulations. The uncertainty caused by the velocity randomness is controlled by
splitting the FoG term into two pieces, 1) the "onepoint" FoG term being
independent of the separation vector between two different points, and 2) the
"correlated" FoG term appearing as an indefinite polynomials which is expanded
in the same order as all other perturbative polynomials. Using 100 realizations
of simulations, we find that the Gaussian FoG function with only one
scaleindependent free parameter works quite well, and that our new mapping
formulation accurately reproduces the observed 2dimensional density power
spectrum in redshift space at the smallest scales by far, up to $k\sim
0.2h$Mpc, considering the resolution of future experiments.

We synthesized nonmagnetic PtBi$_2$ single crystals and observed a giant
linear magnetoresistance (MR) up to 684\% under a magnetic field $\mu_0H$ = 15
T at $T$ = 2 K. The linear MR decreases with increasing temperature, but it is
still as large as 61\% under $\mu_0H$ of 15 T at room temperature. Such a giant
linear MR is unlikely to be described by the quantum model as the quantum
condition is not satisfied. Instead, we found that the slope of MR scales with
the Hall mobility, and it can be well explained by a classical disorder model.

Negative magnetoresistance (NMR) induced by the AdlerBellJackiw anomaly is
regarded as the most prominent quantum signature of Weyl semimetals when
electrical field $E$ is collinear with the external magnetic field $B$. In this
article, we report universal NMR in nonmagnetic, centrosymmetric transition
metal dipnictides MPn$_{2}$ (M=Nb and Ta; Pn=As and Sb), in which the existence
of Weyl fermions can be explicitly excluded. Using temperaturedependent
magnetoresistance, Hall and thermoelectric coefficients of Nernst and Seebeck
effects, we determine that the emergence of the NMR phenomena in MPn$_{2}$ is
coincident with a Lifshitz transition, corresponding to the formation of unique
electronholeelectron ($e$$h$$e$) pockets along the $ILI'$ direction.
Firstprinciples calculations reveal that, along the $ILI'$ line, the
$d_{xy}$ and $d_{x^{2}y^{2}}$ orbitals of the transition metal form tilted
nodal rings of band crossing well below the Fermi level. Strong spinorbital
coupling gaps all the crossing points and creates the characteristic
$e$$h$$e$ structure, making MPn$_{2}$ a topological semimetal with
$\mathbb{Z}_2$ indices of [0;(111)]. By excluding the weak localization
contribution of the bulk states, we conclude that the universal NMR in
MPn$_{2}$ may have an exotic origin in topological surface states, which
appears in pairs with opposite spinmomentum locking on nontrivial surfaces.

Noncentrosymmetric transition metal monopnictides, including TaAs, TaP,
NbAs, and NbP, are emergent topological Weyl semimetals (WSMs) hosting exotic
relativistic Weyl fermions. In this letter, we elucidate the physical origin of
the unprecedented charge carrier mobility of NbP, which can reach
$1\times10^{7}$ cm $^{2}$V$^{1}$s$^{1}$ at 1.5 K. Angle and
temperaturedependent quantum oscillations, supported by density function
theory calculations, reveal that NbP has the coexistence of p and ntype WSM
pockets in the $k_{z}$=1.16$\pi$/c plane (W1WSM) and in the $k_{z}$=0 plane
near the high symmetry points $\Sigma$ (W2WSM), respectively. Uniquely, each
W2WSM pocket forms a large dumbbellshaped Fermi surface (FS) enclosing two
neighboring Weyl nodes with the opposite chirality. The magnetotransport in
NbP is dominated by these highly anisotropic W2WSM pockets, in which Weyl
fermions are well protected from defect backscattering by real spin
conservation associated to the chiral nodes. However, with a minimal doping of
$\sim$1\% Cr, the mobility of NbP is degraded by more than two order of
magnitude, due to the invalid of helicity protection to magnetic impurities.
Helicity protected Weyl fermion transport is also manifested in chiral anomaly
induced negative magnetoresistance, controlled by the W1WSM states. In the
quantum regime below 10 K, the intervalley scattering time by impurities
becomes a large constant, producing the sharp and nearly identical conductivity
enhancement at low magnetic field.

We report extremely large positive magnetoresistance of 1.72 million percent
in single crystal TaSb$_{2}$ at moderate conditions of 1.5 K and 15 T. The
quadratic growth of magnetoresistance (MR $\propto\,B^{1.96}$) is not
saturating up to 15 T, a manifestation of nearly perfect compensation with
$<0.1\%$ mismatch between electron and hole pockets in this semimetal. The
compensation mechanism is confirmed by temperaturedependent MR, Hall and
thermoelectric coefficients of Nernst and Seebeck, revealing two pronounced
Fermi surface reconstruction processes without spontaneous symmetry breaking,
\textit{i.e.} Lifshitz transitions, at around 20 K and 60 K, respectively.
Using quantum oscillations of magnetoresistance and magnetic susceptibility,
supported by densityfunctional theory calculations, we determined that the
main hole Fermi surface of TaSb$_{2}$ forms a unique shoulder structure along
the $FL$ line. The flat band top of this shoulder pocket is just a few meV
above the Fermi level, leading to the observed topological phase transition at
20 K when the shoulder pocket disappears. Further increase in temperature
pushes the Fermi level to the band top of the main hole pocket, induced the
second Lifshitz transition at 60 K when hole pocket vanishes completely.

For the development of new 5G systems to operate in bands up to 100 GHz,
there is a need for accurate radio propagation models at these bands that
currently are not addressed by existing channel models developed for bands
below 6 GHz. This document presents a preliminary overview of 5G channel models
for bands up to 100 GHz. These have been derived based on extensive measurement
and ray tracing results across a multitude of frequencies from 6 GHz to 100
GHz, and this document describes an initial 3D channel model which includes: 1)
typical deployment scenarios for urban microcells (UMi) and urban macrocells
(UMa), and 2) a baseline model for incorporating path loss, shadow fading, line
of sight probability, penetration and blockage models for the typical
scenarios. Various processing methodologies such as clustering and antenna
decoupling algorithms are also presented.

In this paper, we report the intercone transport of Weyl fermions in NbAs
with external magnetic field in parallel to electric field, a quantum
phenomenon known as the AdlerBellJackiw anomaly. Surprisingly, the resulting
negative magnetoresistance (MR) in NbAs shows significant difference from NbP.
The observed lowfield positive MR dip, which is missing in NbP at low
temperatures, indicates that the spinorbital coupling (SOC) is significantly
stronger in NbAs than in the former. The results imply that the contribution of
arsenic to SOC in TaAs and NbAs is not negligible.

A profound assumption in peculiar velocity cosmology is $b_v=1$ at
sufficiently large scales, where $b_v$ is the volume weighted halo(galaxy)
velocity bias with respect to the matter velocity field. However, this
fundamental assumption has not been robustly verified in numerical simulations.
Furthermore, it is challenged by structure formation theory (BBKS, 1986, ApJ;
Desjacques and Sheth, 2010, PRD), which predicts the existence of velocity bias
(at least for protohalos) due to the fact that halos reside in special regions
(local density peaks). The major obstacle to measure the volume weighted
velocity from Nbody simulations is an unphysical sampling artifact. It is
entangled in the measured velocity statistics and becomes significant for
sparse populations. With recently improved understanding of the sampling
artifact (Zhang, Zheng and Jing, 2015, PRD; Zheng, Zhang and Jing, 2015, PRD),
for the first time we are able to {\it appropriately correct this sampling
artifact and then robustly measure the volume weighted halo velocity bias}. (1)
We verify $b_v=1$ within $2\%$ model uncertainty at $k\lesssim 0.1h/$Mpc and
$z=0$$2$ for halos of mass $\sim 10^{12}$$10^{13} h^{1} M_\odot$, and,
therefore, consolidates a foundation of the peculiar velocity cosmology. (2) We
also find statistically significant signs of $b_v\neq 1$ at $k\gtrsim
0.1h/$Mpc. Unfortunately, whether this is real or caused by residual sampling
artifact requires further investigation. Nevertheless, cosmology based on
$k\gtrsim 0.1h/$Mpc velocity data shall be careful this potential velocity
bias.

Measuring the volume weighted velocity power spectrum suffers from a severe
systematic error, due to imperfect sampling of the velocity field from
inhomogeneous distribution of dark matter particles/halos in simulations or
galaxies with velocity measurement. This "sampling artifact" depends on both
the mean particle number density $\bar{n}_P$ and the intrinsic large scale
structure (LSS) fluctuation in the particle distribution. (1) We report robust
detection of this sampling artifact in Nbody simulations. It causes $\sim 12$%
underestimation of the velocity power spectrum at $k=0.1$h/Mpc for samples with
$\bar{n}_P=6\times10^{3}$ (Mpc/h)$^{3}$. This systematic underestimation
increases with decreasing $\bar{n}_P$ and increasing $k$. Its dependence on the
intrinsic LSS fluctuations is also robustly detected. (2) All these findings
are expected by our theoretical modelling in paper I \cite{Zhang14}. In
particular, the leading order theoretical approximation agrees quantitatively
well with simulation result for $\bar{n}_P\gtrsim6\times
10^{4}$(Mpc/h)$^{3}$. Furthermore, we provide an ansatz to take high order
terms into account. It improves the model accuracy to $\lesssim1$% at
$k\lesssim0.1$h/Mpc over 3 orders of magnitude in $\bar{n}_P$ and over typical
LSS clustering from $z=0$ to $z=2$. (3) The sampling artifact is determined by
the deflection ${\bf D}$ field, which is straightforwardly available in both
simulations and data of galaxy velocity. Hence the sampling artifact in the
velocity power spectrum measurement can be selfcalibrated within our
framework. By applying such selfcalibration in simulations, it becomes
promising to determine the {\it real} large scale velocity bias of
$10^{13}M_\odot$ halos with $\sim 1$% accuracy, and that of lower mass halos by
better accuracy. ...[abridged]

We study the properties of the quantum states in the onedimensional system
with a shifted periodic potential in both the discrete model and the continuous
model. With open boundary conditions, the edge states appear in the energy gaps
which indicate nontrivial topological structures. The Chern numbers with
respect to the Bloch vector and the potential shift angle are computed. In the
limit of the continuous model, the Chern number of each band is exactly one. We
demonstrate the particle number pumped by the adiabatically shift of the
potential is directly related to the topological invariants.

Quantum and thermal fluctuations of electromagnetic fields, which give rise
to Planck's law of blackbody radiation, are also responsible for van der Waals
and Casimir forces, as well as nearfield radiative energy transfer between
objects. Electromagnetic waves transport energy, momentum, and entropy. For
classical thermal radiation, the dependence of the above mentioned quantities
on the temperature is wellknown mainly due to Planck's work. When nearfield
effects, namely the collective influence of diffraction, interference, and
tunneling of waves, become important, Planck's theory is no longer valid. Of
momentum, energy, and entropy transfer, the role of nearfield effects on
momentum transfer between two halfspaces separated by a vacuum gap (van der
Waals pressure in the vacuum gap) was first determined by Lifshitz, using
Rytov's theory of fluctuational electrodynamics in 1956. Subsequently,
Dzyaloshinskii, Lifshitz, and Pitaevskii, employing sophisticated methods from
quantum field theory, generalized Lifshitz' result for van der Waals pressure
in a vacuum layer to the case of van der Waals pressure in a dissipative layer
between two halfspaces. The influence of nearfield effects on radiative
transfer was appreciated only in the late 1960s and, subsequently, in the last
two decades because of the enhancement in radiative transfer due to
electromagnetic surface waves. The role played by nearfield effects on entropy
transfer has not been investigated so far, at least when the temperature
distribution is nonuniform.

We use fluctuational electrodynamics to determine emissivity and van der
Waals contribution to surface energy for various spherical shapes, such as a
sphere, a bubble, a spherical shell and a coated sphere, in a homogeneous and
isotropic medium. Nearfield radiative transfer and momentum transfer between
flat plates and curved surfaces have been studied for the past decades, however
the investigation of radiative heat transfer and van der Waals stress due to
fluctuations of electromagnetic fields for a single object is missing from
literature. The dyadic Green's function formalism of radiative energy and
fluctuationinduced van der Waals stress for different spherical configurations
have been developed. We show (1) emission spectra of micro and nanosized
spheres display several emissivity sharp peaks as the size of object reduces,
and (2) surface energy becomes size dependent due to van der Waals phenomena
when size of object is reduced to a nanoscopic length scale.

Printed electronics is quickly emerging with tremendous value in a wide
variety of latest electrical engineering areas. However, restricted to the
rather limited conductive inks and printing strategies, the currently existing
electronics manufacturing tools are generally complicated, expensive, time,
material, water and energy consuming, and thus mainly restricted to the
industrial use. Here from an alternative way, the present article demonstrated
for the first time an entirely automatic printing system for personal
electronics manufacturing through introducing a composite liquid metal ink
delivery and printing mechanism to overcome the large surface tension facing
the solution, and integrating it with the notebook computer controlling
algorithm. With the developed printer, any desired electronically conductive
patterns spanning from single wires to various complex structures like
integrated circuit (IC), antenna, PCB, RFID, electronic greeting cards,
decoration arts, classical buildings (White House, Great Wall etc.) or more
doityourself (DIY) circuits were demonstrated to be printed out in a moment
with high precision. And the total cost for the whole system has reached
personal affordability, which is hard to offer by so far the state of the art
technologies. Some fundamental fluid dynamics mechanisms related to the
proposed tapping mode enabled reliable printing and adhesion of the liquid
metal electronics on the flexible substrate was systematically disclosed
through theoretical interpretation and experimental measurements. This clearly
beyondthelab technology and pervasively available liquid metal printer opens
the way for large scale home level electronics making in the coming time.

A rollerball pen enabled direct writing electronics via room temperature
liquid metal ink was proposed. With the rolling to print mechanism, the
metallic inks were smoothly written on flexible polymer substrate to form
conductive tracks and electronic devices. The contact angle analyzer and
scanning electron microscope were implemented to probe the inner property of
the obtained electronics. An ever high writing resolution with line width and
thickness as 200{\mu}m and 80{\mu}m, respectively was realized. Further, with
the administration of external writing pressure, GaIn24.5 droplets embody
increasing wettability on polymer which demonstrates the pervasive adaptability
of the rollerball pen electronics.