Source author record

Fred Cooper

Fred Cooper appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

49works
21topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

49 published item(s)

preprint2022arXiv

Stability of exact solutions of the $(2+1)$-dimensional nonlinear Schrödinger equation with arbitrary nonlinearity parameter $κ$

In this work, we consider the nonlinear Schrödinger equation (NLSE) in $2+1$ dimensions with arbitrary nonlinearity exponent $κ$ in the presence of an external confining potential. Exact solutions to the system are constructed, and their stability over their "mass" (i.e., the $L^2$ norm) and the parameter $κ$ is explored. We observe both theoretically and numerically that the presence of the confining potential leads to wider domains of stability over the parameter space compared to the unconfined case. Our analysis suggests the existence of a stable regime of solutions for all $κ$ as long as their mass is less than a critical value $M^{\ast}(κ)$. Furthermore, we find that there are two different critical masses, one corresponding to width perturbations and the other one to translational perturbations. The results of Derrick's theorem are also obtained by studying the small amplitude regime of a four-parameter collective coordinate (4CC) approximation. A numerical stability analysis of the NLSE shows that the instability curve $M^{\ast}(κ)$ vs. $κ$ lies below the two curves found by Derrick's theorem and the 4CC approximation. In the absence of the external potential, $κ=1$ demarcates the separation between the blowup regime and the stable regime. In this 4CC approximation, for $κ<1$, when the mass is above the critical mass for the translational instability, quite complicated motions of the collective coordinates are possible. Energy conservation prevents the blowup of the solution as well as confines the center of the solution to a finite spatial domain. We call this regime the "frustrated" blowup regime and give some illustrations. In an appendix, we show how to extend these results to arbitrary initial ground state solution data and arbitrary spatial dimension $d$.

preprint2020arXiv

Neutrino condensation from a New Higgs Interaction

We study the consequences of having a new interaction between neutrinos and a Higgs scalar. We find that there are two possible attractive channels in the resulting effective 4-fermi theory which lead to a neutrino condensation in cosmic neutrinos and the creation of a neutrino superfluid at low temperatures and finite density. We find that at the minimum of the effective potential V the condensates are mostly made up of pairs of left-left+ right-right composites, with a slight admixture of left-right + right-left composites.

preprint2019arXiv

Parametrically driven nonlinear Dirac equation with arbitrary nonlinearity

The damped and parametrically driven nonlinear Dirac equation with arbitrary nonlinearity parameter $κ$ is analyzed, when the external force is periodic in space and given by $f(x) =r\cos(K x)$, both numerically and in a variational approximation using five collective coordinates (time dependent shape parameters of the wave function). Our variational approximation satisfies exactly the low-order moment equations. Because of competition between the spatial period of the external force $λ=2 π/K$, and the soliton width $l_s$, which is a function of the nonlinearity $κ$ as well as the initial frequency $ω_0$ of the solitary wave, there is a transition (at fixed $ω_0$) from trapped to unbound behavior of the soliton, which depends on the parameters $r$ and $K$ of the external force and the nonlinearity parameter $κ$. We previously studied this phenomena when $κ=1$ (2019 J. Phys. A: Math. Theor. {\bf 52} 285201) where we showed that for $λ\gg l_s$ the soliton oscillates in an effective potential, while for $λ\ll l_s$ it moves uniformly as a free particle. In this paper we focus on the $κ$ dependence of the transition from oscillatory to particle behavior and explicitly compare the curves of the transition regime found in the collective coordinate approximation as a function of $r$ and $K$ when $κ=1/2,1,2$ at fixed value of the frequency $ω_0$. Since the solitary wave gets narrower for fixed $ω_0$ as a function of $κ$, we expect and indeed find that the regime where the solitary wave is trapped is extended as we increase $κ$.

preprint2016arXiv

Stability of exact solutions of the nonlinear Schroedinger equation in an external potential having supersymmetry and parity-time symmetry

We discuss the stability properties of the solutions of the general nonlinear Schroedinger equation (NLSE) in 1+1 dimensions in an external potential derivable from a parity-time (PT) symmetric superpotential $W(x)$ that we considered earlier [Kevrekedis et al Phys. Rev. E 92, 042901 (2015)]. In particular we consider the nonlinear partial differential equation $\{ i \partial_t + \partial_x^2 - V^{-}(x) +| ψ(x,t) |^{2κ} \} \, ψ(x,t) = 0$, for arbitrary nonlinearity parameter $κ$. We study the bound state solutions when $V^{-}(x) = (1/4- b^2)$ sech$^2(x)$, which can be derived from two different superpotentials $W(x)$, one of which is complex and $PT$ symmetric. Using Derrick's theorem, as well as a time dependent variational approximation, we derive exact analytic results for the domain of stability of the trapped solution as a function of the depth $b^2$ of the external potential. We compare the regime of stability found from these analytic approaches with a numerical linear stability analysis using a variant of the Vakhitov-Kolokolov (V-K) stability criterion. The numerical results of applying the V-K condition give the same answer for the domain of stability as the analytic result obtained from applying Derrick's theorem. Our main result is that for $κ>2$ a new regime of stability for the exact solutions appears as long as $b > b_{crit}$, where $b_{crit}$ is a function of the nonlinearity parameter $κ$. In the absence of the potential the related solitary wave solutions of the NLSE are unstable for $κ>2$.

preprint2016arXiv

Variational Approach to studying solitary waves in the nonlinear Schrodinger equation with Complex Potentials

We discuss the behavior of solitary wave solutions of the nonlinear Schr{ö}dinger equation (NLSE) as they interact with complex potentials, using a four parameter variational approximation based on a dissipation functional formulation of the dynamics. We concentrate on spatially periodic potentials with the periods of the real and imaginary part being either the same or different. Our results for the time evolution of the collective coordinates of our variational ansatz are in good agreement with direct numerical simulation of the NLSE. We compare our method with a collective coordinate approach of Kominis and give examples where the two methods give qualitatively different answers. In our variational approach, we are able to give analytic results for the small oscillation frequency of the solitary wave oscillating parameters which agree with the numerical solution of the collective coordinate equations. We also verify that instabilities set in when the slope of $dp(t)/dv(t)$ becomes negative when plotted parametrically as a function of time, where $p(t)$ is the momentum of the solitary wave and $v(t)$ the velocity.

