Showing 1–35 of 35 results for author: Shu, C

Algebraic non-Hermitian skin effect and unified non-Bloch band theory in arbitrary dimensions

Authors: Kai Zhang, Chang Shu, Kai Sun

Abstract: The non-Hermitian skin effect, characterized by a proliferation of exponentially-localized edge modes, has led to numerous novel physical phenomena that challenge the limits of conventional band theory. In sharp contrast to the traditional exponential localization, this manuscript reports a new kind of non-Hermitian skin effect, which we term the ``algebraic non-Hermitian skin effect." This effect… ▽ More The non-Hermitian skin effect, characterized by a proliferation of exponentially-localized edge modes, has led to numerous novel physical phenomena that challenge the limits of conventional band theory. In sharp contrast to the traditional exponential localization, this manuscript reports a new kind of non-Hermitian skin effect, which we term the ``algebraic non-Hermitian skin effect." This effect emerges across a diverse spectrum of non-Hermitian systems in both two- and higher space dimensions. For 2D systems with algebraic non-Hermitian skin effect, on geometries such as a torus or cylinder, these systems exhibit behavior reminiscent of the conventional non-Hermitian skin effect, where eigenmodes are either bulk Bloch waves (on a torus) or exponentially localized edge modes (on a cylinder). However, if the same system is placed on a disk or any geometrical shape featuring open boundaries in all directions, the skin modes immediately transform into the algebraic form, with amplitude decaying as a power-law function of the distance from the boundary. To explore these novel effects, we formulate a unified generalized Brillouin zone (GBZ) framework that is universally applicable to all variations of non-Hermitian skin effects across any spatial dimension, developed through the usage of a generalized transfer-matrix approach. We find that in a $d$-dimensional non-Hermitian system, in general, the GBZ manifold's dimensionality must fall into the range from $d$ to $2d-1$, denoted by ${d \leq \dim\text{GBZ} \leq 2d-1}$. In 1D, this inequality is trivial because the upper and lower bounds converge, forcing the GBZ's dimensionality to match with that of the physical space. However, in 2D and above, this inequality indicates that there is no obligation for the GBZ's dimensionality to concur with the physical space's dimensionality, which gives rise to a new class of non-Hermitian skin effects. △ Less

Submitted 10 June, 2024; originally announced June 2024.

Comments: 19 pages, 5 figures

Explicit form of simplified Grad's 13 moments distribution function-based moment gas kinetic solver with unstructured meshes for the multiscale rarefied flow

Authors: W. Liu, C. Shu, Z. J. Liu

Abstract: It is essential to efficiently solve multiscale flows covering the continuum regime to the rarefied regime. The explicit form of Grad's 13 moments distribution function-based moment gas kinetic solver (G13-MGKS) has been proposed in our previous work [Comput. Math. Appl., 137 (2023), pp. 112-125], which demonstrates the potential for efficiently simulating continuum flows accurately and presenting… ▽ More It is essential to efficiently solve multiscale flows covering the continuum regime to the rarefied regime. The explicit form of Grad's 13 moments distribution function-based moment gas kinetic solver (G13-MGKS) has been proposed in our previous work [Comput. Math. Appl., 137 (2023), pp. 112-125], which demonstrates the potential for efficiently simulating continuum flows accurately and presenting reasonable predictions for rarefied flows at moderate Knudsen numbers on structured meshes. To further extend the solver's applicability to unstructured meshes, we propose the simplified version of the Grad's 13 moments distribution function-based moment gas kinetic solver (SG13-MGKS) with an explicit form of the numerical flux in the present paper. The Shakhov collision model has been adopted and validated within the framework of SG13-MGKS to ensure the correct Prandtl number in the simulation. Additionally, a simplified treatment for the numerical fluxes has been adopted to minimize the need for complex calculations of the gradient of integral coefficients. The performance of SG13-MGKS has been evaluated in numerical cases of Couette flow with temperature differences, flow passing through a NACA0012 airfoil, and pressure-driven flow in a variable-diameter circular pipe. Our results demonstrate that SG13-MGKS can achieve reasonably accurate computational results at Knudsen numbers below 0.2. Benefiting from the avoidance of discretization in velocity space, G13-MGKS is able to be two orders of magnitude faster compared to the conventional discrete velocity method. Furthermore, the simplified form of numerical fluxes and the fewer gradients of integration coefficients enable the performance of SG13-MGKS on unstructured grids with a saving of about 4 times the computation time and 3 times the memory cost compared to the previous version of G13-MGKS. △ Less

Submitted 24 March, 2024; originally announced March 2024.

Universal spectral moment theorem and its applications in non-Hermitian systems

Authors: Nan Cheng, Chang Shu, Kai Zhang, Xiaoming Mao, Kai Sun

Abstract: The high sensitivity of the spectrum and wavefunctions to boundary conditions, termed the non-Hermitian skin effect, represents a fundamental aspect of non-Hermitian systems. While it endows non-Hermitian systems with unprecedented physical properties, it presents notable obstacles in grasping universal properties that are robust against microscopic details and boundary conditions. In this Letter,… ▽ More The high sensitivity of the spectrum and wavefunctions to boundary conditions, termed the non-Hermitian skin effect, represents a fundamental aspect of non-Hermitian systems. While it endows non-Hermitian systems with unprecedented physical properties, it presents notable obstacles in grasping universal properties that are robust against microscopic details and boundary conditions. In this Letter, we introduce a pivotal theorem: in the thermodynamic limit, for any non-Hermitian systems with finite-range interactions, all spectral moments are invariant quantities, independent of boundary conditions, posing strong constraints on the spectrum. Utilizing this invariance, we propose a new criterion for bulk dynamical phases based on experimentally observable features and applicable to any dimensions and any boundary conditions. Based on this criterion, we define the bulk dispersive-to-proliferative phase transition, which is distinct from the real-to-complex spectral transition and contrary to traditional expectations. We verify these findings in 1D and 2D lattice models. △ Less

Submitted 15 May, 2024; v1 submitted 3 March, 2024; originally announced March 2024.

Comments: 8 Pages, 4 Figures

Efficient Alternative Finite Difference WENO Schemes for Hyperbolic Systems with Non-Conservative Products

Authors: Dinshaw S. Balsara, Deepak Bhoriya, Chi-Wang Shu, Harish Kumar

Abstract: Higher order finite difference Weighted Essentially Non-Oscillatory (WENO) schemes for conservation laws represent a technology that has been reasonably consolidated. They are extremely popular because, when applied to multidimensional problems, they offer high order accuracy at a fraction of the cost of finite volume WENO or DG schemes. They come in two flavors. There is the classical finite diff… ▽ More Higher order finite difference Weighted Essentially Non-Oscillatory (WENO) schemes for conservation laws represent a technology that has been reasonably consolidated. They are extremely popular because, when applied to multidimensional problems, they offer high order accuracy at a fraction of the cost of finite volume WENO or DG schemes. They come in two flavors. There is the classical finite difference WENO (FD-WENO) method (Shu and Osher, J. Comput. Phys., 83 (1989) 32-78). However, in recent years there is also an alternative finite difference WENO (AFD-WENO) method which has recently been formalized into a very useful general-purpose algorithm for conservation laws (Balsara et al., Efficient Alternative Finite Difference WENO Schemes for Hyperbolic Conservation Laws, submitted to CAMC (2023)). However, the FD-WENO algorithm has only very recently been formulated for hyperbolic systems with non-conservative products (Balsara et al., Efficient Finite Difference WENO Scheme for Hyperbolic Systems with Non-Conservative Products, to appear CAMC (2023)). In this paper we show that there are substantial advantages in obtaining an AFD-WENO algorithm for hyperbolic systems with non-conservative products. Such an algorithm is documented in this paper. We present an AFD-WENO formulation in fluctuation form that is carefully engineered to retrieve the flux form when that is warranted and nevertheless extends to non-conservative products. The method is flexible because it allows any Riemann solver to be used. The formulation we arrive at is such that when non-conservative products are absent it reverts exactly to the formulation in the second citation above which is in exact flux conservation form. The ability to transition to a precise conservation form when non-conservative products are absent ensures, via the Lax-Wendroff theorem, that shock locations will be exactly ... △ Less

