A multidimensionally consistent version of Hirota's discrete KdV equation

A multidimensionally consistent generalisation of Hirota's discrete KdV equation is proposed, it is a quad equation defined by a polynomial that is quadratic in each variable. Soliton solutions and interpretation of the model as superposition principle are given. It is discussed how an important property of the defining polynomial, a factorisation of discriminants, appears also in the few other known discrete integrable multi-quadratic models.Comment: 11 pages, 2 figure

Towards exponentially-convergent simulations of extreme-mass-ratio inspirals: A time-domain solver for the scalar Teukolsky equation with singular source terms

Gravitational wave signals from extreme mass ratio inspirals are a key target for space-based gravitational wave detectors. These systems are typically modeled as a distributionally-forced Teukolsky equation, where the smaller black hole is treated as a Dirac delta distribution. Time-domain solvers often use regularization approaches that approximate the Dirac distribution that often introduce small length scales and are a source of systematic error, especially near the smaller black hole. We describe a multi-domain discontinuous Galerkin method for solving the distributionally-forced Teukolsky equation that describes scalar fields evolving on a Kerr spacetime. To handle the Dirac delta, we expand the solution in spherical harmonics and recast the sourced Teukolsky equation as a first-order, one-dimensional symmetric hyperbolic system. This allows us to derive the method's numerical flux to correctly account for the Dirac delta. As a result, our method achieves global spectral accuracy even at the source's location. To connect the near field to future null infinity, we use the hyperboloidal layer method, allowing us to supply outer boundary conditions and providing direct access to the far-field waveform. We document several numerical experiments where we test our method, including convergence tests against exact solutions, energy luminosities for circular orbits, the scheme's superconvergence properties at future null infinity, and the late-time tail behavior of the scalar field. We also compare two systems that arise from different choices of the first-order reduction variables, finding that certain choices are numerically problematic in practice. The methods developed here may be beneficial when computing gravitational self-force effects, where the regularization procedure has been developed for the spherical harmonic modes and high accuracy is needed at the Dirac delta's location.Comment: 20 pages, 7 figures and 2 table

Validation of gyrokinetic modelling of light impurity transport including rotation in ASDEX Upgrade

Upgraded spectroscopic hardware and an improved impurity concentration calculation allow accurate determination of boron density in the ASDEX Upgrade tokamak. A database of boron measurements is compared to quasilinear and nonlinear gyrokinetic simulations including Coriolis and centrifugal rotational effects over a range of H-mode plasma regimes. The peaking of the measured boron profiles shows a strong anti-correlation with the plasma rotation gradient, via a relationship explained and reproduced by the theory. It is demonstrated that the rotodiffusive impurity flux driven by the rotation gradient is required for the modelling to reproduce the hollow boron profiles at higher rotation gradients. The nonlinear simulations validate the quasilinear approach, and, with the addition of perpendicular flow shear, demonstrate that each symmetry breaking mechanism that causes momentum transport also couples to rotodiffusion. At lower rotation gradients, the parallel compressive convection is required to match the most peaked boron profiles. The sensitivities of both datasets to possible errors is investigated, and quantitative agreement is found within the estimated uncertainties. The approach used can be considered a template for mitigating uncertainty in quantitative comparisons between simulation and experiment.Comment: 19 pages, 11 figures, accepted in Nuclear Fusio

Hydrodynamic reductions of the heavenly equation

We demonstrate that Pleba\'nski's first heavenly equation decouples in infinitely many ways into a triple of commuting (1+1)-dimensional systems of hydrodynamic type which satisfy the Egorov property. Solving these systems by the generalized hodograph method, one can construct exact solutions of the heavenly equation parametrized by arbitrary functions of a single variable. We discuss explicit examples of hydrodynamic reductions associated with the equations of one-dimensional nonlinear elasticity, linearly degenerate systems and the equations of associativity.Comment: 14 page

Third-order nonlinear optical properties of stacked bacteriochlorophylls in bacterial photosynthetic light-harvesting proteins

Enhancement of the nonresonant second order molecular hyperpolarizabilities {gamma} were observed in stacked macrocyclic molecular systems, previously in a {micro}-oxo silicon phthalocyanine (SiPcO) monomer, dimer and trimer series, and now in bacteriochlorophyll a (BChla) arrays of light harvesting (LH) proteins. Compared to monomeric BChla in a tetrahydrofuran (THF) solution, the <{gamma}> for each macrocycle was enhanced in naturally occurring stacked macrocyclic molecular systems in the bacterial photosynthetic LH proteins where BChla`s are arranged in tilted face-to-face arrays. In addition, the {gamma} enhancement is more significant in B875 of LH1 than in B850 in LH2. Theoretical modeling of the nonresonant {gamma} enhancement using simplified molecular orbitals for model SiPcO indicated that the energy level of the two photon state is crucial to the {gamma} enhancement when a two photon process is involved, whereas the charge transfer between the monomers is largely responsible when one photon near resonant process is involved. The calculated results can be extended to {gamma} enhancement in B875 and B850 arrays, suggesting that BChla in B875 are more strongly coupled than in B850. In addition, a 50--160 fold increase in <{gamma}> for the S{sub 1} excited state of relative to S{sub 0} of bacteriochlorophyll in vivo was observed which provides an alternative method for probing excited state dynamics and a potential application for molecular switching

Computer simulation and analysis of hemodynamic changes in abdominal aortic aneurysms treated with fenestrated endovascular grafts

The purpose of this study was to perform a simulation of blood flow and analyze the hemodynamic changes in patients with abdominal aortic aneurysms (AAA) treated with fenestrated stent grafts. Four patients with AAA undergoing multislice computed tomography angiography pre-and post-fenestrated stent graft implantation were selected for inclusion in the study. Geometric models and hexahedral volume meshes were successfully generated for pre- and post-stent fenestrated implantation. The blood flow pattern was simulated inside the abdominal aortic aneurysm and arterial branches, as well as with a stentgraft in situ. Flow visualization showed that flow disturbances inside the aneurysm were apparently decreased and flow rate was not affected significantly at the renal arteries after deployment of the fenestrated stents into these branches. The wall pressure was found to reduce inside the aneurysm sac following implantation of stent grafts. In this preliminary study, we successfully simulated the flow characteristics in abdominal aortic aneurysm before and after fenestrated endovascular repair