• ### High-Order Localized Dissipation Weighted Compact Nonlinear Scheme for Shock- and Interface-Capturing in Compressible Flows(1701.08905)

Simulations of single- and multi-species 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 upwind-biased nonlinear interpolation. In capturing material interfaces between species of different densities, a quasi-conservative five equation model that can conserve mass of each species is used to prevent pressure oscillations across the interfaces. Compared to upwind-biased interpolations with classical nonlinear weights and improved weights, and the interpolation with adaptive central-upwind weights for scale-separation, it is shown that WCNS with the proposed localized dissipative interpolation has better performance to simultaneously capture discontinuities and resolve smooth features.
• ### Magnetohydrodynamic Simulation Code CANS+: Assessments and Applications(1611.01775)

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 fifth-order-monotonicity-preserving 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 two-dimensional problems showed the advantages of using the high-order scheme by comparing with results from a standard second-order TVD MUSCL scheme. The present code enabled us to explore long-term evolution of a three-dimensional accretion disk around a black hole, in which compressible MHD turbulence caused continuous mass accretion via nonlinear growth of the magneto-rotational 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.
• ### Determination of the Electric Dipole Moment of a Molecule from Density Functional Theory Calculations(1708.03834)

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.
• ### Machine Learning by Unitary Tensor Network of Hierarchical Tree Structure(1710.04833)

The resemblance between the methods used in quantum-many 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 one-dimensional TNs in image recognition, showing limited scalability and flexibilities. In this work, we train two-dimensional hierarchical TNs to solve image recognition problems, using a training algorithm derived from the multi-scale entanglement renormalization ansatz. This approach introduces mathematical connections among quantum many-body 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 many-body 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.
• ### TDDFT+$U$: Hubbard corrected approximate density-functional theory in the excited-state regime(1703.05056)

