25,916 research outputs found

    Quantum Mechanics Lecture Notes. Selected Chapters

    Full text link
    These are extended lecture notes of the quantum mechanics course which I am teaching in the Weizmann Institute of Science graduate physics program. They cover the topics listed below. The first four chapter are posted here. Their content is detailed on the next page. The other chapters are planned to be added in the coming months. 1. Motion in External Electromagnetic Field. Gauge Fields in Quantum Mechanics. 2. Quantum Mechanics of Electromagnetic Field 3. Photon-Matter Interactions 4. Quantization of the Schr\"odinger Field (The Second Quantization) 5. Open Systems. Density Matrix 6. Adiabatic Theory. The Berry Phase. The Born-Oppenheimer Approximation 7. Mean Field Approaches for Many Body Systems -- Fermions and Boson

    Full trajectory optimizing operator inference for reduced-order modeling using differentiable programming

    Full text link
    Accurate and inexpensive Reduced Order Models (ROMs) for forecasting turbulent flows can facilitate rapid design iterations and thus prove critical for predictive control in engineering problems. Galerkin projection based Reduced Order Models (GP-ROMs), derived by projecting the Navier-Stokes equations on a truncated Proper Orthogonal Decomposition (POD) basis, are popular because of their low computational costs and theoretical foundations. However, the accuracy of traditional GP-ROMs degrades over long time prediction horizons. To address this issue, we extend the recently proposed Neural Galerkin Projection (NeuralGP) data driven framework to compressibility-dominated transonic flow, considering a prototypical problem of a buffeting NACA0012 airfoil governed by the full Navier-Stokes equations. The algorithm maintains the form of the ROM-ODE obtained from the Galerkin projection; however coefficients are learned directly from the data using gradient descent facilitated by differentiable programming. This blends the strengths of the physics driven GP-ROM and purely data driven neural network-based techniques, resulting in a computationally cheaper model that is easier to interpret. We show that the NeuralGP method minimizes a more rigorous full trajectory error norm compared to a linearized error definition optimized by the calibration procedure. We also find that while both procedures stabilize the ROM by displacing the eigenvalues of the linear dynamics matrix of the ROM-ODE to the complex left half-plane, the NeuralGP algorithm adds more dissipation to the trailing POD modes resulting in its better long-term performance. The results presented highlight the superior accuracy of the NeuralGP technique compared to the traditional calibrated GP-ROM method

    Chiral active fluids: Odd viscosity, active turbulence, and directed flows of hydrodynamic microrotors

    Get PDF
    While the number of publications on rotating active matter has rapidly increased in recent years, studies on purely hydrodynamically interacting rotors on the microscale are still rare, especially from the perspective of particle based hydrodynamic simulations. The work presented here targets to fill this gap. By means of high-performance computer simulations, performed in a highly parallelised fashion on graphics processing units, the dynamics of ensembles of up to 70,000 rotating colloids immersed in an explicit mesoscopic solvent consisting out of up to 30 million fluid particles, are investigated. Some of the results presented in this thesis have been worked out in collaboration with experimentalists, such that the theoretical considerations developed in this thesis are supported by experiments, and vice versa. The studied system, modelled in order to resemble the essential physics of the experimentally realisable system, consists out of rotating magnetic colloidal particles, i.e., (micro-)rotors, rotating in sync to an externally applied magnetic field, where the rotors solely interact via hydrodynamic and steric interactions. Overall, the agreement between simulations and experiments is very good, proving that hydrodynamic interactions play a key role in this and related systems. While already an isolated rotating colloid is driven out of equilibrium, only collections of two or more rotors have experimentally shown to be able to convert the rotational energy input into translational dynamics in an orbital rotating fashion. The rotating colloids inject circular flows into the fluid, such that detailed balance is broken, and it is not a priori known whether equilibrium properties of colloids can be extended to isolated rotating colloids. A joint theoretical and experimental analysis of isolated, pairs, and small groups of hydrodynamically interacting rotors is given in chapter 2. While the translational dynamics of isolated rotors effectively resemble the dynamics of non-rotating colloids, the orbital rotation of pairs of rotors can be described with leading order hydrodynamics and a two-dimensional analogy of Faxén’s law is derived. In chapter 3, a homogeneously distributed ensemble of rotors (bulk) as a realisation of a chiral active fluid is studied and it is explicitly shown computationally and experimentally that it carries odd viscosity. The mutual orbital translation of rotors and an increase of the effective solvent viscosity with rotor density lead to a non-monotonous behaviour of the average translational velocity. Meanwhile, the rotor suspension bears a finite osmotic compressibility resulting from the long-ranged nature of hydrody- namic interactions such that rotational and odd stresses are transmitted through the solvent also at small and intermediate rotor densities. Consequently, density inhomogeneities predicted for chiral active fluids with odd viscosity can be found and allow for an explicit measurement of odd viscosity in simulations and experiments. At intermediate densities, the collective dynamics shows the emergence of multi-scale vortices and chaotic motion which is identified as active turbulence with a self-similar power-law decay in the energy spectrum, showing that the injected energy on the rotor scale is transported to larger scales, similar to the inverse energy cascade of clas- sical two-dimensional turbulence. While either odd viscosity or active turbulence have been reported in chiral active matter previously, the system studied here shows that the emergence of both simultaneously is possible resulting from the osmotic compressibility and hydrodynamic mediation of odd and active stresses. The collective dynamics of colloids rotating out of phase, i.e., where a constant torque instead of a constant angular velocity is applied, is shown to be qualitatively very similar. However, at smaller densities, local density inhomogeneities imply position dependent angular velocities of the rotors resulting from inter-rotor friction. While the friction of a quasi-2D layer of active colloids with the substrate is often not easily modifiable in experiments, the incorporation of substrate friction into the simulation models typically implies a considerable increase in computational effort. In chapter 4, a very efficient way of incorporating the friction with a substrate into a two-dimensional multiparticle collision dynamics solvent is introduced, allowing for an explicit investigation of the influences of substrate on active dynamics. For the rotor fluid, it is explicitly shown that the influence of the substrate friction results in a cutoff of the hydrodynamic interaction length, such that the maximum size of the formed vortices is controlled by the substrate friction, also resulting in a cutoff in the energy spectrum, because energy is taken out of the system at the respective length. These findings are in agreement with the experiments. Since active particles in confinement are known to organise in states of collective dynamics, ensembles of rotationally actuated colloids are studied in circular confinement and in the presence of periodic obstacle lattices in chapters 5 and 6, respectively. The results show that the chaotic active turbulent transport of rotors in suspension can be enhanced and guided resulting from edge flows generated at the boundaries, as has recently been reported for a related chiral active system. The consequent collective rotor dynamics can be regarded as a superposition of active turbulent and imposed flows, leading to on average stationary flows. In contrast to the bulk dynamics, the imposed flows inject additional energy into the system on the long length scales, and the same scaling behaviour of the energy spectrum as in bulk is only obtained if the energy injection scales, due to the mutual generation of rotor translational dynamics throughout the system and the edge flows, are well separated. The combination of edge flow and entropic layering at the boundaries leads to oscillating hydrodynamic stresses and consequently to an oscillating vorticity profile. In the presence of odd viscosity, this consequently leads to non-trivial steady-state density modulations at the boundary, resulting from a balance of osmotic pressure and odd stresses. Relevant for the efficient dispersion and mixing of inert particles on the mesoscale by means of active turbulent mixing powered by rotors, a study of the dynamics of a binary mixture consisting out of rotors and passive particles is presented in chapter 7. Because the rotors are not self-propelled, but the translational dynamics is induced by the surrounding rotors, the passive particles, which do not inject further energy into the system, are transported according to the same mechanism as the rotors. The collective dynamics thus resembles the pure rotor bulk dynamics at the respective density of only rotors. However, since no odd stresses act between the passive particles, only mutual rotor interactions lead to odd stresses leading to the accumulation of rotors in the regions of positive vorticity. This density increase is associated with a pressure increase, which balances the odd stresses acting on the rotors. However, the passive particles are only subject to the accumulation induced pressure increase such that these particles are transported into the areas of low rotor concentration, i.e., the regions of negative vorticity. Under conditions of sustained vortex flow, this results in segregation of both particle types. Since local symmetry breaking can convert injected rotational into translational energy, microswimmers can be constructed out of rotor materials when a suitable breaking of symmetry is kept in the vicinity of a rotor. One hypothetical realisation, i.e., a coupled rotor pair consisting out of two rotors of opposite angular velocity and of fixed distance, termed a birotor, are studied in chapter 8. The birotor pumps the fluid into one direction and consequently translates into the opposite direction, and creates a flow field reminiscent of a source doublet, or sliplet flow field. Fixed in space the birotor might be an interesting realisation of a microfluidic pump. The trans- lational dynamics of a birotor can be mapped onto the active Brownian particle model for single swimmers. However, due to the hydrodynamic interactions among the rotors, the birotor ensemble dynamics do not show the emergence of stable motility induced clustering. The reason for this is the flow created by birotor in small aggregates which effectively pushes further arriving birotors away from small aggregates, which eventually are all dispersed by thermal fluctuations

    Open-source high-performance software packages for direct and inverse solving of horizontal capillary flow

    Get PDF
    This work introduces Fronts, a set of open-source numerical software packages for nonlinear horizontal capillary-driven flow problems in unsaturated porous media governed by the Richards equation. The software uses the Boltzmann transformation to solve such problems in semi-infinite domains. The scheme adopted by Fronts allows it to be faster and easier to use than other tools, and provide continuous functions for all involved fields. The software is capable of solving problems that appear in hydrology, but also in other particular domains of interest such as paper-based microfluidics. As the first known open-source implementation to adopt this approach, Fronts has been validated against analytical solutions as well as existing software achieving remarkable results in terms of computational costs and numerical precision, and is meant to aid the study and modeling of capillary flow. Fronts can be freely downloaded and installed, and offers a friendly environment for new users with its complete documentation and tutorial cases.Cited as: Gerlero, G. S., Berli, C. L. A., Kler, P. A. Open-source high-performance software packages for direct and inverse solving of horizontal capillary flow. Capillarity, 2023, 6(2): 31-40. https://doi.org/10.46690/capi.2023.02.0

    Worldtube excision method for intermediate-mass-ratio inspirals: scalar-field model in 3+1 dimensions

    Full text link
    Binary black hole simulations become increasingly more computationally expensive with smaller mass ratios, partly because of the longer evolution time, and partly because the lengthscale disparity dictates smaller time steps. The program initiated by Dhesi et al. (arXiv:2109.03531) explores a method for alleviating the scale disparity in simulations with mass ratios in the intermediate astrophysical range (10−4≲q≲10−210^{-4} \lesssim q \lesssim 10^{-2}), where purely perturbative methods may not be adequate. A region ("worldtube") much larger than the small black hole is excised from the numerical domain, and replaced with an analytical model approximating a tidally deformed black hole. Here we apply this idea to a toy model of a scalar charge in a fixed circular geodesic orbit around a Schwarzschild black hole, solving for the massless Klein-Gordon field. This is a first implementation of the worldtube excision method in full 3+1 dimensions. We demonstrate the accuracy and efficiency of the method, and discuss the steps towards applying it for evolving orbits and, ultimately, in the binary black-hole scenario. Our implementation is publicly accessible in the SpECTRE numerical relativity code.Comment: 19 pages, 10 figure

    Robust Exponential Runge-Kutta Embedded Pairs

    Full text link
    Exponential integrators are explicit methods for solving ordinary differential equations that treat linear behaviour exactly. The stiff-order conditions for exponential integrators derived in a Banach space framework by Hochbruck and Ostermann are solved symbolically by expressing the Runge--Kutta weights as unknown linear combinations of phi functions. Of particular interest are embedded exponential pairs that efficiently generate both a high- and low-order estimate, allowing for dynamic adjustment of the time step. A key requirement is that the pair be robust: if the nonlinear source function has nonzero total time derivatives, the order of the low-order estimate should never exceed its design value. Robust exponential Runge--Kutta (3,2) and (4,3) embedded pairs that are well-suited to initial value problems with a dominant linearity are constructed.Comment: 24 pages, 8 figures. The Mathematica scripts mentioned in the paper can be found in: https://github.com/stiffode/expint

    Stability of space-time isogeometric methods for wave propagation problems

    Full text link
    This thesis aims at investigating the first steps toward an unconditionally stable space-time isogeometric method, based on splines of maximal regularity, for the linear acoustic wave equation. The unconditional stability of space-time discretizations for wave propagation problems is a topic of significant interest, by virtue of the advantages of space-time methods compared with more standard time-stepping techniques. In the case of continuous finite element methods, several stabilizations have been proposed. Inspired by one of these works, we address the stability issue by studying the isogeometric method for an ordinary differential equation closely related to the wave equation. As a result, we provide a stabilized isogeometric method whose effectiveness is supported by numerical tests. Motivated by these results, we conclude by suggesting an extension of this stabilization tool to the space-time isogeometric formulation of the acoustic wave equation.Comment: Masters thesi

    Path integrals and stochastic calculus

    Full text link
    Path integrals are a ubiquitous tool in theoretical physics. However, their use is sometimes hindered by the lack of control on various manipulations -- such as performing a change of the integration path -- one would like to carry out in the light-hearted fashion that physicists enjoy. Similar issues arise in the field of stochastic calculus, which we review to prepare the ground for a proper construction of path integrals. At the level of path integration, and in arbitrary space dimension, we not only report on existing Riemannian geometry-based approaches that render path integrals amenable to the standard rules of calculus, but also bring forth new routes, based on a fully time-discretized approach, that achieve the same goal. We illustrate these various definitions of path integration on simple examples such as the diffusion of a particle on a sphere.Comment: 96 pages, 4 figures. New title, expanded introduction and additional references. Version accepted in Advandes in Physic

    Effective Numerical Simulations of Synchronous Generator System

    Full text link
    Synchronous generator system is a complicated dynamical system for energy transmission, which plays an important role in modern industrial production. In this article, we propose some predictor-corrector methods and structure-preserving methods for a generator system based on the first benchmark model of subsynchronous resonance, among which the structure-preserving methods preserve a Dirac structure associated with the so-called port-Hamiltonian descriptor systems. To illustrate this, the simplified generator system in the form of index-1 differential-algebraic equations has been derived. Our analyses provide the global error estimates for a special class of structure-preserving methods called Gauss methods, which guarantee their superior performance over the PSCAD/EMTDC and the predictor-corrector methods in terms of computational stability. Numerical simulations are implemented to verify the effectiveness and advantages of our methods

    Machine Learning Methods for Autonomous Ordinary Differential Equations

    Full text link
    Ordinary Differential Equations are generally too complex to be solved analytically. Approximations thereof can be obtained by general purpose numerical methods. However, even though accurate schemes have been developed, they remain computationally expensive: In this paper, we resort to the theory of modified equations in order to obtain ''on the fly'' cheap numerical approximations. The recipe consists in approximating, prior to that, the modified field associated to the modified equation by neural networks. Elementary convergence results are then established and the efficiency of the technique is demonstrated on experiments
    • …
    corecore