STEADY-STATE SOLUTIONS IN NONLINEAR DIFFUSIVE SHOCK ACCELERATION

Stationary solutions to the equations of nonlinear diffusive shock acceleration play a fundamental role in the theory of cosmic-ray acceleration. Their existence usually requires that a fraction of the accelerated particles be allowed to escape from the system. Because the scattering mean free path is thought to be an increasing function of energy, this condition is conventionally implemented as an upper cutoff in energy space-particles are then permitted to escape from any part of the system, once their energy exceeds this limit. However, because accelerated particles are responsible for the substantial amplification of the ambient magnetic field in a region upstream of the shock front, we examine an alternative approach in which particles escape over a spatial boundary. We use a simple iterative scheme that constructs stationary numerical solutions to the coupled kinetic and hydrodynamic equations. For parameters appropriate for supernova remnants, we find stationary solutions with efficient acceleration when the escape boundary is placed at the point where growth and advection of strongly driven nonresonant waves are in balance. We also present the energy dependence of the distribution function close to the energy where it cuts off-a diagnostic that is in principle accessible to observation.

Nonlinear optics from an ab initio approach by means of the dynamical Berry phase: Application to second- and third-harmonic generation in semiconductors

We present an ab initio real-time-based computational approach to study nonlinear optical properties in condensed matter systems that is especially suitable for crystalline solids and periodic nanostructures. The equations of motion and the coupling of the electrons with the external electric field are derived from the Berry-phase formulation of the dynamical polarization [Souza et al., Phys. Rev. B 69, 085106 (2004)]. Many-body effects are introduced by adding single-particle operators to the independent-particle Hamiltonian. We add a Hartree operator to account for crystal local effects and a scissor operator to correct the independent particle band structure for quasiparticle effects. We also discuss the possibility of accurately treating excitonic effects by adding a screened Hartree-Fock self-energy operator. The approach is validated by calculating the second-harmonic generation of SiC and AlAs bulk semiconductors: an excellent agreement is obtained with existing ab initio calculations from response theory in frequency domain [Luppi et al., Phys. Rev. B 82, 235201 (2010)]. We finally show applications to the second-harmonic generation of CdTe and the third-harmonic generation of Si. 

Sound Synthesis for Contact-Driven Musical Instruments via Discretisation of Hamilton's Equations

Physical modelling of musical instruments involves studying nonlinear interactions between parts of the instrument. These can pose several difficulties concerning the accuracy and stability of numerical algorithms. In particular, when the underlying forces are non-analytic functions of the phase-space variables, a stability proof can only be obtained in limited cases. An approach has been recently presented by the authors, leading to unconditionally stable simulations for lumped collision models. In that study, discretisation of Hamilton’s equations instead of the usual Newton’s equation of motion yields a numerical scheme that can be proven to be energy conserving. In this paper, the above approach is extended to collisions of distributed objects. Namely, the interaction of an ideal string with a flat barrier is considered. The problem is formulated within the Hamiltonian framework and subsequently discretised. The resulting nonlinearmatrix equation can be shown to possess a unique solution, that enables the update of the algorithm. Energy conservation and thus numerical stability follows in a way similar to the lumped collision model. The existence of an analytic description of this interaction allows the validation of the model’s accuracy. The proposed methodology can be used in sound synthesis applications involving musical instruments where collisions occur either in a confined (e.g. hammer-string interaction, mallet impact) or in a distributed region (e.g. string-bridge or reed-mouthpiece interaction).

Transonic Nonlinear Aeroelastic Simulations Using An Harmonic Balance Method

This paper presents an approach to compute transonic Limit Cycle O
scillations using a coupled Harmonic Balance formulation based on the Euler equations for fluid dynamics and finite element models. The paper will investigate the role of aerodynamic (shocks) and structural nonlinearities in driving the limit cycle behaviour. Part icular attention will be given to nonlinear interactions for subcritical LCOs. The Aero elastic Harmonic Balance formulation, allows for solutions of the coupled structural dynamics and CFD system at a reduced cost.

Performance of polarisation functionals for linear and nonlinear optical properties of bulk zinc chalcogenides ZnX (X=S,Se,Te).

We calculated the frequency dependent macroscopic dielectric function and second-harmonic generation of cubic ZnS, ZnSe and ZnTe within time-dependent density-polarisation functional theory. The macroscopic dielectric function is calculated in a linear response framework, and second-harmonic generation in a real-time framework. The macroscopic exchange–correlation electric field that enters the time-dependent Kohn–Sham equations and accounts for long range correlation is approximated as a simple polarisation functional αP, where P is the macroscopic polarisation. Expressions for α are taken from the recent literature. The performance of the resulting approximations for the exchange–correlation electric field is analysed by comparing the theoretical spectra with experimental results and results obtained at the levels of the independent particle approximation and the random-phase approximation. For the dielectric function we also compare with state-of-the art calculations at the level of the Bethe–Salpeter equation.

Symmetries and exact solutions of KP equation with an arbitrary nonlinear term

The generalized KP (GKP) equations with an arbitrary nonlinear term model and characterize many nonlinear physical phenomena. The symmetries of GKP equation with an arbitrary nonlinear term are obtained. The condition that must satisfy for existence the symmetries group of GKP is derived and also the obtained symmetries are classified according to different forms of the nonlinear term. The resulting similarity reductions are studied by performing the bifurcation and the phase portrait of GKP and also the corresponding solitary wave solutions of GKP
equation are constructed.