
Simulations of single and multispecies compressible flows with shock waves
and discontinuities are conducted using a weighted compact nonlinear scheme
(WCNS) with a newly developed sixth order localized dissipative interpolation.
In smooth regions, the scheme applies the central nonlinear interpolation with
minimum dissipation to well resolve fluctuating flow features while in regions
containing discontinuities and high wavenumber features, the scheme suppresses
spurious numerical oscillations by hybridizing the central interpolation with
the more dissipative upwindbiased nonlinear interpolation. In capturing
material interfaces between species of different densities, a
quasiconservative five equation model that can conserve mass of each species
is used to prevent pressure oscillations across the interfaces. Compared to
upwindbiased interpolations with classical nonlinear weights and improved
weights, and the interpolation with adaptive centralupwind weights for
scaleseparation, it is shown that WCNS with the proposed localized dissipative
interpolation has better performance to simultaneously capture discontinuities
and resolve smooth features.

We present a new magnetohydrodynamic (MHD) simulation code with the aim of
providing accurate numerical solutions to astrophysical phenomena where
discontinuities, shock waves, and turbulence are inherently important. The code
implements the HLLD approximate Riemann solver, the
fifthordermonotonicitypreserving interpolation (MP5) scheme, and the
hyperbolic divergence cleaning method for a magnetic field. This choice of
schemes significantly improved numerical accuracy and stability, and saved
computational costs in multidimensional problems. Numerical tests of one and
twodimensional problems showed the advantages of using the highorder scheme
by comparing with results from a standard secondorder TVD MUSCL scheme. The
present code enabled us to explore longterm evolution of a threedimensional
accretion disk around a black hole, in which compressible MHD turbulence caused
continuous mass accretion via nonlinear growth of the magnetorotational
instability (MRI). Numerical tests with various computational cell sizes
exhibited a convergent picture of the early nonlinear growth of the MRI in a
global model, and indicated that the MP5 scheme has more than twice the
resolution of the MUSCL scheme in practical applications.

Density functional theory (DFT) calculation has had huge success as a tool
capable of predicting important physical and chemical properties of condensed
matter systems. We calculate the electric dipole moment of a molecule by using
the differential electron density with respect to the superimposed electron
density of the free atoms, exploiting the cancellation of DFT errors. Our
results on a range of molecules show an excellent agreement with experiments.

The resemblance between the methods used in quantummany body physics and in
machine learning has drawn considerable attention. In particular, tensor
networks (TNs) and deep learning architectures bear striking similarities to
the extent that TNs can be used for machine learning. Previous results used
onedimensional TNs in image recognition, showing limited scalability and
flexibilities. In this work, we train twodimensional hierarchical TNs to solve
image recognition problems, using a training algorithm derived from the
multiscale entanglement renormalization ansatz. This approach introduces
mathematical connections among quantum manybody physics, quantum information
theory, and machine learning. While keeping the TN unitary in the training
phase, TN states are defined, which encode classes of images into quantum
manybody states. We study the quantum features of the TN states, including
quantum entanglement and fidelity. We find these quantities could be properties
that characterize the image classes, as well as the machine learning tasks.