preprint2015arXiv

Auxiliary Field Loop Expansion of the Effective Action for Stochastic Partial Differential Equations

We present an alternative to the perturbative diagrammatic approach for studying stochastic dynamics. Our approach is based on an auxiliary field loop expansion for the path integral representation for the generating functional of the noise induced correlation functions. We derive two different effective actions, one based on the Onsager-Machlup (OM) approach, and the other on the Martin-Siggia-Rose (MSR) response function approach. In particular we determine the leading order approximation for the effective action and effective potential for arbitrary spatial dimensions for several simple systems. These include the Kardar-Parisi-Zhang (KPZ) equation, the chemical reaction annihilation and diffusion process $A+A \rightarrow 0$, and the Ginzburg-Landau (GL) model for spin relaxation. We show how to obtain the effective potential of the OM approach from the effective potential in the MSR approach. For the KPZ equation we find that our approximation, which is non-perturbative and obeys broken symmetry Ward identities, does not lead to the appearance of a fluctuation induced symmetry breakdown. This contradicts the results of earlier studies. We also obtain some of the renormalization group flows directly from the effective potential and compare our results with exact and perturbative results.

preprint2015arXiv

Nonlinear Dirac equation solitary waves in the presence of external driving forces

We consider the nonlinear Dirac (NLD) equation in 1+1 dimension with scalar-scalar self-interaction in the presence of external forces as well as damping of the form $ f(x,t) - i μγ^0 Ψ$, where both $f$ and $Ψ$ are two-component spinors. We develop an approximate variational approach using collective coordinates (CC) for studying the time dependent response of the solitary waves to these external forces. This approach predicts intrinsic oscillations of the solitary waves, i.e. the amplitude, width and phase all oscillate with the same frequency. The translational motion is also affected, because the soliton position oscillates around a mean trajectory. We then compare the results of the variational approximation with numerical simulations of the NLD equation, and find a good agreement, if we take into account a certain linear excitation with specific wavenumber that is excited together with the intrinsic oscillations such that the momentum in a transformed NLD equation is conserved. We also solve explicitly the CC equations of the variational approximation in the non-relativistic regime for a homogeneous external force and obtain excellent agreement with the numerical solution of the CC equations.

preprint2015arXiv

Solitary waves of a PT-symmetric Nonlinear Dirac equation

In the present work, we consider a prototypical example of a PT-symmetric Dirac model. We discuss the underlying linear limit of the model and identify the threshold of the PT-phase transition in an analytical form. We then focus on the examination of the nonlinear model. We consider the continuation in the PT-symmetric model of the solutions of the corresponding Hamiltonian model and find that the solutions can be continued robustly as stable ones all the way up to the PT-transition threshold. In the latter, they degenerate into linear waves. We also examine the dynamics of the model. Given the stability of the waveforms in the PT-exact phase we consider them as initial conditions for parameters outside of that phase. We find that both oscillatory dynamics and exponential growth may arise, depending on the size of the corresponding "quench". The former can be characterized by an interesting form of bi-frequency solutions that have been predicted on the basis of the SU(1,1) symmetry. Finally, we explore some special, analytically tractable, but not PT-symmetric solutions in the massless limit of the model.

preprint2014arXiv

Auxiliary Field Loop Expansion for the Effective Action for Stochastic Partial Differential Equations I

Using a path integral formulation for correlation functions of stochastic partial differential equations based on the Onsager-Machlup approach, we show how, by introducing a composite auxiliary field one can generate an auxiliary field loop expansion for the correlation functions which is similar to the one used in the $1/N$ expansion for an $O(N)$ scalar quantum field theory. We apply this formalism to the Kardar Parisi Zhang (KPZ) equation, and introduce the composite field $σ= \fracλ{2} \nabla ϕ\cdot \nabla ϕ$ by inserting a representation of the unit operator into the path integral which enforces this constraint. In leading order we obtain a self-consistent mean field approximation for the effective action similar to that used for the Bardeen-Cooper-Schrieffer (BCS) and Bose-Einstein Condensate (BEC) theories of dilute Fermi and Bose gases. This approximation, though related to a self-consistent Gaussian approximation, preserves all symmetries and broken symmetries. We derive the leading order in the auxiliary field (LOAF) effective potential and compare our results to the one loop in the fluctuation strength ${\cal A}$ approximation. We find, contrary to what is found in the one loop and self-consistent Gaussian approximation schemes that in the LOAF approximation there is no fluctuation induced symmetry breaking as a function of the coupling constant in any dimension $d$.

preprint2014arXiv

Auxiliary Field Loop expansion for the Effective Action for Stochastic Partial Differential equations II

We extend our discussion of effective actions for stochastic partial differential equations to systems that give rise to a Martin-Siggia-Rose (MSR) type of action. This type of action naturally arises when one uses the many-body formalism of Doi and Peliti to describe reaction-diffusion models which undergo transitions into the absorbing state and which are described by a Master equation. These models include predator prey models, and directed percolation models as well as chemical kinetic models. For classical dynamical systems with external noise it is always possible to construct an MSR action. Using a path integral representation for the generator of the correlation functions, we show how, by introducing a composite auxiliary field, one can generate an auxiliary field loop expansion for the effective action for both types of systems. As a specific example of the Doi-Peliti formalism we determine the effective action for the chemical reaction annihilation and diffusion process $A+A \rightarrow 0$. For the external noise problem we evaluate the effective action for the Cole-Hopf form of the Kardar-Parisi Zhang (KPZ) equation as well as for the Ginzburg Landau model of spin relaxation. We determine for arbitrary spatial dimension $d$, the renormalized effective potential in leading order in the auxiliary field loop expansion (LOAF) and also determine the renormalization group equation for the running of the reaction rate (coupling constant) for arbitrary $d$. We compare our results with known perturbative and non-perturbative results for the renormalization group equations.

preprint2014arXiv

Effects of intrinsic noise on a cubic autocatalytic reaction diffusion system

