1,292 research outputs found

    Sparse Bayesian Nonlinear System Identification using Variational Inference

    Get PDF
    IEEE Bayesian nonlinear system identification for one of the major classes of dynamic model, the nonlinear autoregressive with exogenous input (NARX) model, has not been widely studied to date. Markov chain Monte Carlo (MCMC) methods have been developed, which tend to be accurate but can also be slow to converge. In this contribution, we present a novel, computationally efficient solution to sparse Bayesian identification of the NARX model using variational inference, which is orders of magnitude faster than MCMC methods. A sparsity-inducing hyper-prior is used to solve the structure detection problem. Key results include: 1. successful demonstration of the method on low signal-to-noise ratio signals (down to 2dB); 2. successful benchmarking in terms of speed and accuracy against a number of other algorithms: Bayesian LASSO, reversible jump MCMC, forward regression orthogonalisation, LASSO and simulation error minimisation with pruning; 3. accurate identification of a real world system, an electroactive polymer; and 4. demonstration for the first time of numerically propagating the estimated nonlinear time-domain model parameter uncertainty into the frequency-domain

    Sparse Bayesian Nonlinear System Identification using Variational Inference

    Get PDF
    IEEE Bayesian nonlinear system identification for one of the major classes of dynamic model, the nonlinear autoregressive with exogenous input (NARX) model, has not been widely studied to date. Markov chain Monte Carlo (MCMC) methods have been developed, which tend to be accurate but can also be slow to converge. In this contribution, we present a novel, computationally efficient solution to sparse Bayesian identification of the NARX model using variational inference, which is orders of magnitude faster than MCMC methods. A sparsity-inducing hyper-prior is used to solve the structure detection problem. Key results include: 1. successful demonstration of the method on low signal-to-noise ratio signals (down to 2dB); 2. successful benchmarking in terms of speed and accuracy against a number of other algorithms: Bayesian LASSO, reversible jump MCMC, forward regression orthogonalisation, LASSO and simulation error minimisation with pruning; 3. accurate identification of a real world system, an electroactive polymer; and 4. demonstration for the first time of numerically propagating the estimated nonlinear time-domain model parameter uncertainty into the frequency-domain

    Transient HDO rovibrational satellite peaks in solid parahydrogen: evidence of hydrogen atoms or vacancies?

    No full text
    We present FTIR studies of the 193 nm photolysis of fully deuterated formic acid (DCOOD) isolated in solid parahydrogen at 1.9 K which show evidence of transient HDO rovibrational satellite peaks. The S1 and S2 satellite peaks are readily detected for α-type (1₀₁ ← 0₀₀) rovibrational transitions of HDO either during or immediately after photolysis. Intensity measurements show the HDO b-type (1₁₁ ← 0₀₀) rovibrational transitions have satellite peaks as well, but due to the greater linewidth of these absorptions, the satellite peaks cannot be spectroscopically resolved from the monomer transition and are therefore difficult to detect. These newly identified HDO satellite peaks may result from the HDO photoproduct being formed next to an H atom or a vacancy in the parahydrogen solid. The development of the infrared spectroscopy of these satellite peaks can provide a new means to study radiation effects on low-temperature hydrogen solids doped with chemical species

    Frequency-domain analysis for nonlinear systems with time-domain model parameter uncertainty

    Get PDF
    Frequency-domain analysis of dynamic systems is important across many areas of engineering. However, whilst there are many analysis methods for linear systems, the problem is much less widely studied for nonlinear systems. Frequency-domain analysis of nonlinear systems using frequency response functions (FRFs) is particularly important to reveal resonances, super/sub-harmonics and energy transfer across frequencies. In this paper the novel contribution is a time-domain model-based approach to describing the uncertainty of nonlinear systems in the frequency-domain. The method takes a nonlinear input-output model that has normally distributed parameters, and propagates that uncertainty into the frequency-domain using analytic expressions based on FRFs. We demonstrate the approach on both synthetic examples of nonlinear systems and a real-world nonlinear system identified from experimental data. We benchmark the proposed approach against a brute-force technique based on Monte Carlo sampling and show that there is good agreement between the methods

    Trapping of Projectiles in Fixed Scatterer Calculations

    Full text link
    We study multiple scattering off nuclei in the closure approximation. Instead of reducing the dynamics to one particle potential scattering, the scattering amplitude for fixed target configurations is averaged over the target groundstate density via stochastic integration. At low energies a strong coupling limit is found which can not be obtained in a first order optical potential approximation. As its physical explanation, we propose it to be caused by trapping of the projectile. We analyse this phenomenon in mean field and random potential approximations. (PACS: 24.10.-i)Comment: 15 page

    Radiative charge transfer lifetime of the excited state of (NaCa)+^+

    Get PDF
    New experiments were proposed recently to investigate the regime of cold atomic and molecular ion-atom collision processes in a special hybrid neutral-atom--ion trap under high vacuum conditions. The collisional cooling of laser pre-cooled Ca+^+ ions by ultracold Na atoms is being studied. Modeling this process requires knowledge of the radiative lifetime of the excited singlet A1Σ+^1\Sigma^+ state of the (NaCa)+^+ molecular system. We calculate the rate coefficient for radiative charge transfer using a semiclassical approach. The dipole radial matrix elements between the ground and the excited states, and the potential curves were calculated using Complete Active Space Self-Consistent field and M\"oller-Plesset second order perturbation theory (CASSCF/MP2) with an extended Gaussian basis, 6-311+G(3df). The semiclassical charge transfer rate coefficient was averaged over a thermal Maxwellian distribution. In addition we also present elastic collision cross sections and the spin-exchange cross section. The rate coefficient for charge transfer was found to be 2.3×10162.3\times 10^{-16} cm3^3/sec, while those for the elastic and spin-exchange cross sections were found to be several orders of magnitude higher (1.1×1081.1\times 10^{-8} cm3^3/sec and 2.3×1092.3\times 10^{-9} cm3^3/sec, respectively). This confirms our assumption that the milli-Kelvin regime of collisional cooling of calcium ions by sodium atoms is favorable with the respect to low loss of calcium ions due to the charge transfer.Comment: 4 pages, 5 figures; v.2 - conceptual change

    Bcc 4^4He as a Coherent Quantum Solid

    Full text link
    In this work we investigate implications of the quantum nature of bcc 4^{4}% He. We show that it is a unique solid phase with both a lattice structure and an Off-Diagonal Long Range Order of coherently oscillating local electric dipole moments. These dipoles arise from the local motion of the atoms in the crystal potential well, and oscillate in synchrony to reduce the dipolar interaction energy. The dipolar ground-state is therefore found to be a coherent state with a well defined global phase and a three-component complex order parameter. The condensation energy of the dipoles in the bcc phase stabilizes it over the hcp phase at finite temperatures. We further show that there can be fermionic excitations of this ground-state and predict that they form an optical-like branch in the (110) direction. A comparison with 'super-solid' models is also discussed.Comment: 12 pages, 8 figure

    Cardiovascular disease diagnoses among older women with endometrial cancer

    Get PDF
    Background: Endometrial cancer (EC) shares risk factors (e.g. obesity) with cardiovascular disease (CVD), yet little research has investigated CVD diagnoses among EC survivors. We aimed to describe the burden of CVD diagnoses among older women with EC compared to women without a cancer history. Methods: Women aged 66+ years with an EC diagnosis during 2004–2017 (N = 44,386) and matched women without cancer (N = 221,219) were identified in the SEER-Medicare linked data. An index date was defined as the cancer diagnosis date of the EC case in that matched set. ICD-9/10 diagnosis codes were used to define CVD outcomes in the Medicare claims. Prevalent CVD was identified using diagnosis codes in the year before the index date. Hazard ratios (HRs) for incident CVD diagnoses after the index date were estimated using multivariable Cox proportional hazards regression. Women with a prevalent CVD were excluded from incidence analyses for that outcome. Results: Compared to women without cancer, women with EC had a higher prevalence of CVD diagnoses at the index date. In analyses beginning follow-up at 1 year post-index date, EC survivors had an increased risk of incident CVD diagnoses including ischemic heart diseases (HR = 1.73; 95% CI: 1.69–1.78), pulmonary heart disease (HR = 1.95; 95% CI: 1.88–2.02), and diseases of the veins and lymphatics (HR = 2.71; 95% CI: 95% CI: 2.64–2.78). Risk of CVD diagnoses among women with EC was also elevated within the first year post-index date. Conclusions: Management of pre-existing CVD and monitoring for incident CVD may be critical during EC treatment and throughout long-term survivorship

    Species abundance dynamics under neutral assumptions: a Bayesian approach to the controversy

    Get PDF
    1. Hubbell's 'Unified Neutral Theory of Biodiversity and Biogeography' (UNTB) has generated much controversy about both the realism of its assumptions and how well it describes the species abundance dynamics in real communities. 2. We fit a discrete-time version of Hubbell's neutral model to long-term macro-moth (Lepidoptera) community data from the Rothamsted Insect Survey (RIS) light-traps network in the United Kingdom. 3. We relax the assumption of constant community size and use a hierarchical Bayesian approach to show that the model does not fit the data well as it would need parameter values that are impossible. 4. This is because the ecological communities fluctuate more than expected under neutrality. 5. The model, as presented here, can be extended to include environmental stochasticity, density-dependence, or changes in population sizes that are correlated between different species
    corecore