Submitted 2 March, 2024; originally announced March 2024.

Comments: Accepted in Communications on Applied Mathematics and Computation

MSC Class: 65M06; 65M22; 65N06; 35Q35; 76T10; 76M20

Efficient Alternative Finite Difference WENO Schemes for Hyperbolic Conservation Laws

Authors: Dinshaw S. Balsara, Deepak Bhoriya, Chi-Wang Shu, Harish Kumar

Abstract: Higher order finite difference Weighted Essentially Non-Oscillatory (WENO) schemes for conservation laws are extremely popular because, for multidimensional problems, they offer high order accuracy at a fraction of the cost of finite volume WENO or DG schemes. Such schemes come in two formulations. The very popular classical finite difference WENO (FD-WENO) method (Shu and Osher, J. Comput. Phys.,… ▽ More Higher order finite difference Weighted Essentially Non-Oscillatory (WENO) schemes for conservation laws are extremely popular because, for multidimensional problems, they offer high order accuracy at a fraction of the cost of finite volume WENO or DG schemes. Such schemes come in two formulations. The very popular classical finite difference WENO (FD-WENO) method (Shu and Osher, J. Comput. Phys., 83 (1989) 32-78) relies two reconstruction steps applied to two split fluxes. However, the method cannot accommodate different types of Riemann solvers and cannot preserve free stream boundary conditions on curvilinear meshes. This limits its utility. The alternative finite difference WENO (AFD-WENO) method can overcome these deficiencies, however, much less work has been done on this method. The reasons are three-fold. First, it is difficult for the casual reader to understand the intricate logic that requires higher order derivatives of the fluxes to be evaluated at zone boundaries. The analytical methods for deriving the update equation for AFD-WENO schemes are somewhat recondite. To overcome that difficulty, we provide an easily accessible script that is based on a computer algebra system in Appendix A of this paper. Second, the method relies on interpolation rather than reconstruction, and WENO interpolation formulae have not been documented in the literature as thoroughly as WENO reconstruction formulae. In this paper, we explicitly provide all necessary WENO interpolation formulae that are needed for implementing AFD-WENO up to ninth order. The third reason is that AFD-WENO requires higher order derivatives of the fluxes to be available at zone boundaries. Since those derivatives are usually obtained by finite differencing the zone-centered fluxes, they become susceptible to a Gibbs phenomenon when the solution ... △ Less

Submitted 2 March, 2024; originally announced March 2024.

Comments: Accepted in Communications on Applied Mathematics and Computation

MSC Class: 65M06; 65M22; 65N06; 35Q35; 76T10; 76M20

Increased Atom-Cavity Coupling through Cooling-Induced Atomic Reorganization

Authors: Chi Shu, Simone Colombo, Zeyang Li, Albert Adiyatullin, Enrique Mendez, Edwin Pedrozo-Peñafiel, Vladan Vuletić

Abstract: The strong coupling of atoms to optical cavities can improve optical lattice clocks as the cavity enables metrologically useful collective atomic entanglement and high-fidelity measurement. To this end, it is necessary to cool the ensemble to suppress motional broadening, and advantageous to maximize and homogenize the atom-cavity coupling. We demonstrate resolved Raman sideband cooling via the ca… ▽ More The strong coupling of atoms to optical cavities can improve optical lattice clocks as the cavity enables metrologically useful collective atomic entanglement and high-fidelity measurement. To this end, it is necessary to cool the ensemble to suppress motional broadening, and advantageous to maximize and homogenize the atom-cavity coupling. We demonstrate resolved Raman sideband cooling via the cavity as a method that can simultaneously achieve both goals. In 200 ms, we cool 171Yb atoms to an average vibration number <nx> = 0.23(7) in the tightly binding direction, resulting in 93% optical π-pulse fidelity on the clock transition 1S0 -> 3P0. During cooling, the atoms self-organize into locations with maximal atom-cavity-coupling, which will improve quantum metrology applications. △ Less

Submitted 15 October, 2023; originally announced October 2023.

Comments: 5 + 2 pages, 5 + 2 figures

A variational quantum algorithm-based numerical method for solving potential and Stokes flows

Authors: Yangyang Liu, Zhen Chen, Chang Shu, Patrick Rebentrost, Yaguang Liu, S. C. Chew, B. C. Khoo, Y. D. Cui

Abstract: This paper presents a numerical method based on the variational quantum algorithm to solve potential and Stokes flow problems. In this method, the governing equations for potential and Stokes flows can be respectively written in the form of Laplace's equation and Stokes equations using velocity potential, stream function and vorticity formulations. Then the finite difference method and the general… ▽ More This paper presents a numerical method based on the variational quantum algorithm to solve potential and Stokes flow problems. In this method, the governing equations for potential and Stokes flows can be respectively written in the form of Laplace's equation and Stokes equations using velocity potential, stream function and vorticity formulations. Then the finite difference method and the generalised differential quadrature (GDQ) method are applied to discretize the governing equations. For the prescribed boundary conditions, the corresponding linear systems of equations can be obtained. These linear systems are solved by using the variational quantum linear solver (VQLS), which resolves the potential and Stokes flow problems equivalently. To the best of authors' knowledge, this is the first study that incorporates the GDQ method which is inherently a high-order discretization method with the VQLS algorithm. Since the GDQ method can utilize much fewer grid points than the finite difference method to approximate derivatives with a higher order of accuracy, the size of the input matrix for the VQLS algorithm can be smaller. In this way, the computational cost may be saved. The performance of the present method is comprehensively assessed by two representative examples, namely, the potential flow around a circular cylinder and Stokes flow in a lid-driven cavity. Numerical results validate the applicability and accuracy of the present VQLS-based method. Furthermore, its time complexity is evaluated by the heuristic scaling, which demonstrates that the present method scales efficiently in the number of qubits and the precision. This work brings quantum computing to the field of computational fluid dynamics. By virtue of quantum advantage over classical methods, promising advances in solving large-scale fluid mechanics problems of engineering interest may be prompted. △ Less

Submitted 3 March, 2023; originally announced March 2023.

Improving Metrology with Quantum Scrambling

Authors: Zeyang Li, Simone Colombo, Chi Shu, Gustavo Velez, Saúl Pilatowsky-Cameo, Roman Schmied, Soonwon Choi, Mikhail Lukin, Edwin Pedrozo-Peñafiel, Vladan Vuletić

Abstract: Quantum scrambling describes the spreading of local information into many degrees of freedom in quantum systems. This provides the conceptual connection among diverse phenomena ranging from thermalizing quantum dynamics to models of black holes. Here we experimentally probe the exponential scrambling of a multi-particle system near a bistable point in phase space and utilize it for entanglement-en… ▽ More Quantum scrambling describes the spreading of local information into many degrees of freedom in quantum systems. This provides the conceptual connection among diverse phenomena ranging from thermalizing quantum dynamics to models of black holes. Here we experimentally probe the exponential scrambling of a multi-particle system near a bistable point in phase space and utilize it for entanglement-enhanced metrology. We use a time-reversal protocol to observe a simultaneous exponential growth of both the metrological gain and the out-of-time-order correlator, thereby experimentally verifying the relation between quantum metrology and quantum information scrambling. Our experiments demonstrate that fast-scrambling dynamics capable of exponentially fast entanglement generation are useful for practical metrology, resulting in 6.8(4) dB gain beyond the Standard Quantum Limit. △ Less

Submitted 4 February, 2023; v1 submitted 24 December, 2022; originally announced December 2022.

Comments: 5 figures

Journal ref: Science 380, 1381-1384 (2023)

Grad's Distribution Function for 13 Moments based Moment Gas Kinetic Solver for Steady and Unsteady Rarefied flows: Discrete and Explicit Forms

Authors: W. Liu, Z. J. Liu, Z. L. Zhang, C. J. Teo, C. Shu

Abstract: Efficient modeling of rarefied flow has drawn widespread interest for practical engineering applications. In the present work, we proposed the Grad's distribution function for 13 moments-based moment gas kinetic solver (G13-MGKS) and the macroscopic governing equations are derived based on the moment integral of discrete Boltzmann equation in the finite volume framework. Numerical fluxes at the ce… ▽ More Efficient modeling of rarefied flow has drawn widespread interest for practical engineering applications. In the present work, we proposed the Grad's distribution function for 13 moments-based moment gas kinetic solver (G13-MGKS) and the macroscopic governing equations are derived based on the moment integral of discrete Boltzmann equation in the finite volume framework. Numerical fluxes at the cell interface related to the macroscopic variables, stress and heat flux can be reconstructed from the Boltzmann integration equation at surrounding points of the cell interface directly, so the complicated partial differential equations with tedious implementation of boundary conditions in the moment method can be avoided. Meanwhile, the explicit expression of numerical fluxes is proposed, which could release the present solver the from the discretization and numerical summation in molecular velocity space. To evaluate the Grad's distribution function for 13 moments in the present framework, the G13-MGKS with the discrete and explicit form of numerical fluxes are examined by several test cases covering the steady and unsteady rarefied flows. Numerical results indicate that the G13-MGKS could simulate continuum flows accurately and present reasonable prediction for rarefied flows at moderate Knudsen number. Moreover, the tests of computations and memory costs demonstrate that the present framework could preserve the highly efficient feature. △ Less

Submitted 19 July, 2022; originally announced July 2022.

Accurate near wall steady flow field prediction using Physics Informed Neural Network (PINN)

Authors: Vinothkumar Sekar, Qinghua Jiang, Chang Shu, Boo Cheong Khoo

Abstract: In this paper, Physics Informed Neural Network (PINN) is explored in order to obtain flow predictions near the wall region accurately with measurements (or sampling points) away from the wall. Often, in fluid mechanics experiments, it is difficult to perform velocity measurements near the wall accurately. Therefore, the present study reveals a new and elegant approach to recover the flow solutions… ▽ More In this paper, Physics Informed Neural Network (PINN) is explored in order to obtain flow predictions near the wall region accurately with measurements (or sampling points) away from the wall. Often, in fluid mechanics experiments, it is difficult to perform velocity measurements near the wall accurately. Therefore, the present study reveals a new and elegant approach to recover the flow solutions near the wall. Laminar boundary layer flow over a flat plate case is considered for this study in order to explore the ability of PINN to accurately predict the flow field. All the required sampling data for this study is obtained from CFD simulations. A wide range of Reynolds number cases from Re=500 to 100000 has been investigated. First, using PINN, the boundary layer solution is obtained with three different types of boundary conditions. Further, the influence of the location of the sampling points on the accuracy is analysed. From the velocity profiles and the skin friction coefficient distribution, it is clear that PINN results are reasonably accurate near the wall with only a few sampling points away from the wall. This approach has potential application in experiments to obtain the near wall solutions accurately with measurements away from the wall. △ Less

Submitted 7 April, 2022; originally announced April 2022.

Provably Positive Central DG Schemes via Geometric Quasilinearization for Ideal MHD Equations

Authors: Kailiang Wu, Haili Jiang, Chi-Wang Shu

Abstract: In the numerical simulation of ideal MHD, keeping the pressure and density positive is essential for both physical considerations and numerical stability. This is a challenge, due to the underlying relation between such positivity-preserving (PP) property and the magnetic divergence-free (DF) constraint as well as the strong nonlinearity of the MHD equations. This paper presents the first rigorous… ▽ More In the numerical simulation of ideal MHD, keeping the pressure and density positive is essential for both physical considerations and numerical stability. This is a challenge, due to the underlying relation between such positivity-preserving (PP) property and the magnetic divergence-free (DF) constraint as well as the strong nonlinearity of the MHD equations. This paper presents the first rigorous PP analysis of the central discontinuous Galerkin (CDG) methods and constructs arbitrarily high-order PP CDG schemes for ideal MHD. By the recently developed geometric quasilinearization (GQL) approach, our analysis reveals that the PP property of standard CDG methods is closely related to a discrete DF condition, whose form was unknown and differs from the non-central DG and finite volume cases in [K. Wu, SIAM J. Numer. Anal. 2018]. This result lays the foundation for the design of our PP CDG schemes. In 1D case, the discrete DF condition is naturally satisfied, and we prove the standard CDG method is PP under a condition that can be enforced with a PP limiter. However, in the multidimensional cases, the discrete DF condition is highly nontrivial yet critical, and we prove the the standard CDG method, even with the PP limiter, is not PP in general, as it fails to meet the discrete DF condition. We address this issue by carefully analyzing the structure of the discrete divergence and then constructing new locally DF CDG schemes for Godunov's modified MHD equations with an additional source. The key point is to find out the suitable discretization of the source term such that it exactly offsets all the terms in the discrete DF condition. Based on the GQL approach, we prove the PP property of the new multidimensional CDG schemes. The robustness and accuracy of PP CDG schemes are validated by several demanding examples, including the high-speed jets and blast problems with very low plasma beta. △ Less

Submitted 28 March, 2022; originally announced March 2022.

Geometric Quasilinearization Framework for Analysis and Design of Bound-Preserving Schemes

Authors: Kailiang Wu, Chi-Wang Shu

Abstract: Solutions to many partial differential equations satisfy certain bounds or constraints. For example, the density and pressure are positive for equations of fluid dynamics, and in the relativistic case the fluid velocity is upper bounded by the speed of light, etc. As widely realized, it is crucial to develop bound-preserving numerical methods that preserve such intrinsic constraints. Exploring pro… ▽ More Solutions to many partial differential equations satisfy certain bounds or constraints. For example, the density and pressure are positive for equations of fluid dynamics, and in the relativistic case the fluid velocity is upper bounded by the speed of light, etc. As widely realized, it is crucial to develop bound-preserving numerical methods that preserve such intrinsic constraints. Exploring provably bound-preserving schemes has attracted much attention and is actively studied in recent years. This is however still a challenging task for many systems especially those involving nonlinear constraints. Based on some key insights from geometry, we systematically propose an innovative and general framework, referred to as geometric quasilinearization (GQL), which paves a new effective way for studying bound-preserving problems with nonlinear constraints. The essential idea of GQL is to equivalently transfer all nonlinear constraints into linear ones, through properly introducing some free auxiliary variables. We establish the fundamental principle and general theory of GQL via the geometric properties of convex regions, and propose three simple effective methods for constructing GQL. We apply the GQL approach to a variety of partial differential equations, and demonstrate its effectiveness and remarkable advantages for studying bound-preserving schemes, by diverse challenging examples and applications which cannot be easily handled by direct or traditional approaches. △ Less

Submitted 8 November, 2021; originally announced November 2021.

Collective Spin-Light and Light-Mediated Spin-Spin Interactions in an Optical Cavity

Authors: Zeyang Li, Boris Braverman, Simone Colombo, Chi Shu, Akio Kawasaki, Albert Adiyatullin, Edwin Pedrozo-Peñafiel, Enrique Mendez, Vladan Vuletić

Abstract: The interaction between an atomic ensemble and a light mode in a high-finesse optical cavity can easily reach the strong-coupling regime, where quantum effects dominate. In this regime, the interaction can be used to generate both atom-light and atom-atom entanglement. We analyze the dominant effects on the collective atomic state and the light field, and derive a unified approach that can account… ▽ More The interaction between an atomic ensemble and a light mode in a high-finesse optical cavity can easily reach the strong-coupling regime, where quantum effects dominate. In this regime, the interaction can be used to generate both atom-light and atom-atom entanglement. We analyze the dominant effects on the collective atomic state and the light field, and derive a unified approach that can account for atomic entanglement induced both by measurements on the light field, and by ignoring the state of the light field altogether. We present analytical expressions for the entanglement induced by the interaction, and determine the conditions that maximize the entanglement-induced gain over the standard quantum limit in quantum sensors and atomic clocks. △ Less

Submitted 25 March, 2022; v1 submitted 22 June, 2021; originally announced June 2021.

Comments: 25 pages, accepted for publication in PRX Quantum

Journal ref: PRX Quantum 3, 020308, 2022

Time-Reversal-Based Quantum Metrology with Many-Body Entangled States

Authors: Simone Colombo, Edwin Pedrozo-Peñafiel, Albert F. Adiyatullin, Zeyang Li, Enrique Mendez, Chi Shu, Vladan Vuletic

Abstract: In quantum metrology, entanglement represents a valuable resource that can be used to overcome the Standard Quantum Limit (SQL) that bounds the precision of sensors that operate with independent particles. Measurements beyond the SQL are typically enabled by relatively simple entangled states (squeezed states with Gaussian probability distributions), where quantum noise is redistributed between di… ▽ More In quantum metrology, entanglement represents a valuable resource that can be used to overcome the Standard Quantum Limit (SQL) that bounds the precision of sensors that operate with independent particles. Measurements beyond the SQL are typically enabled by relatively simple entangled states (squeezed states with Gaussian probability distributions), where quantum noise is redistributed between different quadratures. However, due to both fundamental limitations and the finite measurement resolution achieved in practice, sensors based on squeezed states typically operate far from the true fundamental limit of quantum metrology, the Heisenberg Limit. Here, by implementing an effective time-reversal protocol through a controlled sign change in an optically engineered many-body Hamiltonian, we demonstrate atomic-sensor performance with non-Gaussian states beyond the limitations of spin squeezing, and without the requirement of extreme measurement resolution. Using a system of 350 neutral $^{171}$Yb atoms, this signal amplification through time-reversed interaction (SATIN) protocol achieves the largest sensitivity improvement beyond the SQL ($11.8 \pm 0.5$~dB) demonstrated in any interferometer to date. Furthermore, we demonstrate a precision improving in proportion to the particle number (Heisenberg scaling), at fixed distance of 12.6~dB from the Heisenberg Limit. These results pave the way for quantum metrology using complex entangled states, with potential broad impact in science and technology. Potential applications include searches for dark matter and for physics beyond the standard model, tests of the fundamental laws of physics, timekeeping, and geodesy. △ Less

Submitted 27 September, 2021; v1 submitted 7 June, 2021; originally announced June 2021.

On the evolution of fuel droplet evaporation zone and its interaction with the flame front in ignition of spray flames

Authors: Qiang Li, Chang Shu, Huangwei Zhang

Abstract: Evolution of fuel droplet evaporation zone and its interaction with the propagating flame front are studied in this work. A general theory is developed to describe the evolutions of flame propagation speed, flame temperature, droplet evaporation onset and completion locations in ignition and propagation of spherical flames. The influences of liquid droplet mass loading, heat exchange coefficient (… ▽ More Evolution of fuel droplet evaporation zone and its interaction with the propagating flame front are studied in this work. A general theory is developed to describe the evolutions of flame propagation speed, flame temperature, droplet evaporation onset and completion locations in ignition and propagation of spherical flames. The influences of liquid droplet mass loading, heat exchange coefficient (or evaporation rate) and Lewis number on spherical spray flame ignition are studied. Two flame regimes are considered, i.e., heterogeneous and homogeneous flames, based on the mixture condition near the flame front. The results indicate that the spray flame trajectories are considerably affected by the ignition energy addition. The critical condition for successful ignition for the fuel-rich mixture is coincidence of inner and outer flame balls from igniting kernel and propagating flame. The flame balls always exist in homogeneous mixtures, indicating that ignition failure and critical successful events occur only in purely gaseous mixture. The fuel droplets have limited effects on minimum ignition energy, which however increases monotonically with the Lewis number. Moreover, flame kernel originates from heterogeneous mixtures due to the initially dispersed droplets near the spark. The evaporative heat loss in the burned and unburned zones of homogeneous and heterogeneous spray flames is also evaluated, and the results show that for the failed flame kernels, evaporative heat loss behind and before the flame front first increases and then decreases. The evaporative heat loss before the flame front generally increases, although non-monotonicity exists, when the flame is successfully ignited and propagate outwardly. For heterogeneous flames, the ratio of the heat loss from the burned zone to the total one decreases as the flame expands. △ Less

Submitted 25 March, 2021; originally announced March 2021.

Propagation of weakly stretched premixed spherical spray flames in localized homogeneous and heterogeneous reactants

Authors: Qiang Li, Huangwei Zhang, Chang Shu

Abstract: Propagation of weakly stretched spherical flames in partially pre-vaporized fuel sprays is theoretically investigated in this work. A general theory is developed to describe flame propagation speed, flame temperature, droplet evaporation onset and completion locations. The influences of liquid fuel and gas mixture properties on spherical spray flame propagation are studied. The results indicate th… ▽ More Propagation of weakly stretched spherical flames in partially pre-vaporized fuel sprays is theoretically investigated in this work. A general theory is developed to describe flame propagation speed, flame temperature, droplet evaporation onset and completion locations. The influences of liquid fuel and gas mixture properties on spherical spray flame propagation are studied. The results indicate that the spray flame propagation speed is enhanced with increased droplet mass loading and/or evaporation heat exchange coefficient (or evaporation rate). Opposite trends are found when the latent heat is high, due to strong evaporation heat absorption. Fuel vapor and temperature gradients are observed in the post-flame evaporation zone of heterogeneous flames. Evaporation completion front location considerably changes with flame radius, but the evaporation onset location varies little relative to the flame front when the flame propagates. For larger droplet loading and smaller evaporation rate, the fuel droplet tends to complete evaporation behind the flame front. Flame bifurcation occurs with high droplet mass loading under large latent heat, leading to multiplicity of flame propagation speed, droplet evaporation onset and completion fronts. The flame enhancement or weakening effects by the fuel droplet sprays are revealed by enhanced or suppressed heat and mass diffusion process in the pre-flame zone. Besides, for heterogeneous flames, heat and mass diffusion in the post-flame zone also exists. The mass diffusion for both homogeneous and heterogeneous flames is enhanced with decreased Lewis number. The magnitude of Markstein length is considerably reduced with increased droplet loading. Moreover, post-flame droplet burning behind heterogeneous flame influences the flame propagation speed and Markstein length when the liquid fuel loading is relatively low. △ Less

Submitted 27 September, 2020; originally announced September 2020.

Entanglement-Enhanced Optical Atomic Clock

Authors: Edwin Pedrozo-Peñafiel, Simone Colombo, Chi Shu, Albert F. Adiyatullin, Zeyang Li, Enrique Mendez, Boris Braverman, Akio Kawasaki, Daisuke Akamatsu, Yanhong Xiao, Vladan Vuletić

Abstract: State-of-the-art atomic clocks are based on the precise detection of the energy difference between two atomic levels, measured as a quantum phase accumulated in a given time interval. Optical-lattice clocks (OLCs) now operate at or near the standard quantum limit (SQL) that arises from the quantum noise associated with discrete measurement outcomes. While performance beyond the SQL has been achiev… ▽ More State-of-the-art atomic clocks are based on the precise detection of the energy difference between two atomic levels, measured as a quantum phase accumulated in a given time interval. Optical-lattice clocks (OLCs) now operate at or near the standard quantum limit (SQL) that arises from the quantum noise associated with discrete measurement outcomes. While performance beyond the SQL has been achieved in microwave clocks and other atomic sensors by engineering quantum correlations (entanglement) between the atoms, the generation of entanglement on an optical-clock transition and operation of such a clock beyond the SQL represent major goals in quantum metrology that have never been demonstrated. Here we report creation of a many-atom entangled state on an optical transition, and demonstrate an OLC with an Allan deviation below the SQL. We report a metrological gain of $4.4^{+0.6}_{-0.4}$ dB over the SQL using an ensemble consisting of a few hundred 171Yb atoms, allowing us to reach a given stability $2.8{\pm}0.3$ times faster than the same clock operated at the SQL. Our results should be readily applicable to other systems, thus enabling further advances in timekeeping precision and accuracy. Entanglement-enhanced OLCs will have many scientific and technological applications, including precision tests of the fundamental laws of physics, geodesy, or gravitational wave detection. △ Less

Submitted 16 June, 2020; v1 submitted 12 June, 2020; originally announced June 2020.

Comments: 4 figures, 7 pages

Journal ref: Nature 588, 414 418 (2020)

Trapping $^{171}$Yb Atoms into a One-Dimensional Optical Lattice with a Small Waist

Authors: Akio Kawasaki, Boris Braverman, Edwin Pedrozo-Peñafiel, Chi Shu, Simone Colombo, Zeyang Li, Vladan Vuletić

Abstract: In most experiments with atoms trapped in optical lattices, the transverse size of the optical lattice beams is on the order of tens of micrometers, and loading many atoms into smaller optical lattices has not been carefully investigated. We report trapping 1500 $^{171}$Yb atoms in a one-dimensional optical lattice generated by a narrow cavity mode at a distance of 0.14 mm from a mirror surface. T… ▽ More In most experiments with atoms trapped in optical lattices, the transverse size of the optical lattice beams is on the order of tens of micrometers, and loading many atoms into smaller optical lattices has not been carefully investigated. We report trapping 1500 $^{171}$Yb atoms in a one-dimensional optical lattice generated by a narrow cavity mode at a distance of 0.14 mm from a mirror surface. The simplest approach of loading atoms from a mirror magneto-optical trap overlapped with the cavity mode allows the adjustment of the loading position by tuning a uniform bias magnetic field. The number of atoms trapped in the optical lattice exhibits two local maxima for different lattice depths, with a global maximum in the deeper lattice. These results open a way to quantum mechanical manipulation of atoms based on strong interaction with a tightly focused light field. △ Less

Submitted 6 August, 2020; v1 submitted 23 February, 2020; originally announced February 2020.

Comments: 9 pages, 8 figures

Journal ref: Phys. Rev. A 102, 013114 (2020)

Provably Physical-Constraint-Preserving Discontinuous Galerkin Methods for Multidimensional Relativistic MHD Equations

Authors: Kailiang Wu, Chi-Wang Shu

Abstract: We propose and analyze a class of robust, uniformly high-order accurate discontinuous Galerkin (DG) schemes for multidimensional relativistic magnetohydrodynamics (RMHD) on general meshes. A distinct feature of the schemes is their physical-constraint-preserving (PCP) property, i.e., they are proven to preserve the subluminal constraint on the fluid velocity and the positivity of density, pressure… ▽ More We propose and analyze a class of robust, uniformly high-order accurate discontinuous Galerkin (DG) schemes for multidimensional relativistic magnetohydrodynamics (RMHD) on general meshes. A distinct feature of the schemes is their physical-constraint-preserving (PCP) property, i.e., they are proven to preserve the subluminal constraint on the fluid velocity and the positivity of density, pressure, and specific internal energy. Developing PCP high-order schemes for RMHD is highly desirable but remains a challenging task, especially in the multidimensional cases, due to the inherent strong nonlinearity in the constraints and the effect of the magnetic divergence-free condition. Inspired by some crucial observations at the PDE level, we construct the provably PCP schemes by using the locally divergence-free DG schemes of the recently proposed symmetrizable RMHD equations as the base schemes, a limiting technique to enforce the PCP property of the DG solutions, and the strong-stability-preserving methods for time discretization. We rigorously prove the PCP property by using a novel "quasi-linearization" approach to handle the highly nonlinear physical constraints, technical splitting to offset the influence of divergence error, and sophisticated estimates to analyze the beneficial effect of the additional source term in the symmetrizable RMHD system. Several two-dimensional numerical examples are provided to confirm the PCP property and to demonstrate the accuracy, effectiveness and robustness of the proposed PCP schemes. △ Less

Submitted 9 February, 2020; originally announced February 2020.

Entropy Symmetrization and High-Order Accurate Entropy Stable Numerical Schemes for Relativistic MHD Equations

Authors: Kailiang Wu, Chi-Wang Shu

Abstract: This paper presents entropy symmetrization and high-order accurate entropy stable schemes for the relativistic magnetohydrodynamic (RMHD) equations. It is shown that the conservative RMHD equations are not symmetrizable and do not admit a thermodynamic entropy pair. To address this issue, a symmetrizable RMHD system, equipped with a convex thermodynamic entropy pair, is proposed by adding a source… ▽ More This paper presents entropy symmetrization and high-order accurate entropy stable schemes for the relativistic magnetohydrodynamic (RMHD) equations. It is shown that the conservative RMHD equations are not symmetrizable and do not admit a thermodynamic entropy pair. To address this issue, a symmetrizable RMHD system, equipped with a convex thermodynamic entropy pair, is proposed by adding a source term into the equations, providing an analogue to the non-relativistic Godunov--Powell system. Arbitrarily high-order accurate entropy stable finite difference schemes are developed on Cartesian meshes based on the symmetrizable RMHD system. The crucial ingredients of these schemes include (i) affordable explicit entropy conservative fluxes which are technically derived through carefully selected parameter variables, (ii) a special high-order discretization of the source term in the symmetrizable RMHD system, and (iii) suitable high-order dissipative operators based on essentially non-oscillatory reconstruction to ensure the entropy stability. Several numerical tests demonstrate the accuracy and robustness of the proposed entropy stable schemes. △ Less

Submitted 9 July, 2021; v1 submitted 17 July, 2019; originally announced July 2019.

Comments: 32 pages; 11 figures; add some references on entropy symmetrization

Journal ref: SIAM Journal on Scientific Computing, 42(4): A2230-A2261, 2020

Near-Unitary Spin Squeezing in $^{171}$Yb

Authors: Boris Braverman, Akio Kawasaki, Edwin Pedrozo-Peñafiel, Simone Colombo, Chi Shu, Zeyang Li, Enrique Mendez, Megan Yamoah, Leonardo Salvi, Daisuke Akamatsu, Yanhong Xiao, Vladan Vuletić

Abstract: Spin squeezing can improve atomic precision measurements beyond the standard quantum limit (SQL), and unitary spin squeezing is essential for improving atomic clocks. We report substantial and nearly unitary spin squeezing in $^{171}$Yb, an optical lattice clock atom. The collective nuclear spin of $\sim 10^3$ atoms is squeezed by cavity feedback, using light detuned from the system's resonances t… ▽ More Spin squeezing can improve atomic precision measurements beyond the standard quantum limit (SQL), and unitary spin squeezing is essential for improving atomic clocks. We report substantial and nearly unitary spin squeezing in $^{171}$Yb, an optical lattice clock atom. The collective nuclear spin of $\sim 10^3$ atoms is squeezed by cavity feedback, using light detuned from the system's resonances to attain unitarity. The observed precision gain over the SQL is limited by state readout to 6.5(4) dB, while the generated states offer a gain of 12.9(6) dB, limited by the curvature of the Bloch sphere. Using a squeezed state within 30% of unitarity, we demonstrate an interferometer that improves the averaging time over the SQL by a factor of 3.7(2). In the future, the squeezing can be simply transferred onto the optical clock transition of $^{171}$Yb. △ Less

Submitted 9 May, 2019; v1 submitted 29 January, 2019; originally announced January 2019.

Comments: 5 pages, 4 figures

Journal ref: Phys. Rev. Lett. 122, 223203 (2019)

Well-balanced finite volume schemes for hydrodynamic equations with general free energy

Authors: José A. Carrillo, Serafim Kalliadasis, Sergio P. Perez, Chi-Wang Shu

Abstract: Well balanced and free energy dissipative first- and second-order accurate finite volume schemes are proposed for a general class of hydrodynamic systems with linear and nonlinear damping. The natural Liapunov functional of the system, given by its free energy, allows for a characterization of the stationary states by its variation. An analog property at the discrete level enables us to preserve s… ▽ More Well balanced and free energy dissipative first- and second-order accurate finite volume schemes are proposed for a general class of hydrodynamic systems with linear and nonlinear damping. The natural Liapunov functional of the system, given by its free energy, allows for a characterization of the stationary states by its variation. An analog property at the discrete level enables us to preserve stationary states at machine precision while keeping the dissipation of the discrete free energy. These schemes allow for analysing accurately the stability properties of stationary states in challeging problems such as: phase transitions in collective behavior, generalized Euler-Poisson systems in chemotaxis and astrophysics, and models in dynamic density functional theories; having done a careful validation in a battery of relevant test cases. △ Less

Submitted 4 November, 2020; v1 submitted 3 December, 2018; originally announced December 2018.

Comments: Videos from the simulations of this work are available at https://figshare.com/projects/Well-balanced_finite_volume_schemes_for_hydrodynamic_equations_with_general_free_energy/60122

Geometrically asymmetric optical cavity for strong atom-photon coupling

Authors: Akio Kawasaki, Boris Braverman, Edwin Pedrozo-Peñafiel, Chi Shu, Simone Colombo, Zeyang Li, Özge Özel, Wenlan Chen, Leonardo Salvi, André Heinz, David Levonian, Daisuke Akamatsu, Yanhong Xiao, Vladan Vuletić

Abstract: Optical cavities are widely used to enhance the interaction between atoms and light. Typical designs using a geometrically symmetric structure in the near-concentric regime face a tradeoff between mechanical stability and high single-atom cooperativity. To overcome this limitation, we design and implement a geometrically asymmetric standing-wave cavity. This structure, with mirrors of very differe… ▽ More Optical cavities are widely used to enhance the interaction between atoms and light. Typical designs using a geometrically symmetric structure in the near-concentric regime face a tradeoff between mechanical stability and high single-atom cooperativity. To overcome this limitation, we design and implement a geometrically asymmetric standing-wave cavity. This structure, with mirrors of very different radii of curvature, allows strong atom-light coupling while exhibiting good stability against misalignment. We observe effective cooperativities ranging from $η_{\rm eff}=10$ to $η_{\rm eff}=0.2$ by shifting the location of the atoms in the cavity mode. By loading $^{171}$Yb atoms directly from a mirror magneto-optical trap into a one-dimensional optical lattice along the cavity mode, we produce atomic ensembles with collective cooperativities up to $Nη=2\times 10^4$. This system opens a way to preparing spin squeezing for an optical lattice clock and to accessing a range of nonclassical collective states. △ Less

Submitted 5 February, 2019; v1 submitted 20 November, 2018; originally announced November 2018.

Comments: 8 pages, 7 figures, published version

Journal ref: Phys. Rev. A 99, 013437 (2019)

Novel development of dissipative-soliton-resonance pulses with pump power in an all-normal-dispersion fiber laser

Authors: Yufei Wang, Lei Li, Shuai Wang, Liming Hua, Chaojie Shu, Lei Su, D. Y. Tang, D. Y. Shen, Luming Zhao

Abstract: Evolution of dissipative-soliton-resonance (DSR) pulses in an all-normal-dispersion fiber laser with pump power under various operation conditions are experimentally studied. The fiber laser is mode-locked by using a nonlinear fiber loop mirror. Apart from the typical pulse broadening due to the peak power clamping effect, pulse breaking was observed. In addition, pulse narrowing with pump power i… ▽ More Evolution of dissipative-soliton-resonance (DSR) pulses in an all-normal-dispersion fiber laser with pump power under various operation conditions are experimentally studied. The fiber laser is mode-locked by using a nonlinear fiber loop mirror. Apart from the typical pulse broadening due to the peak power clamping effect, pulse breaking was observed. In addition, pulse narrowing with pump power increasing was observed. All the incompatible evolution of DSR pulses could be attributed to the multiple parameters change due to pump power increasing under specific operation conditions. △ Less

Submitted 7 October, 2018; originally announced October 2018.

Comments: 5 pages, 7 figures

Provably Positive High-Order Schemes for Ideal Magnetohydrodynamics: Analysis on General Meshes

Authors: Kailiang Wu, Chi-Wang Shu

Abstract: This paper proposes and analyzes arbitrarily high-order discontinuous Galerkin (DG) and finite volume methods which provably preserve the positivity of density and pressure for the ideal MHD on general meshes. Unified auxiliary theories are built for rigorously analyzing the positivity-preserving (PP) property of MHD schemes with a HLL type flux on polytopal meshes in any space dimension. The main… ▽ More This paper proposes and analyzes arbitrarily high-order discontinuous Galerkin (DG) and finite volume methods which provably preserve the positivity of density and pressure for the ideal MHD on general meshes. Unified auxiliary theories are built for rigorously analyzing the positivity-preserving (PP) property of MHD schemes with a HLL type flux on polytopal meshes in any space dimension. The main challenges overcome here include establishing relation between the PP property and discrete divergence of magnetic field on general meshes, and estimating proper wave speeds in the HLL flux to ensure the PP property. In 1D case, we prove that the standard DG and finite volume methods with the proposed HLL flux are PP, under condition accessible by a PP limiter. For multidimensional conservative MHD system, standard DG methods with a PP limiter are not PP in general, due to the effect of unavoidable divergence-error. We construct provably PP high-order DG and finite volume schemes by proper discretization of symmetrizable MHD system, with two divergence-controlling techniques: locally divergence-free elements and a penalty term. The former leads to zero divergence within each cell, while the latter controls the divergence error across cell interfaces. Our analysis reveals that a coupling of them is important for positivity preservation, as they exactly contribute the discrete divergence-terms absent in standard DG schemes but crucial for ensuring the PP property. Numerical tests confirm the PP property and the effectiveness of proposed PP schemes. Unlike conservative MHD system, the exact smooth solutions of symmetrizable MHD system are proved to retain the positivity even if the divergence-free condition is not satisfied. Our analysis and findings further the understanding, at both discrete and continuous levels, of the relation between the PP property and the divergence-free constraint. △ Less

Submitted 30 July, 2018; originally announced July 2018.

Comments: 49 pages, 11 figures

Journal ref: Numerische Mathematik 142, 995-1047, 2019

Provably Positive Discontinuous Galerkin Methods for Multidimensional Ideal Magnetohydrodynamics

Authors: Kailiang Wu, Chi-Wang Shu

Abstract: The density and pressure are positive physical quantities in magnetohydrodynamics (MHD). Design of provably positivity-preserving (PP) numerical schemes for ideal compressible MHD is highly desirable, but remains a challenge especially in the multidimensional cases. In this paper, we first develop uniformly high-order discontinuous Galerkin (DG) schemes which provably preserve the positivity of de… ▽ More The density and pressure are positive physical quantities in magnetohydrodynamics (MHD). Design of provably positivity-preserving (PP) numerical schemes for ideal compressible MHD is highly desirable, but remains a challenge especially in the multidimensional cases. In this paper, we first develop uniformly high-order discontinuous Galerkin (DG) schemes which provably preserve the positivity of density and pressure for multidimensional ideal MHD. The schemes are constructed by using the locally divergence-free DG schemes for the symmetrizable ideal MHD equations as the base schemes, a PP limiter to enforce the positivity of the DG solutions, and the strong stability preserving methods for time discretization. The significant innovation is that we discover and rigorously prove the PP property of the proposed DG schemes by using a novel equivalent form of the admissible state set and very technical estimates. Several two-dimensional numerical examples further confirm the PP property, and demonstrate the accuracy, effectiveness and robustness of the proposed PP methods. △ Less

Submitted 30 June, 2018; originally announced July 2018.

Comments: Submitted for publication on Jan. 31, 2018

Journal ref: SIAM Journal on Scientific Computing, 40(5), B1302-B1329, 2018

Identifying Strong-Field Effects in Indirect Photofragmentation Reactions

Authors: Chuan-Cun Shu, Kai-Jun Yuan, Daoyi Dong, Ian R. Petersen, Andre D. Bandrauk

Abstract: Exploring molecular breakup processes induced by light-matter interactions has both fundamental and practical implications. However, it remains a challenge to elucidate the underlying reaction mechanism in the strong field regime, where the potentials of the reactant are modified dramatically. Here, we perform a theoretical analysis combined with a time-dependent wavepacket calculation to show how… ▽ More Exploring molecular breakup processes induced by light-matter interactions has both fundamental and practical implications. However, it remains a challenge to elucidate the underlying reaction mechanism in the strong field regime, where the potentials of the reactant are modified dramatically. Here, we perform a theoretical analysis combined with a time-dependent wavepacket calculation to show how a strong ultrafast laser field affects the photofragment products. As an example, we examine the photochemical reaction of breaking up the molecule NaI into the neutral atoms Na and I, which due to inherent nonadiabatic couplings is indirectly formed in a stepwise fashion via the reaction intermediate NaI. By analyzing the angular dependencies of fragment distributions, we are able to identify the reaction intermediate NaI from the weak to the strong field-induced nonadiabatic regimes. Furthermore, the energy levels of NaI can be extracted from the quantum interference patterns of the transient photofragment momentum distribution. △ Less

Submitted 8 June, 2018; originally announced June 2018.

Comments: 17 pages, 5 figures

Journal ref: Journal of Physical Chemistry Letters, vol 8, p1-6, 2017

A third-order discrete unified gas kinetic scheme for continuum and rarefied flows: low-speed isothermal case

Authors: Chen Wu, Chang Shu, Baochang Shi, Zhen Chen

Abstract: An efficient third-order discrete unified gas kinetic scheme (DUGKS) with efficiency is presented in this work for simulating continuum and rarefied flows. By employing two-stage time-stepping scheme and the high-order DUGKS flux reconstruction strategy, third-order of accuracy in both time and space can be achieved in the present method. It is also analytically proven that the second-order DUGKS… ▽ More An efficient third-order discrete unified gas kinetic scheme (DUGKS) with efficiency is presented in this work for simulating continuum and rarefied flows. By employing two-stage time-stepping scheme and the high-order DUGKS flux reconstruction strategy, third-order of accuracy in both time and space can be achieved in the present method. It is also analytically proven that the second-order DUGKS is a special case of the present method. Compared with the high-order lattice Boltzmann equation {LBE} based methods, the present method is capable to deal with the rarefied flows by adopting the Newton-Cotes quadrature to approximate the integrals of moments. Instead of being constrained by the second-order (or lower-order) of accuracy in time splitting scheme as in the conventional high-order Runge-Kutta (RK) based kinetic methods, the present method solves the original BE, which overcomes the limitation in time accuracy. Typical benchmark tests are carried out for comprehensive evaluation of the present method. It is observed in the tests that the present method is advantageous over the original DUGKS in accuracy and capturing delicate flow structures. Moreover, the efficiency of the present third-order method is also shown in simulating rarefied flows. △ Less

Submitted 16 September, 2017; originally announced September 2017.

Journal ref: Phys. Rev. E 97, 023306 (2018)

Achieving robust and high-fidelity quantum control via spectral phase optimization

Authors: Yu Guo, Daoyi Dong, Chuan-Cun Shu

Abstract: Achieving high-fidelity control of quantum systems is of fundamental importance in physics, chemistry and quantum information sciences. However, the successful implementation of a high-fidelity quantum control scheme also requires robustness against control field fluctuations. Here, we demonstrate a robust optimization method for control of quantum systems by optimizing the spectral phase of an ul… ▽ More Achieving high-fidelity control of quantum systems is of fundamental importance in physics, chemistry and quantum information sciences. However, the successful implementation of a high-fidelity quantum control scheme also requires robustness against control field fluctuations. Here, we demonstrate a robust optimization method for control of quantum systems by optimizing the spectral phase of an ultrafast laser pulse, which is accomplished in the framework of frequency domain quantum optimal control theory. By incorporating a filtering function of frequency into the optimization algorithm, our numerical simulations in an abstract two-level quantum system as well as in a three-level atomic rubidium show that the optimization procedure can be enforced to search optimal solutions while achieving remarkable robustness against the control field fluctuations, providing an efficient approach to optimize the spectral phase of the ultrafast laser pulse to achieve a desired final quantum state of the system. △ Less

Submitted 2 September, 2017; originally announced September 2017.

Comments: 17 pages, 8 figures

Journal ref: Phys. Chem. Chem. Phys., 2018, 20, 9498-9506

Atomically thin mirrors made of monolayer semiconductors

Authors: Giovanni Scuri, You Zhou, Alexander A. High, Dominik S. Wild, Chi Shu, Kristiaan De Greve, Luis A. Jauregui, Takashi Taniguchi, Kenji Watanabe, Philip Kim, Mikhail D. Lukin, Hongkun Park

Abstract: Transition metal dichalcogenide monolayers are promising candidates for exploring new electronic and optical phenomena and for realizing atomically thin optoelectronic devices. They host tightly bound electron-hole pairs (excitons) that can be efficiently excited by resonant light fields. Here, we demonstrate that a single monolayer of molybdenum diselenide (MoSe2) can dramatically modify light tr… ▽ More Transition metal dichalcogenide monolayers are promising candidates for exploring new electronic and optical phenomena and for realizing atomically thin optoelectronic devices. They host tightly bound electron-hole pairs (excitons) that can be efficiently excited by resonant light fields. Here, we demonstrate that a single monolayer of molybdenum diselenide (MoSe2) can dramatically modify light transmission near the excitonic resonance, acting as an electrically switchable mirror that reflects up to 85% of incident light at cryogenic temperatures. This high reflectance is a direct consequence of the excellent coherence properties of excitons in this atomically thin semiconductor, encapsulated by hexagonal boron nitride. Furthermore, we show that the MoSe2 monolayer exhibits power- and wavelength-dependent nonlinearities that stem from exciton-based lattice heating in the case of continuous-wave excitation and exciton-exciton interactions when fast, pulsed laser excitation is used. These observations open up new possibilities for studying quantum nonlinear optical phenomena and topological photonics, and for miniaturizing optical devices. △ Less

Submitted 19 May, 2017; originally announced May 2017.

Comments: 4 figures

Journal ref: Phys. Rev. Lett. 120, 037402 (2018)

Complete elimination of nonlinear light-matter interactions with broadband ultrafast laser pulses

Authors: Chuan-Cun Shu, Daoyi Dong, Ian R. Petersen, Niels E. Henriksen

Abstract: The absorption of a single photon that excites a quantum system from a low to a high energy level is an elementary process of light-matter interaction, and a route towards realizing pure single-photon absorption has both fundamental and practical implications in quantum technology. Due to nonlinear optical effects, however, the probability of pure single-photon absorption is usually very low, whic… ▽ More The absorption of a single photon that excites a quantum system from a low to a high energy level is an elementary process of light-matter interaction, and a route towards realizing pure single-photon absorption has both fundamental and practical implications in quantum technology. Due to nonlinear optical effects, however, the probability of pure single-photon absorption is usually very low, which is particularly pertinent in the case of strong ultrafast laser pulses with broad bandwidth. Here we demonstrate theoretically a counterintuitive coherent single-photon absorption scheme by eliminating nonlinear interactions of ultrafast laser pulses with quantum systems. That is, a completely linear response of the system with respect to the spectral energy density of the incident light at the transition frequency can be obtained for all transition probabilities between 0 and 100% in a multi-level quantum systems. To that end, a new multi-objective optimization algorithm is developed to find an optimal spectral phase of an ultrafast laser pulse, which is capable of eliminating all possible nonlinear optical responses while maximizing the probability of single-photon absorption between quantum states. This work not only deepens our understanding of light-matter interactions, but also offers a new way to study photophysical and photochemical processes in the "absence" of nonlinear optical effects. △ Less

Submitted 5 May, 2017; originally announced May 2017.

Comments: 11 pages, 5 figures

Journal ref: Physical Review A 95, 033809 (2017)

Positivity-preserving flux limiters for high-order conservative schemes

Authors: X. Y. Hu, N. A. Adams, C. -W. Shu

Abstract: In this work a simple method to enforce the positivity-preserving property for general high-order conservative schemes is proposed. The method keeps the original scheme unchanged and detects critical numerical fluxes which may lead to negative density and pressure, and then imposes a cut-off flux limiter to satisfy a sufficient condition for preserving positivity. Though an extra time-step size co… ▽ More In this work a simple method to enforce the positivity-preserving property for general high-order conservative schemes is proposed. The method keeps the original scheme unchanged and detects critical numerical fluxes which may lead to negative density and pressure, and then imposes a cut-off flux limiter to satisfy a sufficient condition for preserving positivity. Though an extra time-step size condition is required to maintain the formal order of accuracy, it is less restrictive than those in previous works. A number of numerical examples suggest that this method, when applied on an essentially non-oscillatory base scheme, can be used to prevent positivity failure when the flow involves vacuum or near vacuum and very strong discontinuities. △ Less

Submitted 17 September, 2012; v1 submitted 7 March, 2012; originally announced March 2012.

Comments: 27 pages, 7 figures

Bilinear effect in complex systems

Authors: Lui Lam, David C. Bellavia, Xiao-Pu Han, Chih-Hui Alston Liu, Chang-Qing Shu, Zhengjin Wei, Tao Zhou, Jichen Zhu

Abstract: The distribution of the lifetime of Chinese dynasties (as well as that of the British Isles and Japan) in a linear Zipf plot is found to consist of two straight lines intersecting at a transition point. This two-section piecewise-linear distribution is different from the power law or the stretched exponent distribution, and is called the Bilinear Effect for short. With assumptions mimicking the or… ▽ More The distribution of the lifetime of Chinese dynasties (as well as that of the British Isles and Japan) in a linear Zipf plot is found to consist of two straight lines intersecting at a transition point. This two-section piecewise-linear distribution is different from the power law or the stretched exponent distribution, and is called the Bilinear Effect for short. With assumptions mimicking the organization of ancient Chinese regimes, a 3-layer network model is constructed. Numerical results of this model show the bilinear effect, providing a plausible explanation of the historical data. Bilinear effect in two other social systems is presented, indicating that such a piecewise-linear effect is widespread in social systems. △ Less

Submitted 14 September, 2010; v1 submitted 27 March, 2010; originally announced March 2010.

Comments: 5 pages, 5 figures

Journal ref: EPL, 91 (2010) 68004

Investigation of Gd3N@C2n (40 < n < 44) family by Raman and inelastic electron tunneling spectroscopy

Authors: Brian G. Burke, Jack Chan, Keith A. Williams, Jiechao Ge, Chunying Shu, Wujun Fu, Harry C. Dorn, James G. Kushmerick, Alexander A. Puretzky, David B. Geohegan

Abstract: The structure and vibrational spectrum of Gd3N@C80 is studied through Raman and inelastic electron tunneling spectroscopy (IETS) as well as density functional theory (DFT) and universal force field (UFF) calculations. Hindered rotations, shown by both theory and experiment, indicate the formation of a Gd3N-C80 bond which reduces the ideal icosahedral symmetry of the C80 cage. The vibrational mod… ▽ More The structure and vibrational spectrum of Gd3N@C80 is studied through Raman and inelastic electron tunneling spectroscopy (IETS) as well as density functional theory (DFT) and universal force field (UFF) calculations. Hindered rotations, shown by both theory and experiment, indicate the formation of a Gd3N-C80 bond which reduces the ideal icosahedral symmetry of the C80 cage. The vibrational modes involving the movement of the encapsulated species are a fingerprint of the interaction between the fullerene cage and the core complex. We present Raman data for the Gd3N@C2n (40 < n < 44) family as well as Y3N@C80, Lu3N@C80, and Y3N@C88 for comparison. Conductance measurements have been performed on Gd3N@C80 and reveal a Kondo effect similar to that observed in C60. △ Less

Submitted 27 October, 2009; originally announced October 2009.

Comments: 7 pages, 8 figures, paper

Low-Redshift Cosmic Baryon Fluid on Large Scales and She-Leveque Universal Scaling

Authors: Ping He, Jiren Liu, Long-Long Feng, Chi-Wang Shu, Li-Zhi Fang

Abstract: We investigate the statistical properties of cosmic baryon fluid in the nonlinear regime, which is crucial for understanding the large-scale structure formation of the universe. With the hydrodynamic simulation sample of the Universe in the cold dark matter model with a cosmological constant, we show that the intermittency of the velocity field of cosmic baryon fluid at redshift z=0 in the scale… ▽ More We investigate the statistical properties of cosmic baryon fluid in the nonlinear regime, which is crucial for understanding the large-scale structure formation of the universe. With the hydrodynamic simulation sample of the Universe in the cold dark matter model with a cosmological constant, we show that the intermittency of the velocity field of cosmic baryon fluid at redshift z=0 in the scale range from the Jeans length to about 16 Mpc/h can be extremely well described by She-Leveque's universal scaling formula. The baryon fluid also possesses the following features: (1) for volume weight statistics, the dissipative structures are dominated by sheets, and (2) the relation between the intensities of fluctuations is hierarchical. These results imply that the evolution of highly evolved cosmic baryon fluid is similar to a fully developed turbulence. △ Less

Submitted 9 February, 2006; v1 submitted 10 January, 2006; originally announced January 2006.

Comments: 4 pages, 3 figures

Journal ref: Phys.Rev.Lett. 96 (2006) 051302