Starting from our recent chemical master equation derivation of the model of an autocatalytic reaction-diffusion chemical system with reactions $U+2V {\stackrel {λ_0}{\rightarrow}}~ 3 V;$ and $V {\stackrel μ{\rightarrow}}~P$, $U {\stackrel ν{\rightarrow}}~ Q$, we determine the effects of intrinsic noise on the momentum-space behavior of its kinetic parameters and chemical concentrations. We demonstrate that the intrinsic noise induces $n \rightarrow n$ molecular interaction processes with $n \geq 4$, where $n$ is the number of molecules participating of type $U$ or $V$. The momentum dependences of the reaction rates are driven by the fact that the autocatalytic reaction (inelastic scattering) is renormalized through the existence of an arbitrary number of intermediate elastic scatterings, which can also be interpreted as the creation and subsequent decay of a three body composite state $σ= ϕ_u ϕ_v^2$, where $ϕ_i$ corresponds to the fields representing the densities of $U$ and $V$. Finally, we discuss the difference between representing $σ$ as a composite or an elementary particle (molecule) with its own kinetic parameters. In one dimension we find that while they show markedly different behavior in the short spatio-temporal scale, high momentum (UV) limit, they are formally equivalent in the large spatio-temporal scale, low momentum (IR) regime. On the other hand in two dimensions and greater, due to the effects of fluctuations, there is no way to experimentally distinguish between a fundamental and composite $σ$. Thus in this regime $σ$ behave as an entity unto itself suggesting that it can be effectively treated as an independent chemical species.

preprint2014arXiv

Stability of solitary waves in the nonlinear Dirac equation with arbitrary nonlinearity

We consider the nonlinear Dirac equation in 1+1 dimension with scalar-scalar self interaction $ \frac{g^2}{κ+1} ({\bar Ψ} Ψ)^{κ+1}$ and with mass $m$. Using the exact analytic form for rest frame solitary waves of the form $Ψ(x,t) = ψ(x) e^{-i ωt}$ for arbitrary $ κ$, we discuss the validity of various approaches to understanding stability that were successful for the nonlinear Schrödinger equation. In particular we study the validity of a version of Derrick's theorem, the criterion of Bogolubsky as well as the Vakhitov-Kolokolov criterion, and find that these criteria yield inconsistent results. Therefore, we study the stability by numerical simulations using a recently developed 4th-order operator splitting integration method. For different ranges of $κ$ we map out the stability regimes in $ω$. We find that all stable nonlinear Dirac solitary waves have a one-hump profile, but not all one-hump waves are stable, while all waves with two humps are unstable. We also find that the time $t_c$, it takes for the instability to set in, is an exponentially increasing function of $ω$ and $t_c$ decreases monotonically with increasing $κ$.

preprint2013arXiv

Composite Bound States and Broken U(1) symmetry in the Chemical Master Equation derivation of the Gray-Scott Model

We give a first principles derivation of the stochastic partial differential equations that describe the chemical reactions of the Gray-Scott model (GS): $U+2V {\stackrel λ{\rightarrow}} 3 V;$ and $V {\stackrel μ{\rightarrow}} P$, $U {\stackrel ν{\rightarrow}} Q$, with a constant feed rate for $U$. We find that the conservation of probability ensured by the chemical master equation leads to a modification of the usual differential equations for the GS model which now involves two composite fields and also intrinsic noise terms. One of the composites is $ψ_1 = ϕ_v^2$, where $ < ϕ_v >_η = v$ is the concentration of the species $V$ and the averaging is over the internal noise $η_{u,v,ψ_1}$. The second composite field is the product of three fields $ χ= λϕ_u ϕ_v^2$ and requires a noise source to ensure probability conservation. A third composite $ψ_2 = ϕ_{u} ϕ_{v}$ can be also be identified from the noise-induced reactions. The Hamiltonian that governs the time evolution of the many-body wave function, associated with the master equation, has a broken U(1) symmetry related to particle number conservation. By expanding around the (broken symmetry) zero energy solution of the Hamiltonian (by performing a Doi shift) one obtains from our path integral formulation the usual reaction diffusion equation, at the classical level. The Langevin equations that are derived from the chemical master equation have multiplicative noise sources for the density fields $ϕ_u, ϕ_v, χ$ that induce higher order processes such as $n \rightarrow n$ scattering for $n > 3$. The amplitude of the noise acting on $ ϕ_v$ is itself stochastic in nature.

preprint2013arXiv

Internal composite bound states in deterministic reaction diffusion models

By identifying potential composite states that occur in the Sel'kov-Gray-Scott (GS) model, we show that it can be considered as an effective theory at large spatio-temporal scales, arising from a more \textit{fundamental} theory (which treats these composite states as fundamental chemical species obeying the diffusion equation) relevant at shorter spatio-temporal scales. When simulations in the latter model are performed as a function of a parameter $M = λ^{-1}$, the generated spatial patterns evolve at late times into those of the GS model at large $M$, implying that the composites follow their own unique dynamics at short scales. This separation of scales is an example of \textit{dynamical} decoupling in reaction diffusion systems.

preprint2013arXiv

Leading-Order Auxiliary Field Theory of the Bose-Hubbard Model

We discuss the phase diagram of the Bose-Hubbard (BH) model in the leading-order auxiliary field (LOAF) theory. LOAF is a conserving non-perturbative approximation that treats on equal footing the normal and anomalous density condensates. The mean-field solutions in LOAF correspond to first-order and second-order phase transition solutions with two critical temperatures corresponding to a vanishing Bose-Einstein condensate, $T_c$, and a vanishing diatom condensate, $T^\star$. The \emph{second-order} phase transition solution predicts the correct order of the transition in continuum Bose gases. For either solution, the superfluid state is tied to the presence of the diatom condensate related to the anomalous density in the system. In ultracold Bose atomic gases confined on a three-dimensional lattice, the critical temperature $T_c$ exhibits a quantum phase transition, where $T_c$ goes to zero at a finite coupling. The BH phase diagram in LOAF features a line of first-order transitions ending in a critical point beyond which the transition is second order while approaching the quantum phase transition. We identify a region where a diatom condensate is expected for temperatures higher than $T_c$ and less than $T_0$, the critical temperature of the non-interacting system. The LOAF phase diagram for the BH model compares qualitatively well with existing experimental data and results of \emph{ab initio} Monte Carlo simulations.

preprint2013arXiv

Nonlinear Dirac equation solitary waves in external fields