We develop a generalization of the Kohn-Sham density functional theory (KS-DFT) + Hubbard $U$ (DFT+$U$) method to the excited-state regime. This has the form of Hubbard $U$ corrected linear-response time-dependent DFT, or `TDDFT+$U$'. Combined with calculated linear-response Hubbard $U$ parameters, it may provide a computationally light, first-principles method for the simulation of tightly-bound excitons on transition-metal ions. Our presented implementation combines linear-scaling DFT+$U$ and linear-scaling TDDFT, but the approach is broadly applicable. In detailed benchmark tests on two Ni-centred diamagnetic coordination complexes with variable $U$ values, it is shown that the Hubbard $U$ correction to an approximate adiabatic semi-local exchange-correlation interaction kernel lowers the excitation energies of transitions exclusively within the targeted localized subspace, by increasing the exciton binding of the corresponding electron-hole pairs. This partially counteracts the Hubbard $U$ correction to the exchange-correlation potential in KS-DFT, which increases excitation energies into, out of, and within the targeted localised subspace by modifying the underlying KS-DFT 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 transition-metal orbitals but not to ligand ones, or both.
• ### Sparse blind deconvolution of seismic data via spectral projected-gradient(1611.03754)

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 projected-gradient algorithm - faithfulness to the recorded traces is constrained by the measured noise level. Wavelet re-estimation 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.
• ### Variational integrators for anelastic and pseudo-incompressible flows(1701.06448)

Feb. 2, 2019 math.NA, physics.comp-ph
The anelastic and pseudo-incompressible equations are two well-known 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 structure-preserving numerical schemes for these two systems. The derivations are based on a discrete version of the Euler-Poincar\'e variational method. This approach relies on a finite dimensional approximation of the (Lie) group of diffeomorphisms that preserve weighted-volume forms. These weights describe the background stratification of the fluid and correspond to the weighed velocity fields for anelastic and pseudo-incompressible 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 divergence-free 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.
• ### Reliable thermodynamic estimators for screening multicaloric materials(1702.03042)

Reversible, diffusionless, first-order solid-solid phase transitions accompanied by caloric effects are critical for applications in the solid-state cooling and heat-pumping devices. Accelerated discovery of caloric materials requires reliable but faster estimators for predictions and high-throughput screening of system-specific 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 well-studied 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 baro-caloric) responses. For lattice entropy contributions, we find that the commonly-used linear-response and small-displacement 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.
• ### Diagrammatic Monte Carlo algorithm for the resonant Fermi gas(1305.3901)

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 self-consistent bold scheme, we also describe a non-self-consistent scheme, for which the ultraviolet treatment is more involved.
• ### A spectral/hp element MHD solver(1707.08957)

Jan. 21, 2019 physics.comp-ph
A new MHD solver, based on the Nektar++ spectral/hp element framework, is presented in this paper. The velocity and electric potential quasi-static 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.
• ### Mechanism of the double heterostructure TiO2/ZnO/TiO2 for photocatalytic and photovoltaic applications: A theoretical study(1702.02545)

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 first-principle 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.
• ### Lithium Adsorbtion on Polyacenes $\&$ Zig-zag-edge Graphene Strips(1603.09308)

The effect of increased electron-density (from adsorbed Li atoms) in polyacenes and in nano-ribbons with zig-zag edge is discussed in terms of resonance theoretical considerations and in terms edge-localized 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 zig-zag edge graphene.
• ### Speeding Up Simulations By Slowing Down Particles: Speed-Limited Particle-In-Cell Simulation(1511.08225)

Based on the particle-in-cell (PIC) plasma simulation method, the speed-limited 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), two-fluid, 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, two-fluid, and Boltzmann electron approaches, SLPIC does not fundamentally neglect any first-principles 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 first-principles 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 speed-up of a factor of 160 over the corresponding PIC simulation, without loss of accuracy.
• ### A monolithic fluid-structure interaction formulation for solid and liquid membranes including free-surface contact(1710.02128)

Dec. 27, 2018 physics.comp-ph, cs.CE
A unified fluid-structure interaction (FSI) formulation is presented for solid, liquid and mixed membranes. Nonlinear finite elements (FE) and the generalized-alpha 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 Navier-Stokes equations, and its discretization is based on stabilized Petrov-Galerkin FE. The coupling between fluid and structure uses a conforming sharp interface discretization, and the resulting non-linear FE equations are solved monolithically within the Newton-Raphson scheme. An arbitrary Lagrangian-Eulerian 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 Reynolds-number range from 0.001 to 2000. One of the examples considers the extension to rotation-free shells using isogeometric FE.
• ### Linear stability of periodic three-body orbits with zero angular momentum and topological dependence of Kepler's third law: a numerical test(1705.03728)

We test numerically the recently proposed linear relationship between the scale-invariant period $T_{\rm s.i.} = T |E|^{3/2}$, and the topology of an orbit, on several hundred planar Newtonian periodic three-body orbits. Here $T$ is the period of an orbit, $E$ is its energy, so that $T_{\rm s.i.}$ is the scale-invariant (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 well-defined 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 Birkhoff-Lewis 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.
• ### Symbolic and Numerical Analysis in General Relativity with Open Source Computer Algebra Systems(1703.09738)

Dec. 18, 2018 hep-th, gr-qc, physics.comp-ph
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 Klein-Gordon equation and geodesic motion with Hamilton-Jacobi formalism. We compare our numerical result of the Klein-Gordon equation with the asymptotic form of the analytical solution to see that they agree.
• ### Field-only integral equation method for time domain scattering of electromagnetic pulses(1708.07934)

The scattering of electromagnetic pulses is described using a non-singular boundary integral method to solve directly for the field components in the frequency domain, and Fourier transform is then used to obtain the complete space-time 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.
• ### A Robust Multi-Scale Field-Only Formulation of Electromagnetic Scattering(1611.01978)

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 well-known 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.
• ### Non-singular field-only surface integral equations for electromagnetic scattering(1608.02058)

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 non-singular 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.
• ### On the nonreflecting boundary operators for the general two dimensional Schr\"odinger equation(1502.04519)

Dec. 16, 2018 math.AP, physics.comp-ph
Of the two main objectives we pursue in this paper, the first one consists in the studying operators of the form $(\partial_t-i\triangle_{\Gamma})^{\alpha},\,\,\alpha=1/2,-1/2,-1,\ldots,$ where $\triangle_{\Gamma}$ is the Laplace-Beltrami operator. These operators arise in the context of nonreflecting boundary conditions in the pseudo-differential 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 pseudo-differential approach. Stability and uniqueness of the solution is investigated for each of these novel boundary conditions.
• ### Quantum++: A modern C++ quantum computing library(1412.4704)

Dec. 11, 2018 quant-ph, physics.comp-ph
Quantum++ is a modern general-purpose multi-threaded 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.
• ### An accurate finite element method for the numerical solution of isothermal and incompressible flow of viscous fluid(1709.00913)

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 algorithms---allowing a robust and accurate simulation for any process---are 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 open-source packages and the accuracy as well as the robustness is demonstrated by comparing results to the closed-form solutions and also by solving well-known benchmarking problems.
• ### Predictions of homogeneous nucleation rates for $n$-alkanes accounting for the diffuse phase interface and capillary waves(1705.07588)

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 size-dependent surface tensions and nucleation rates for selected $n$-alkanes. The DGT model was completed by an equation of state (EoS) based on the perturbed-chain statistical associating fluid theory (PC-SAFT) and compared to the classical nucleation theory and the Peng--Robinson 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.
• ### Inferring low-dimensional microstructure representations using convolutional neural networks(1611.02764)

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 pre-trained convolutional neural network to provide a high-dimensional characterization of a set of synthetic microstructural images. Next, we use manifold learning to obtain a low-dimensional embedding of this statistical characterization. We show that the low-dimensional 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 two-point correlations alone.
• ### Efficient Decomposition of High-Rank Tensors(1708.07471)

Nov. 29, 2018 cs.NA, physics.comp-ph
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 near-optimally when the tree decomposition is chosen to reflect the correlation structure in question, but making such a choice is non-trivial and good heuristics remain highly context-specific. 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 brute-force algorithm which produces optimal decompositions up to truncation error but is generally impractical for high-rank 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.