We develop a generalization of the KohnSham density functional theory
(KSDFT) + Hubbard $U$ (DFT+$U$) method to the excitedstate regime. This has
the form of Hubbard $U$ corrected linearresponse timedependent DFT, or
`TDDFT+$U$'. Combined with calculated linearresponse Hubbard $U$ parameters,
it may provide a computationally light, firstprinciples method for the
simulation of tightlybound excitons on transitionmetal ions. Our presented
implementation combines linearscaling DFT+$U$ and linearscaling TDDFT, but
the approach is broadly applicable. In detailed benchmark tests on two
Nicentred diamagnetic coordination complexes with variable $U$ values, it is
shown that the Hubbard $U$ correction to an approximate adiabatic semilocal
exchangecorrelation interaction kernel lowers the excitation energies of
transitions exclusively within the targeted localized subspace, by increasing
the exciton binding of the corresponding electronhole pairs. This partially
counteracts the Hubbard $U$ correction to the exchangecorrelation potential in
KSDFT, which increases excitation energies into, out of, and within the
targeted localised subspace by modifying the underlying KSDFT eigenspectrum.
This compensating effect is most pronounced for optically dark transitions
between localized orbitals of the same angular momentum, for which experimental
observation may be challenging and theoretical approaches are at their most
necessary. Overall, our results point to shortcomings in the contemporary
DFT+$U$ corrective potential, either in its functional form, or when applied to
transitionmetal orbitals but not to ligand ones, or both.

In this work, an efficient numerical scheme is presented for seismic blind
deconvolution in a multichannel scenario. The proposed method iterate with wo
steps: first, wavelet estimation across all channels and second, refinement of
the reflectivity estimate simultaneously in all channels using sparse
deconvolution. The reflectivity update step is formulated as a basis pursuit
denoising problem and a sparse solution is obtained with the spectral
projectedgradient algorithm  faithfulness to the recorded traces is
constrained by the measured noise level. Wavelet reestimation has a closed
form solution when performed in the frequency domain by finding the minimum
energy wavelet common to all channels. Nothing is assumed known about the
wavelet apart from its time duration. In tests with both synthetic and real
data, the method yields sparse reflectivity series and stable wavelet estimates
results compared to existing methods with significantly less computational
effort.

The anelastic and pseudoincompressible equations are two wellknown
soundproof approximations of compressible flows useful for both theoretical and
numerical analysis in meteorology, atmospheric science, and ocean studies. In
this paper, we derive and test structurepreserving numerical schemes for these
two systems. The derivations are based on a discrete version of the
EulerPoincar\'e variational method. This approach relies on a finite
dimensional approximation of the (Lie) group of diffeomorphisms that preserve
weightedvolume forms. These weights describe the background stratification of
the fluid and correspond to the weighed velocity fields for anelastic and
pseudoincompressible approximations. In particular, we identify to these
discrete Lie group configurations the associated Lie algebras such that
elements of the latter correspond to weighted velocity fields that satisfy the
divergencefree conditions for both systems. Defining discrete Lagrangians in
terms of these Lie algebras, the discrete equations follow by means of
variational principles. Descending from variational principles, the schemes
exhibit further a discrete version of Kelvin circulation theorem, are
applicable to irregular meshes, and show excellent long term energy behavior.
We illustrate the properties of the schemes by performing preliminary test
cases.

Reversible, diffusionless, firstorder solidsolid phase transitions
accompanied by caloric effects are critical for applications in the solidstate
cooling and heatpumping devices. Accelerated discovery of caloric materials
requires reliable but faster estimators for predictions and highthroughput
screening of systemspecific dominant caloric contributions. We assess
reliability of the computational methods that provide thermodynamic properties
in relevant solid phases at or near a phase transition. We test the methods
using the wellstudied B2 FeRh alloy as a "fruit fly" in such a materials
genome discovery, as it exhibits a metamagnetic transition which generates
multicaloric (magneto, elasto, and barocaloric) responses. For lattice
entropy contributions, we find that the commonlyused linearresponse and
smalldisplacement phonon methods are invalid near instabilities that arise
from the anharmonicity of atomic potentials, and we offer a more reliable and
precise method for calculating lattice entropy at a fixed temperature. Then, we
apply a set of reliable methods and estimators to the metamagnetic transition
in FeRh (predicted $346 \pm 12$ K, observed $353 \pm 1$ K) and calculate the
associated caloric properties, such as isothermal entropy and isentropic
temperature changes.

We provide a description of a diagrammatic Monte Carlo algorithm for the
resonant Fermi gas in the normal phase. Details are given on diagrammatic
framework, Monte Carlo moves, and incorporation of ultraviolet asymptotics.
Apart from the selfconsistent bold scheme, we also describe a
nonselfconsistent scheme, for which the ultraviolet treatment is more
involved.

A new MHD solver, based on the Nektar++ spectral/hp element framework, is
presented in this paper. The velocity and electric potential quasistatic MHD
model is used. The Hartmann flow in plane channel and its stability, the
Hartmann flow in rectangular duct, and the stability of Hunt's flow are
explored as examples. Exponential convergence is achieved and the resulting
numerical values were found to have an accuracy up to $10^{12}$ for the state
flows compared to an exact solution, and $10^{5}$ for the stability
eigenvalues compared to independent numerical results.

Understanding the mechanism of the heterojunction is an important step
towards controllable and tunable interfaces for photocatalytic and photovoltaic
based devices. To this aim, we propose a thorough study of a double
heterostructure system consisting of two semiconductors with large band gap,
namely, wurtzite ZnO and anatase TiO2. We demonstrate via firstprinciple
calculations two stable configurations of ZnO/TiO2 interfaces. Our structural
analysis provides a key information on the nature of the complex interface and
lattice distortions occurring when combining these materials. The study of the
electronic properties of the sandwich nanostructure TiO2/ZnO/TiO2 reveals that
conduction band arises mainly from Ti3d orbitals, while valence band is
maintained by O2p of ZnO, and that the trapped states within the gap region
frequent in single heterostructure are substantially reduced in the double
interface system. Moreover, our work explains the origin of certain optical
transitions observed in the experimental studies. Unexpectedly, as a
consequence of different bond distortions, the results on the band alignments
show electron accumulation in the left shell of TiO2 rather than the right one.
Such behavior provides more choice for the sensitization and functionalization
of TiO2 surfaces.

The effect of increased electrondensity (from adsorbed Li atoms) in
polyacenes and in nanoribbons with zigzag edge is discussed in terms of
resonance theoretical considerations and in terms edgelocalized frontier
molecular orbitals. The argumentation from simple pictures is finally using the
density functional theory (DFT) for anthracene, polyacene polymer and graphene
strips. Some discussion is made for zigzag edge graphene.

Based on the particleincell (PIC) plasma simulation method, the
speedlimited PIC (SLPIC) method delivers faster kinetic plasma simulation in
cases where the particle distributions evolve slowly compared with the maximum
stable PIC timestep. SLPIC thus offers more feasible, fully kinetic simulation
in regimes that historically have required fluid approaches, such as
magnetohydrodynamic (MHD), twofluid, or Boltzmann electron treatments. In
particular, SLPIC allows an explicit time advance with steps much larger than
the inverse plasma frequency, avoiding the instability explicit PIC faces with
large timesteps. SLPIC avoids this instability by slowing down fast particles
(e.g., electrons) in a way that is rigorously underpinned by an approximate
Vlasov equation; unlike MHD, twofluid, and Boltzmann electron approaches,
SLPIC does not fundamentally neglect any firstprinciples plasma physics,
although the choices of grid cell size, timestep, and number of macroparticles
per cell naturally limit the physical phenomena that can be accurately
represented. SLPIC can be implemented with minor modifications of a standard
PIC code and does not require an implicit time advance. It enables large
timesteps in firstprinciples kinetic plasma simulation of appropriately slow
phenomena, and it can handle many of the same complications as PIC, such as
boundary conditions and collisions. In an argon plasma sheath test problem, a
SLPIC simulation achieved a speedup of a factor of 160 over the corresponding
PIC simulation, without loss of accuracy.

A unified fluidstructure interaction (FSI) formulation is presented for
solid, liquid and mixed membranes. Nonlinear finite elements (FE) and the
generalizedalpha scheme are used for the spatial and temporal discretization.
The membrane discretization is based on curvilinear surface elements that can
describe large deformations and rotations, and also provide a straightforward
description for contact. The fluid is described by the incompressible
NavierStokes equations, and its discretization is based on stabilized
PetrovGalerkin FE. The coupling between fluid and structure uses a conforming
sharp interface discretization, and the resulting nonlinear FE equations are
solved monolithically within the NewtonRaphson scheme. An arbitrary
LagrangianEulerian formulation is used for the fluid in order to account for
the mesh motion around the structure. The formulation is very general and
admits diverse applications that include contact at free surfaces. This is
demonstrated by two analytical and three numerical examples exhibiting strong
coupling between fluid and structure. The examples include balloon inflation,
droplet rolling and flapping flags. They span a Reynoldsnumber range from
0.001 to 2000. One of the examples considers the extension to rotationfree
shells using isogeometric FE.

We test numerically the recently proposed linear relationship between the
scaleinvariant period $T_{\rm s.i.} = T E^{3/2}$, and the topology of an
orbit, on several hundred planar Newtonian periodic threebody orbits. Here $T$
is the period of an orbit, $E$ is its energy, so that $T_{\rm s.i.}$ is the
scaleinvariant (s.i.) period, or, equivalently, the period at unit energy $E
= 1$. All of these orbits have vanishing angular momentum and pass through a
linear, equidistant configuration at least once. Such orbits are classified in
ten algebraically welldefined sequences. Orbits in each sequence follow an
approximate linear dependence of $T_{\rm s.i.}$, albeit with slightly different
slopes and intercepts. The orbit with the shortest period in its sequence is
called the "progenitor": six distinct orbits are the progenitors of these ten
sequences. We have studied linear stability of these orbits, with the result
that 21 orbits are linearly stable, which includes all of the progenitors. This
is consistent with the BirkhoffLewis theorem, which implies existence of
infinitely many periodic orbits for each stable progenitor, and in this way
explains the existence and ensures infinite extension of each sequence.

We study three computer algebra systems, namely SageMath (with SageManifolds
package), Maxima (with ctensor package) and Python language (with GraviPy
module), which allow tensor manipulation for general relativity calculations
along with general algebraic calculations. We present a benchmark of these
systems using simple examples. After the general analysis, we focus on the
SageMath and SageManifolds system to derive, analyze and visualize the
solutions of the massless KleinGordon equation and geodesic motion with
HamiltonJacobi formalism. We compare our numerical result of the KleinGordon
equation with the asymptotic form of the analytical solution to see that they
agree.

The scattering of electromagnetic pulses is described using a nonsingular
boundary integral method to solve directly for the field components in the
frequency domain, and Fourier transform is then used to obtain the complete
spacetime behavior. This approach is stable for wavelengths both small and
large relative to characteristic length scales. Amplitudes and phases of field
values can be obtained accurately on or near material boundaries. Local field
enhancement effects due to multiple scattering of interest to applications in
microphotonics are demonstrated.

We present a boundary integral formulation of electromagnetic scattering by
homogeneous bodies that are characterized by linear constitutive equations in
the frequency domain. By working with the Cartesian components of the electric,
E and magnetic, H fields and with the scalar functions (r*E) and (r*H), the
problem is cast as solving a set of scalar Helmholtz equations for the field
components that are coupled by the usual electromagnetic boundary conditions at
material boundaries. This facilitates a direct solution for E and H rather than
working with surface currents as intermediate quantities in existing methods.
Consequently, our formulation is free of the wellknown numerical instability
that occurs in the zero frequency or long wavelength limit in traditional
surface integral solutions of Maxwell's equations and our numerical results
converge uniformly to the static results in the long wavelength limit.
Furthermore, we use a formulation of the scalar Helmholtz equation that is
expressed as classically convergent integrals and does not require the
evaluation of principal value integrals or any knowledge of the solid angle.
Therefore, standard quadrature and higher order surface elements can readily be
used to improve numerical precision. In addition, near and far field values can
be calculated with equal precision and multiscale problems in which the
scatterers possess characteristic length scales that are both large and small
relative to the wavelength can be easily accommodated. From this we obtain
results for the scattering and transmission of electromagnetic waves at
dielectric boundaries that are valid for any ratio of the local surface
curvature to the wave number. This is a generalization of the familiar Fresnel
formula and Snell's law, valid at planar dielectric boundaries, for the
scattering and transmission of electromagnetic waves at surfaces of arbitrary
curvature.

A boundary integral formulation of electromagnetics that involves only the
components of $\boldsymbol{E}$ and $\boldsymbol{H}$ is derived without the use
of surface currents that appear in the classical PMCHWT formulation. The
kernels of the boundary integral equations for $\boldsymbol{E}$ and
$\boldsymbol{H}$ are nonsingular so that all field quantities at the surface
can be determined to high precision and also geometries with closely spaced
surfaces present no numerical difficulties. Quadratic elements can readily be
used to represent the surfaces so that the surface integrals can be calculated
to higher numerical precision than using planar elements for the same numbers
of degrees of freedom.

Of the two main objectives we pursue in this paper, the first one consists in
the studying operators of the form
$(\partial_ti\triangle_{\Gamma})^{\alpha},\,\,\alpha=1/2,1/2,1,\ldots,$
where $\triangle_{\Gamma}$ is the LaplaceBeltrami operator. These operators
arise in the context of nonreflecting boundary conditions in the
pseudodifferential approach for the general Schr\"odinger equation. The
definition of such operators is discussed in various settings and a formulation
in terms of fractional operators is provided. The second objective consists in
deriving corner conditions for a rectangular domain in order to make such
domains amenable to the pseudodifferential approach. Stability and uniqueness
of the solution is investigated for each of these novel boundary conditions.

Quantum++ is a modern generalpurpose multithreaded quantum computing
library written in C++11 and composed solely of header files. The library is
not restricted to qubit systems or specific quantum information processing
tasks, being capable of simulating arbitrary quantum processes. The main design
factors taken in consideration were the ease of use, portability, and
performance. The library's simulation capabilities are only restricted by the
amount of available physical memory. On a typical machine (Intel i5 8Gb RAM)
Quantum++ can successfully simulate the evolution of 25 qubits in a pure state
or of 12 qubits in a mixed state reasonably fast. The library also includes
support for classical reversible logic, being able to simulate classical
reversible operations on billions of bits. This latter feature may be useful in
testing quantum circuits composed solely of Toffoli gates, such as certain
arithmetic circuits.

Despite its numerical challenges, finite element method is used to compute
viscous fluid flow. A consensus on the cause of numerical problems has been
reached; however, general algorithmsallowing a robust and accurate
simulation for any processare still missing. Either a very high
computational cost is necessary for a direct numerical solution (DNS) or some
limiting procedure is used by adding artificial dissipation to the system.
These stabilization methods are useful; however, they are often applied
relative to the element size such that a local monotonous convergence is
challenging to acquire. We need a computational strategy for solving viscous
fluid flow using solely the balance equations. In this work, we present a
general procedure solving fluid mechanics problems without use of any
stabilization or splitting schemes. Hence, its generalization to multiphysics
applications is straightforward. We discuss emerging numerical problems and
present the methodology rigorously. Implementation is achieved by using
opensource packages and the accuracy as well as the robustness is demonstrated
by comparing results to the closedform solutions and also by solving
wellknown benchmarking problems.

Homogeneous droplet nucleation has been studied for almost a century but has
not yet been fully understood. In this work, we used the density gradient
theory (DGT) and considered the influence of capillary waves (CW) on the
predicted sizedependent surface tensions and nucleation rates for selected
$n$alkanes. The DGT model was completed by an equation of state (EoS) based on
the perturbedchain statistical associating fluid theory (PCSAFT) and compared
to the classical nucleation theory and the PengRobinson EoS. It was found
that the critical clusters are practically free of CW because they are so small
that even the smallest CW wavelengths do not fit into their finite dimensions.
The CW contribute to the entropy of the system and thus decrease the surface
tension. A correction for the effect of CW on the surface tension is presented.
The effect of the different EoSs is relatively small because by a fortuitous
coincidence their predictions are similar in the relevant range of critical
cluster sizes. The difference of the DGT predictions to the classical
nucleation theory computations is important but not decisive. Of the effects
investigated, the most pronounced is the suppression of the CW which causes a
sizable decrease of the predicted nucleation rates. The major difference
between experimental nucleation rate data and theoretical predictions remains
in the temperature dependence. For normal alkanes, this discrepancy is much
stronger than observed, e.g., for water. The theoretical corrections developed
here have a minor influence on the temperature dependency. We provide empirical
equations correcting the predicted nucleation rates to values comparable with
experiments.

We apply recent advances in machine learning and computer vision to a central
problem in materials informatics: The statistical representation of
microstructural images. We use activations in a pretrained convolutional
neural network to provide a highdimensional characterization of a set of
synthetic microstructural images. Next, we use manifold learning to obtain a
lowdimensional embedding of this statistical characterization. We show that
the lowdimensional embedding extracts the parameters used to generate the
images. According to a variety of metrics, the convolutional neural network
method yields dramatically better embeddings than the analogous method derived
from twopoint correlations alone.

Tensors are a natural way to express correlations among many physical
variables, but storing tensors in a computer naively requires memory which
scales exponentially in the rank of the tensor. This is not optimal, as the
required memory is actually set not by the rank but by the mutual information
amongst the variables in question. Representations such as the tensor tree
perform nearoptimally when the tree decomposition is chosen to reflect the
correlation structure in question, but making such a choice is nontrivial and
good heuristics remain highly contextspecific. In this work I present two new
algorithms for choosing efficient tree decompositions, independent of the
physical context of the tensor. The first is a bruteforce algorithm which
produces optimal decompositions up to truncation error but is generally
impractical for highrank tensors, as the number of possible choices grows
exponentially in rank. The second is a greedy algorithm, and while it is not
optimal it performs extremely well in numerical experiments while having
runtime which makes it practical even for tensors of very high rank.