We consider the nonlinear Dirac equations (NLDE's) in 1+1 dimension with scalar-scalar self interaction $\frac{g^2}{κ+1} ({\bPsi} Ψ)^{κ+1}$ in the presence of various external electromagnetic fields. Starting from the exact solutions for the unforced problem we study the behavior of solitary wave solutions to the NLDE in the presence of a wide variety of fields in a variational approximation depending on collective coordinates which allows the position, width and phase of these waves to vary in time. We find that in this approximation the position $q(t)$ of the center of the solitary wave obeys the usual behavior of a relativistic point particle in an external field. For time independent external fields we find that the energy of the solitary wave is conserved but not the momentum which becomes a function of time. We postulate that similar to the nonlinear Schr{ö}dinger equation (NLSE) that a sufficient dynamical condition for instability to arise is that $ dP(t)/d \dq(t) < 0$. Here $P(t)$ is the momentum of the solitary wave, and $\dq$ is the velocity of the center of the wave in the collective coordinate approximation. We found for our choices of external potentials we always have $ dP(t)/d \dq(t) > 0$ so when instabilities do occur they are due to a different source. We investigate the accuracy of our variational approximation using numerical simulations of the NLDE and find that when the forcing term is small and we are in a regime where the solitary wave is stable, that the behavior of the solutions of the collective coordinate equations agrees very well with the numerical simulations.

preprint2013arXiv

Quench Dynamics and Emergence of Phase Separation in Two-Component Atomic Bose Gases at Zero Temperature and above the BEC Critical Temperature

We study the dynamics of two-component atomic Bose gases initially in a mixture encountering a sudden quench of the inter-species interactions. The dynamics above the critical temperature $T_c$ is studied using a leading order large-N approximation that predicted a phase transition from mixing to phase separation as a function of the inter-species coupling. Here we explore the dynamics of this phase transition following a quench and compare our results to those found at zero temperature using the time-dependent Gross-Pitaevskii equations which ignore quantum and thermal fluctuations. In the regime above $T_c$ where no condensate is present, however, the time evolution of the densities following the quench exhibits features similar to that found at zero temperature where only the condensates contribute to the densities. When the inter-species interaction jumps above the critical value, we observe dynamical transitions from a homogeneous mixture to a phase-separated structure for both cases. Our simulations suggest that at temperatures above $T_c$ where no condensate is present this dynamical transition should still be observable in experiments.

preprint2012arXiv

Composite-Field Goldstone States and Higgs Mechanism in Dilute Bose Gases

We show that a composite-field (diatom) Goldstone state is expected in a dilute Bose gas for temperatures between the Bose gas critical temperature where the atom Bose-Einstein condensate appears and the temperature where superfluidity sets in. The presence of superfluidity is tied to the existence of a U(1) charge-two diatom condensate in the system. By promoting the global U(1) symmetry of the theory to a gauge symmetry, we find that the mass of the gauge particle generated through the Anderson-Higgs mechanism is related to the superfluid density via the Meissner effect and the superfluid density is related to the square of the anomalous density in the Bose system.

preprint2012arXiv

Large-N approximation for single- and two-component dilute Bose gases

We discuss the mean-field theories obtained from the leading order in a large-$N$ approximation for one- and two- component dilute Bose gases. For a one-component Bose gas this approximation has the following properties: the Bose-Einstein condensation (BEC) phase transition is second order but the critical temperature $T_c$ is not shifted from the non-interacting gas value $T_0$. The spectrum of excitations in the BEC phase resembles the Bogoliubov dispersion with the usual coupling constant replaced by the running coupling constant which depends on both temperature and momentum. We then study two-component Bose gases with both inter- and intra- species interactions and focus on the stability of the mixture state above $T_c$. Our mean-field approximation predicts an instability from the mixture state to a phase-separated state when the ratio of the inter-species interaction strength to the intra-species interaction strength (assuming equal strength for both species) exceeds a critical value. At high temperature this is a structural transition and the global translational symmetry is broken. Our work complements previous studies on the instability of the mixture phase in the presence of BEC.

preprint2012arXiv

Mean-field description of pairing effects, BKT physics, and superfluidity in 2D Bose gases

We derive a mean-field description for two-dimensional (2D) interacting Bose gases at arbitrary temperatures. We find that genuine Bose-Einstein condensation with long-range coherence only survives at zero temperature. At finite temperatures, many-body pairing effects included in our mean-field theory introduce a finite amplitude for the pairing density, which results in a finite superfluid density. We incorporate Berenzinskii-Kosterlitz-Thouless (BKT) physics into our model by considering the phase fluctuations of our pairing field. This then leads to the result that the superfluid phase is only stable below the BKT temperature due to these phase fluctuations. In the weakly interacting regime at low temperature we compare our theory to previous results from perturbative calculations, renormalization group calculations as well as Monte Carlo simulations. We present a finite-temperature phase diagram of 2D Bose gases. One signature of the finite amplitude of the pairing density field is a two-peak structure in the single-particle spectral function, resembling that of the pseudogap phase in 2D attractive Fermi gases.

preprint2011arXiv

Analytical limits for cold atom Bose gases with tunable interactions

We discuss the equilibrium properties of dilute Bose gases using a non-perturbative formalism based on auxiliary fields related to the normal and anomalous densities. We show analytically that for a dilute Bose gas of weakly-interacting particles at zero temperature, the leading-order auxiliary field (LOAF) approximation leads to well-known analytical results. Close to the critical point the LOAF predictions are the same as those obtained using an effective field theory in the large-N approximation. We also report analytical approximations for the LOAF results in the unitarity limit, which compare favorably with our numerical results. LOAF predicts that the equation of state for the Bose gas in the unitarity limit is E / (p V) = 1, unlike the case of the Fermi gas when E / (p V) = 3/2.

preprint2011arXiv

Auxiliary field approach to dilute Bose gases with tunable interactions

We rewrite the Lagrangian for a dilute Bose gas in terms of auxiliary fields related to the normal and anomalous condensate densities. We derive the loop expansion of the effective action in the composite-field propagators. The lowest-order auxiliary field (LOAF) theory is a conserving mean-field approximation consistent with the Goldstone theorem without some of the difficulties plaguing approximations such as the Hartree and Popov approximations. LOAF predicts a second-order phase transition. We give a set of Feynman rules for improving results to any order in the loop expansion in terms of composite-field propagators. We compare results of the LOAF approximation with those derived using the Popov approximation. LOAF allows us to explore the critical regime for all values of the coupling constant and we determine various parameters in the unitarity limit.

preprint2011arXiv

Auxiliary field formalism for dilute fermionic atom gases with tunable interactions

We develop the auxiliary field formalism corresponding to a dilute system of spin-1/2 fermions. This theory represents the Fermi counterpart of the BEC theory developed recently by F. Cooper et al. [Phys. Rev. Lett. 105, 240402 (2010)] to describe a dilute gas of Bose particles. Assuming tunable interactions, this formalism is appropriate for the study of the crossover from the regime of Bardeen-Cooper-Schriffer (BCS) pairing to the regime of Bose-Einstein condensation (BEC) in ultracold fermionic atom gases. We show that when applied to the Fermi case at zero temperature, the leading-order auxiliary field (LOAF) approximation gives the same equations as those obtained in the standard BCS variational picture. At finite temperature, LOAF leads to the theory discussed by by Sa de Melo, Randeria, and Engelbrecht [Phys. Rev. Lett. 71, 3202(1993); Phys. Rev. B 55, 15153(1997)]. As such, LOAF provides a unified framework to study the interacting Fermi gas. The mean-field results discussed here can be systematically improved upon by calculating the one-particle irreducible (1-PI) action corrections, order by order.

preprint2011arXiv

Forced Nonlinear Schroedinger Equation with Arbitrary Nonlinearity

We consider the nonlinear Schr{ö}dinger equation (NLSE) in 1+1 dimension with scalar-scalar self interaction $\frac{g^2}{κ+1} (ψ^\star ψ)^{κ+1}$ in the presence of the external forcing terms of the form $r e^{-i(kx + θ)} -δψ$. We find new exact solutions for this problem and show that the solitary wave momentum is conserved in a moving frame where $v_k=2 k$. These new exact solutions reduce to the constant phase solutions of the unforced problem when $r \rightarrow 0.$ In particular we study the behavior of solitary wave solutions in the presence of these external forces in a variational approximation which allows the position, momentum, width and phase of these waves to vary in time. We show that the stationary solutions of the variational equations include a solution close to the exact one and we study small oscillations around all the stationary solutions. We postulate that the dynamical condition for instability is that $ dp(t)/d \dot{q} (t) < 0$, where $p(t)$ is the normalized canonical momentum $p(t) = \frac{1}{M(t)} \frac {\partial L}{\partial {\dot q}}$, and $\dot{q}(t)$ is the solitary wave velocity. Here $M(t) = \int dx ψ^\star(x,t) ψ(x,t)$. Stability is also studied using a "phase portrait" of the soliton, where its dynamics is represented by two-dimensional projections of its trajectory in the four-dimensional space of collective coordinates. The criterion for stability of a soliton is that its trajectory is a closed single curve with a positive sense of rotation around a fixed point. We investigate the accuracy of our variational approximation and these criteria using numerical simulations of the NLSE.

preprint2011arXiv

On the properties of compacton-anticompacton collisions

We study the properties of compacton-anticompacton collision processes. We compare and con- trast results for the case of compacton-anticompacton solutions of the K(l, p) Rosenau-Hyman (RH) equation for l = p = 2, with compacton-anticompacton solutions of the L(l,p) Cooper-Shepard- Sodano (CSS) equation for p = 1 and l = 3. This study is performed using a Padé discretization of the RH and CSS equations. We find a significant difference in the behavior of compacton- anticompacton scattering. For the CSS equation, the scattering can be interpreted as "annihila- tion" as the wake left behind dissolves over time. In the RH equation, the numerical evidence is that multiple shocks form after the collision which eventually lead to "blowup" of the resulting waveform.

preprint2011arXiv

The Josephson relation for the superfluid density and the connection to the Goldstone theorem in dilute Bose atomic gasses

We derive the Josephson relation for a dilute Bose gas in the framework of an auxiliary-field resummation of the theory in terms of the normal- and anomalous-density condensates. The mean-field phase diagram of this theory features two critical temperatures, T_c and $T^*, associated with the presence in the system of the Bose-Einstein condensate (BEC) and superfluid state, respectively. In this context, the Josephson relation shows that the superfluid density is related to a second order parameter, the square of the anomalous-density condensate. This is in contrast with the corresponding result in the Bose gas theory without an anomalous condensate, which predicts that the superfluid density is proportional to the BEC condensate density. Our findings are consistent with the prediction that in the temperature range between T_c and T^* a fraction of the system is in the superfluid state in the absence of the BEC condensate. This situation is similar to the case of dilute Fermi gases, where the superfluid density is proportional to the square of the gap parameter. The Josephson relation relies on the existence of zero energy and momentum excitations showing the intimate relationship between superfluidity and the Goldstone theorem.

preprint2010arXiv

Non-perturbative predictions for cold atom Bose gases with tunable interactions

We derive a theoretical description for dilute Bose gases as a loop expansion in terms of composite-field propagators by rewriting the Lagrangian in terms of auxiliary fields related to the normal and anomalous densities. We demonstrate that already in leading order this non-perturbative approach describes a large interval of coupling-constant values, satisfies Goldstone's theorem, yields a Bose-Einstein transition that is second-order, and is consistent with the critical temperature predicted in the weak-coupling limit by the next-to-leading order large-N expansion.

preprint2010arXiv

Solitary wave in the Nonlinear Dirac Equation with arbitrary nonlinearity

We consider the nonlinear Dirac equations (NLDE's) in 1+1 dimension with scalar-scalar self interaction $\frac{g^2}{k+1} ({\bar Ψ} Ψ)^{k+1}$, as well as a vector-vector self interaction $\frac{g^2}{k+1} ({\bar Ψ} γ_μΨ\bPsi γ^μΨ)^{\frac{1}{2}(k+1)}$. We find the exact analytic form for solitary waves for arbitrary $k$ and find that they are a generalization of the exact solutions for the nonlinear Schrödinger equation (NLSE) and reduce to these solutions in a well defined nonrelativistic limit. We perform the nonrelativistic reduction and find the $1/2m$ correction to the NLSE, valid when $|ω-m |\ll 2m$, where $ω$ is the frequency of the solitary wave in the rest frame. We discuss the stability and blowup of solitary waves assuming the modified NLSE is valid and find that they should be stable for $k < 2$.

preprint2010arXiv

Stability and dynamical properties of Cooper-Shepard-Sodano compactons

Extending a Pade approximant method used for studying compactons in the Rosenau-Hyman (RH) equation, we study the numerical stability of single compactons of the Cooper-Shepard-Sodano (CSS) equation and their pairwise interactions. The CSS equation has a conserved Hamiltonian which has allowed various approaches for studying analytically the nonlinear stability of the solutions. We study three different compacton solutions and find they are numerically stable. Similar to the collisions between RH compactons, the CSS compactons reemerge with same coherent shape when scattered. The time evolution of the small-amplitude ripple resulting after scattering depends on the values of the parameters $l$ and $p$ characterizing the corresponding CSS equation. The simulation of the CSS compacton scattering requires a much smaller artificial viscosity to obtain numerical stability, than in the case of RH compactons propagation.

preprint2010arXiv

Stability and dynamical properties of Rosenau-Hyman compactons using Pade approximants

We present a systematic approach for calculating higher-order derivatives of smooth functions on a uniform grid using Padé approximants. We illustrate our findings by deriving higher-order approximations using traditional second-order finite-differences formulas as our starting point. We employ these schemes to study the stability and dynamical properties of K(2,2) Rosenau-Hyman (RH) compactons including the collision of two compactons and resultant shock formation. Our approach uses a differencing scheme involving only nearest and next-to-nearest neighbors on a uniform spatial grid. The partial differential equation for the compactons involves first, second and third partial derivatives in the spatial coordinate and we concentrate on four different fourth-order methods which differ in the possibility of increasing the degree of accuracy (or not) of one of the spatial derivatives to sixth order. A method designed to reduce roundoff errors was found to be the most accurate approximation in stability studies of single solitary waves, even though all derivates are accurate only to fourth order. Simulating compacton scattering requires the addition of fourth derivatives related to artificial viscosity. For those problems the different choices lead to different amounts of "spurious" radiation and we compare the virtues of the different choices.

preprint2009arXiv

Backreaction and Particle Production in (3+1)-dimensional QED

We study the fermion pair production from a strong electric field in boost-invariant coordinates in (3+1) dimensions and exploit the cylindrical symmetry of the problem. This problem has been used previously as a toy model for populating the central-rapidity region of a heavy-ion collision (when we can replace the electric by a chromoelectric field). We derive and solve the renormalized equations for the dynamics of the mean electric field and current of the produced particles, when the field is taken to be a function only of the fluid proper time $τ= \sqrt{t^2-z^2}$. We determine the proper-time evolution of the comoving energy density and pressure of the ensuing plasma and the time evolution of suitable interpolating number operators. We find that unlike in (1+1) dimensions, the energy density closely follows the longitudinal pressure. The transverse momentum distribution of fermion pairs at large momentum is quite different and larger than that expected from the constant field result.

preprint2008arXiv

Casimir dependence of transverse distribution of pairs produced from a strong constant chromo-electric background field

The transverse distribution of gluon and quark-antiquark pairs produced from a strong constant chromo-electric field depends on two gauge invariant quantities, $C_1=E^aE^a$ and $C_2=[d_{abc}E^aE^bE^c]^2$, as shown earlier in [G.C. Nayak and P. van Nieuwenhuizen, Phys. Rev. D 71, 125001 (2005)] for gluons and in [G.C. Nayak, Phys. Rev. D 72, 125010 (2005)] for quarks. Here, we discuss the explicit dependence of the distribution on the second Casimir invariant, C_2, and show the dependence is at most a 15% effect.

preprint2008arXiv

Compactons in PT-symmetric generalized Korteweg-de Vries Equations

In an earlier paper Cooper, Shepard, and Sodano introduced a generalized KdV equation that can exhibit the kinds of compacton solitary waves that were first seen in equations studied by Rosenau and Hyman. This paper considers the PT-symmetric extensions of the equations examined by Cooper, Shepard, and Sodano. From the scaling properties of the PT-symmetric equations a general theorem relating the energy, momentum, and velocity of any solitary-wave solution of the generalized KdV equation is derived, and it is shown that the velocity of the solitons is determined by their amplitude, width, and momentum.

preprint2005arXiv

Exact Elliptic Compactons in Generalized Korteweg-DeVries Equations

We derive a general theorem relating the energy, momentum and velocity of any solitary wave solution of the generalized KdV equation which enables us to relate the amplitude, width, and momentum to the velocity of these solutions. We obtain the general condition for linear and Lyapunov stability. We then obtain a two parameter family of exact solutions to these equations which include elliptic and hyper-elliptic compacton solutions. For this general family we explicitly verify both the theorem and the stability criteria.

preprint2005arXiv

Renormalized broken-symmetry Schwinger-Dyson equations and the 2PI-1/N expansion for the O(N) model

We derive the renormalized Schwinger-Dyson equations for the one- and two-point functions in the auxiliary field formulation of $λϕ^4$ field theory to order 1/N in the 2PI-1/N expansion. We show that the renormalization of the broken-symmetry theory depends only on the counter terms of the symmetric theory with $ϕ= 0$. We find that the 2PI-1/N expansion violates the Goldstone theorem at order 1/N. In using the O(4) model as a low energy effective field theory of pions to study the time evolution of disoriented chiral condensates one has to {\em{explicitly}} break the O(4) symmetry to give the physical pions a nonzero mass. In this effective theory the {\em additional} small contribution to the pion mass due to the violation of the Goldstone theorem in the 2-PI-1/N equations should be numerically unimportant.

preprint2004arXiv

J/ΨProduction in pp Collisions at \sqrt s = 200 GeV at RHIC

We study J/ψproduction in pp collisions at RHIC within the PHENIX detector acceptance range using the color singlet and color octet mechanism which are based on pQCD and NRQCD. Here we show that the color octet mechanism reproduces the RHIC data for J/ψproduction in pp collisions with respect to the p_T distribution, the rapidity distribution and the total cross section at \sqrt s = 200 GeV. The color singlet mechanism leads to a relatively small contribution to the total cross section when compared to the octet contribution.

preprint2003arXiv

Color Octet Contribution to High p_T J/ΨProduction in pp Collisions at \sqrt s = 500 and 200 GeV at RHIC

We compute \frac{dσ}{dp_T} of the J/ψproduction in pp collisions at RHIC at \sqrt s = 500 and 200 GeV by using both the color octet and singlet models in the framework of non-relativistic QCD. The J/ψwe compute here includes the direct J/ψfrom the partonic fusion processes and the J/ψcoming from the radiative decays of χ_J's both in the color octet and singlet channel. The high p_T J/ψproduction cross section is computed within the PHENIX detector acceptance ranges: -0.35 < η< 0.35 and 1.2 < \|η| < 2.4, the central electron and forward muon arms. It is found that the color octet contribution to J/ψproduction is dominant at RHIC energy in comparison to the color singlet contributions. We compare our results with the recent preliminary data obtained by PHENIX detector for the high p_T J/ψmeasurements. While the color singlet model fails to explain the data completely the color octet model is in agreement with the single data point above 2 GeV transverse momentum. A measurement of J/ψproduction at RHIC in the next run with better statistics will allow us to determine the validity of the color octet model of J/ψproduction at RHIC energies. This is very important because it is necessary to know the exact mechanism for J/ψproduction in pp collisions at RHIC if one is to make predictions of J/ψsuppression as a signature of quark-gluon plasma. These mechanisms also play an important role in determining the polarized spin structure function of the proton at RHIC.

preprint2003arXiv

Interaction of a TeV Scale Black Hole with the Quark-Gluon Plasma at LHC

If the fundamental Planck scale is near a TeV, then parton collisions with high enough center-of-mass energy should produce black holes. The production rate for such black holes has been extensively studied for the case of a proton-proton collision at \sqrt s = 14 TeV and for a lead-lead collision at \sqrt s = 5.5 TeV at LHC. As the parton energy density is much higher at lead-lead collisions than in pp collisions at LHC, one natural question is whether the produced black holes will be able to absorb the partons formed in the lead-lead collisions and eventually `eat' the quark-gluon plasma formed at LHC. In this paper, we make a quantitative analysis of this possibility and find that since the energy density of partons formed in lead-lead collisions at LHC is about 500 GeV/fm^3, the rate of absorption for one of these black holes is much smaller than the rate of evaporation. Hence, we argue that black holes formed in such collisions will decay very quickly, and will not absorb very many nearby partons. More precisely, we show that for the black hole mass to increase via parton absorption at the LHC the typical energy density of quarks and gluons should be of the order of 10^{10} GeV/fm^3. As LHC will not be able to produce such a high energy density partonic system, the black hole will not be able to absorb a sufficient number of nearby partons before it decays. The typical life time of the black hole formed at LHC is found to be a small fraction of a fm/c.

preprint2003arXiv

Minijet Initial Conditions For Non-Equilibrium Parton Evolution at RHIC and LHC

An important ingredient for the non-equilibrium evolution of partons at RHIC and LHC is to have some physically reasonable initial conditions for the single particle phase space distribution functions for the partons. We consider several plausible parametrizations of initial conditions for the single particle distribution function f_i(x,p) and fix the parameters by matching \int f(x,p)p^μd σ_μto the invariant momentum space semi-hard parton distributions obtained using perturbative QCD (pQCD), as well as fitting low order moments of the distribution function. We consider parametrizations of f_i(x,p) with both boost invariant and boost non-invariant assumptions. We determine the initial number density, energy density and the corresponding (effective) temperature of the minijet plasma at RHIC and LHC energies. For a boost non-invariant minijet phase-space distribution function we obtain ~ 30(140) /fm^3 as the initial number density, ~ 50(520) GeV/fm^3 as the initial energy density and ~ 520(930) MeV as the corresponding initial effective temperature at RHIC(LHC).

preprint2002arXiv

Quantum dynamics of phase transitions in broken symmetry $λϕ^4$ field theory

We perform a detailed numerical investigation of the dynamics of broken symmetry $λϕ^4$ field theory in 1+1 dimensions using a Schwinger-Dyson equation truncation scheme based on ignoring vertex corrections. In an earlier paper, we called this the bare vertex approximation (BVA). We assume the initial state is described by a Gaussian density matrix peaked around some non-zero value of $<ϕ(0)>$, and characterized by a single particle Bose-Einstein distribution function at a given temperature. We compute the evolution of the system using three different approximations: Hartree, BVA and a related 2PI-1/N expansion, as a function of coupling strength and initial temperature. In the Hartree approximation, the static phase diagram shows that there is a first order phase transition for this system. As we change the initial starting temperature of the system, we find that the BVA relaxes to a new final temperature and exhibits a second order phase transition. We find that the average fields thermalize for arbitrary initial conditions in the BVA, unlike the behavior exhibited by the Hartree approximation, and we illustrate how $<ϕ(t)>$ and $<χ(t)>$ depend on the initial temperature and on the coupling constant. We find that the 2PI-1/N expansion gives dramatically different results for $<ϕ(t)>$.

preprint1997arXiv

Compacton Solutions in a Class of Generalized Fifth Order Korteweg-de Vries Equations

We study a class of generalized fifth order Korteweg-de Vries (KdV) equations which are derivable from a Lagrangian L(p,m,n,l) which has variable powers of the first and second derivatives of the field with powers given by the parameters p,m,n,l. The resulting field equation has solitary wave solutions of both the usual (non-compact) and compact variety ("compactons"). For the particular case that p=m=n+l, the solitary wave solutions have compact support and the feature that their width is independent of the amplitude. We discuss the Hamiltonian structure of these theories and find that mass, momentum, and energy are conserved. We find in general that these are not completely integrable systems. Numerical simulations show that an arbitrary compact initial wave packet whose width is wider than that of a compacton breaks up into several compactons all having the same width. The scattering of two compactons is almost elastic, with the left over wake eventually turning into compacton-anticompacton pairs. When there are two different compacton solutions for a single set of parameters the wider solution is stable, and this solution is a minimum of the Hamiltonian.

preprint1997arXiv

Order 1/N corrections to the time-dependent Hartree approximation for a system of N+1 oscillators

We solve numerically to order 1/N the time evolution of a quantum dynamical system of N oscillators of mass m coupled quadratically to a massless dynamic variable. We use Schwinger's closed time path (CTP) formalism to derive the equations. We compare two methods which differ by terms of order 1/N^2. The first method is a direct perturbation theory in 1/N using the path integral. The second solves exactly the theory defined by the effective action to order 1/N. We compare the results of both methods as a function of N. At N=1, where we expect the expansion to be quite innacurate, we compare our results to an exact numerical solution of the Schroedinger equation. In this case we find that when the two methods disagree they also diverge from the exact answer. We also find at N=1 that the 1/N corrected evolutions track the exact answer for the expectation values much longer than the mean field (N= \infty) result.

preprint1996arXiv

Nonequilibrium Dynamics of Symmetry Breaking in Lambda Phi^4 Field Theory

The time evolution of O(N) symmetric lambda Phi^4 scalar field theory is studied in the large N limit. In this limit the <Phi> mean field and two-point correlation function <Phi Phi> evolve together as a self-consistent closed Hamiltonian system, characterized by a Gaussian density matrix. The static part of the effective Hamiltonian defines the True Effective Potential U_eff for configurations far from thermal equilibrium. Numerically solving the time evolution equations for energy densities corresponding to a quench in the unstable spinodal region, we find results quite different from what might be inferred from the equilibrium free energy ``effective'' potential F. Typical time evolutions show effectively irreversible energy flow from the coherent mean fields to the quantum fluctuating modes, due to the creation of massless Goldstone bosons near threshold. The plasma frequency and collisionless damping rate of the mean fields are calculated in terms of the particle number density by a linear response analysis and compared with the numerical results. Dephasing of the fluctuations leads also to the growth of an effective entropy and the transition from quantum to classical behavior of the ensemble. In addition to casting some light on fundamental issues of nonequilibrium quantum statistical mechanics, the general framework presented in this work may be applied to a study of the dynamics of second order phase transitions in a wide variety of Landau-Ginsburg systems described by a scalar order parameter.

preprint1996arXiv

Time evolution of the chiral phase transition during a spherical expansion

We examine the non-equilibrium time evolution of the hadronic plasma produced in a relativistic heavy ion collision, assuming a spherical expansion into the vacuum. We study the $O(4)$ linear sigma model to leading order in a large-$N$ expansion. Starting at a temperature above the phase transition, the system expands and cools, finally settling into the broken symmetry vacuum state. We consider the proper time evolution of the effective pion mass, the order parameter $\langle σ\rangle$, and the particle number distribution. We examine several different initial conditions and look for instabilities (exponentially growing long wavelength modes) which can lead to the formation of disoriented chiral condensates (DCCs). We find that instabilities exist for proper times which are less than 3 fm/c. We also show that an experimental signature of domain growth is an increase in the low momentum spectrum of outgoing pions when compared to an expansion in thermal equilibrium. In comparison to particle production during a longitudinal expansion, we find that in a spherical expansion the system reaches the ``out'' regime much faster and more particles get produced. However the size of the unstable region, which is related to the domain size of DCCs, is not enhanced.

preprint1995arXiv

SPHERICALLY SYMMETRIC RANDOM WALKS I. REPRESENTATION IN TERMS OF ORTHOGONAL POLYNOMIALS

Spherically symmetric random walks in arbitrary dimension $D$ can be described in terms of Gegenbauer (ultraspherical) polynomials. For example, Legendre polynomials can be used to represent the special case of two-dimensional spherically symmetric random walks. In general, there is a connection between orthogonal polynomials and semibounded one-dimensional random walks; such a random walk can be viewed as taking place on the set of integers $n$, $n=0,~1,~2,~\ldots$, that index the polynomials. This connection allows one to express random-walk probabilities as weighted inner products of the polynomials. The correspondence between polynomials and random walks is exploited here to construct and analyze spherically symmetric random walks in $D$-dimensional space, where $D$ is {\sl not} restricted to be an integer. The weighted inner-product representation is used to calculate exact closed-form spatial and temporal moments of the probability distribution associated with the random walk. The polynomial representation of spherically symmetric random walks is also used to calculate the two-point Green's function for a rotationally symmetric free scalar quantum field theory.

preprint1994arXiv

Semiquantum Chaos and the Large N Expansion

We consider the dynamical system consisting of a quantum degree of freedom $A$ interacting with $N$ quantum oscillators described by the Lagrangian \bq L = {1\over 2}\dot{A}^2 + \sum_{i=1}^{N} \left\{{1\over 2}\dot{x}_i^2 - {1\over 2}( m^2 + e^2 A^2)x_i^2 \right\}. \eq In the limit $N \rightarrow \infty$, with $e^2 N$ fixed, the quantum fluctuations in $A$ are of order $1/N$. In this limit, the $x$ oscillators behave as harmonic oscillators with a time dependent mass determined by the solution of a semiclassical equation for the expectation value $\VEV{A(t)}$. This system can be described, when $\VEV{x(t)}= 0$, by a classical Hamiltonian for the variables $G(t) = \VEV{x^2(t)}$, $\dot{G}(t)$, $A_c(t) = \VEV{A(t)}$, and $\dot{A_c}(t)$. The dynamics of this latter system turns out to be chaotic. We propose to study the nature of this large-$N$ limit by considering both the exact quantum system as well as by studying an expansion in powers of $1/N$ for the equations of motion using the closed time path formalism of quantum dynamics.

preprint1994arXiv

Supersymmetry and Quantum Mechanics

In the past ten years, the ideas of supersymmetry have been profitably applied to many nonrelativistic quantum mechanical problems. In particular, there is now a much deeper understanding of why certain potentials are analytically solvable and an array of powerful new approximation methods for handling potentials which are not exactly solvable. In this report, we review the theoretical formulation of supersymmetric quantum mechanics and discuss many applications. Exactly solvable potentials can be understood in terms of a few basic ideas which include supersymmetric partner potentials, shape invariance and operator transformations. Familiar solvable potentials all have the property of shape invariance. We describe new exactly solvable shape invariant potentials which include the recently discovered self-similar potentials as a special case. The connection between inverse scattering, isospectral potentials and supersymmetric quantum mechanics is discussed and multi-soliton solutions of the KdV equation are constructed. Approximation methods are also discussed within the framework of supersymmetric quantum mechanics and in particular it is shown that a supersymmetry inspired WKB approximation is exact for a class of shape invariant potentials. Supersymmetry ideas give particularly nice results for the tunneling rate in a double well potential and for improving large $N$ expansions. We also discuss the problem of a charged Dirac particle in an external magnetic field and other potentials in terms of supersymmetric quantum mechanics. Finally, we discuss structures more general than supersymmetric quantum mechanics such as parasupersymmetric quantum mechanics in which there is a symmetry between a boson and a para-fermion of order $p$.