Continuum Shape Research Articles

THE MUSCLES TREASURY SURVEY. III. X-RAY TO INFRARED SPECTRA OF 11 M AND K STARS HOSTING PLANETS

ABSTRACT We present a catalog of panchromatic spectral energy distributions (SEDs) for 7 M and 4 K dwarf stars that span X-ray to infrared wavelengths (5 Å –5.5 μm). These SEDs are composites of Chandra or XMM-Newton data from 5–∼50 Å, a plasma emission model from ∼50–100 Å, broadband empirical estimates from 100–1170 Å, Hubble Space Telescope data from 1170–5700 Å, including a reconstruction of stellar Lyα emission at 1215.67 Å, and a PHOENIX model spectrum from 5700–55000 Å. Using these SEDs, we computed the photodissociation rates of several molecules prevalent in planetary atmospheres when exposed to each star’s unattenuated flux (“unshielded” photodissociation rates) and found that rates differ among stars by over an order of magnitude for most molecules. In general, the same spectral regions drive unshielded photodissociations both for the minimally and maximally FUV active stars. However, for O3 visible flux drives dissociation for the M stars whereas near-UV flux drives dissociation for the K stars. We also searched for an far-UV continuum in the assembled SEDs and detected it in 5/11 stars, where it contributes around 10% of the flux in the range spanned by the continuum bands. An ultraviolet continuum shape is resolved for the star ϵ Eri that shows an edge likely attributable to Si ii recombination. The 11 SEDs presented in this paper, available online through the Mikulski Archive for Space Telescopes, will be valuable for vetting stellar upper-atmosphere emission models and simulating photochemistry in exoplanet atmospheres.

Late-time spectroscopy of Type Iax Supernovae

We examine the late-time (t > 200 days after peak brightness) spectra of Type Iax supernovae (SNe Iax), a low-luminosity, low-energy class of thermonuclear stellar explosions observationally similar to, but distinct from, Type Ia supernovae. We present new spectra of SN 2014dt, resulting in the most complete published late-time spectral sequence of a SN Iax. At late times, SNe Iax have generally similar spectra, all with a similar continuum shape and strong forbidden-line emission. However, there is also significant diversity where some late-time SN Iax spectra display narrow P-Cygni features and a continuum indicative of a photosphere in addition to strong narrow forbidden lines, while others have no obvious P-Cygni features, strong broad forbidden lines, and weak narrow forbidden lines. Finally, some SNe Iax have spectra intermediate to these two varieties with weak P-Cygni features and broad/narrow forbidden lines of similar strength. We find that SNe Iax with strong broad forbidden lines also tend to be more luminous and have higher-velocity ejecta at peak brightness. We estimate blackbody and kinematic radii of the late-time photosphere, finding the latter an order of magnitude larger than the former. We propose a two-component model that solves this discrepancy and explains the diversity of the late-time spectra of SNe Iax. In this model, the broad forbidden lines originate from the SN ejecta, while the photosphere, P-Cygni lines, and narrow forbidden lines originate from a wind launched from the remnant of the progenitor white dwarf and is driven by the radioactive decay of newly synthesized material left in the remnant. The relative strength of the two components accounts for the diversity of late-time SN Iax spectra. This model also solves the puzzle of a long-lived photosphere and slow late-time decline of SNe Iax. (Abridged)

The thermal instability of the warm absorber in NGC 3783

We model the observed X-ray spectral continuum shape, ionic column densities, and absorption measure distribution (AMD) of the warm absorber in the Seyfert galaxy NGC 3783. We assume a photo-ionized medium with a uniform total (gas+radiation) pressure. The irradiation causes the wind to be radiation pressure compressed (RPC). We compare the observational characteristics derived from the 900 ksec Chandra observation to radiative transfer computations in pressure equilibrium using the radiative transfer code TITAN. We explore different values of the ionization parameter xi of the incident flux and adjust the hydrogen-equivalent column density, N_H0 of the warm absorber to match the observed soft X-ray continuum. We derive theoretical column densities for a broad range of ionic species of iron and neon and an AMD that we compare to the observations. We find an extension of the degeneracy between xi and N_H0 for the constant pressure models previously discussed for NGC 3783. Including the ionic column densities of iron and neon in the comparison between observations and data we conclude that a range of ionization parameters between 4000 and 8000 ergs cm/s is preferred. For the first time, we present theoretical AMD for a constant pressure wind in NGC 3783 that correctly reproduces the observed level and is in approximate agreement with the observational appearance of an instability region. Using a variety of observational indicators, we confirm that the X-ray outflow of NGC 3783 can be described as an RPC medium in pressure equilibrium. The observed AMD agrees with a uniformly hot or a uniformly cold thermal state. The measured ionic column densities suggest that the wind tends to the uniformly cold thermal state. The occurrence of thermal instability in the warm absorber model may depend on the computational method and the spatial scale on which the radiative transfer is solved.

Hole Sensitivity Analysis for Topology Optimization in Electrostatic System Using Virtual Hole Concept and Shape Sensitivity

In this paper, the hole sensitivity analysis is proposed for topology optimization in the electrostatic system. The hole sensitivity formula is derived using the virtual hole concept and the continuum shape sensitivity. Since the derived formula is expressed as a closed form, it is simple and easy to implement and apply to real problems. The hole sensitivity analysis provides information for topology variation. The geometrical expression of the topology variation is obtained using the level set method. Two numerical examples are tested to show the usefulness of the hole sensitivity analysis.

Two forms of continuum shape sensitivity method for fluid–structure interaction problems

Continuum Sensitivity Equation (CSE) methods for deriving and computing derivatives with respect to shape design variables are developed in two forms and compared in their application to fluid–structure interaction (FSI) problems. The local derivative form poses the CSEs in terms of the partial derivatives of the state variables with respect to shape parameters, while the CSEs in total derivative form are posed in terms of the total derivative, also known as the material or substantial derivative. In the literature CSEs are often posed in local form for fluids and total form for solids. The two forms are compared here for the purpose of applying a single form to both fluid and structure domains. The local form, also known as the boundary velocity method, requires design velocity only at the boundaries and interfaces of the domains to pose the CSEs. In contrast, the total form, also known as the domain velocity method, requires the design velocity in the whole domain. The local form requires higher-order spatial derivatives of the analysis solution than the total form, which affects the accuracy of its results. Higher order p-elements are shown to be a remedy to the inaccuracy of local form CSE seen in the literature for finite element solutions. The practicality, accuracy, and efficiency of these two CSE forms are compared based on the implementation and computed derivatives for three examples: a linear Timoshenko beam subject to a tip force, fluid flow around an airfoil, and an airfoil attached to a nonlinear joined beam subject to a gust load.

Convergence Study of Local Continuum Sensitivity Method Using Spatial Gradient Reconstruction

Gradient-based optimization for large-scale, multidisciplinary design problems requires accurate and efficient sensitivity analysis to compute design derivatives. The local continuum shape sensitivity method with spatial gradient reconstruction is a nonintrusive design sensitivity analysis method that was previously published for linear and nonlinear structural analysis. The numerical behavior, accuracy, and convergence for this method are studied using benchmark problems. The benchmark problems include finite element analysis of an axial bar and finite element analysis of a rectangular membrane using unstructured grids. The results are used to establish general guidelines for conducting mesh refinement and selecting the parameters associated with spatial gradient reconstruction calculations. Lastly, design derivatives are calculated for the potential flow solution around a Joukowsky airfoil. The analytic design derivative solutions are used to calculate the true errors associated with the local continuum design derivative results.

Spitzerinfrared spectrograph point source classification in the Small Magellanic Cloud

The Magellanic clouds are uniquely placed to study the stellar contribution to dust emission. Individual stars can be resolved in these systems even in the mid-infrared, and they are close enough to allow detection of infrared excess caused by dust.We have searched the Spitzer Space Telescope data archive for all Infrared Spectrograph (IRS) staring-mode observations of the Small Magellanic Cloud (SMC) and found that 209 Infrared Array Camera (IRAC) point sources within the footprint of the Surveying the Agents of Galaxy Evolution in the Small Magellanic Cloud (SAGE-SMC) Spitzer Legacy programme were targeted, within a total of 311 staring mode observations. We classify these point sources using a decision tree method of object classification, based on infrared spectral features, continuum and spectral energy distribution shape, bolometric luminosity, cluster membership and variability information. We find 58 asymptotic giant branch (AGB) stars, 51 young stellar objects (YSOs), 4 post-AGB objects, 22 Red Supergiants (RSGs), 27 stars (of which 23 are dusty OB stars), 24 planetary nebulae (PNe), 10Wolf-Rayet (WR) stars, 3 Hii regions, 3 R Coronae Borealis (R CrB) stars, 1 Blue Supergiant and 6 other objects, including 2 foreground AGB stars. We use these classifications to evaluate the success of photometric classification methods reported in the literature.

Inner-shell photoionization and core-hole decay of Xe and XeF2.

Photoionization cross sections and partial ion yields of Xe and XeF2 from Xe 3d(5/2), Xe 3d(3/2), and F 1s subshells in the 660-740 eV range are compared to explore effects of the F ligands. The Xe 3d-ϵf continuum shape resonances dominate the photoionization cross sections of both the atom and molecule, but prominent resonances appear in the XeF2 cross section due to nominal excitation of Xe 3d and F 1s electrons to the lowest unoccupied molecular orbital (LUMO), a delocalized anti-bonding MO. Comparisons of the ion products from the atom and molecule following Xe 3d photoionization show that the charge-state distribution of Xe ions is shifted to lower charge states in the molecule along with production of energetic F(+) and F(2+) ions. This suggests that, in decay of a Xe 3d core hole, charge is redistributed to the F ligands and the system dissociates due to Coulomb repulsion. The ion products from excitation of the F 1s-LUMO resonance are different and show strong increases in the yields of Xe(+) and F(+) ions. The subshell ionization thresholds, the LUMO resonance energies, and their oscillator strengths are calculated by relativistic coupled-cluster methods and agree well with measurements.

GRO J1744−28: an intermediate B-field pulsar in a low-mass X-ray binary

The bursting pulsar, GRO J1744-28, went again in outburst after $\sim$18 years of quiescence in mid-January 2014. We studied the broad-band, persistent, X-ray spectrum using X-ray data from a XMM-Newton observation, performed almost at the peak of the outburst, and from a close INTEGRAL observation, performed 3 days later, thus covering the 1.3-70.0 keV band. The spectrum shows a complex continuum shape that cannot be modelled with standard high-mass X-ray pulsar models, nor by two-components models. We observe broadband and peaked residuals from 4 to 15 keV, and we propose a self-consistent interpretation of these residuals, assuming they are produced by cyclotron absorption features and by a moderately smeared, highly ionized, reflection component. We identify the cyclotron fundamental at $\sim$ 4.7 keV, with hints for two possible harmonics at 10.4 keV and 15.8 keV. The position of the cyclotron fundamental allows an estimate for the pulsar magnetic field of (5.27 $\pm$ 0.06) $\times$ 10$^{11}$ G, if the feature is produced at its surface. From the dynamical and relativistic smearing of the disk reflected component, we obtain a lower limit estimate for the truncated accretion disk inner radius, ($\gtrsim$ 100 R$_g$), and for the inclination angle (18$^{\circ}$-48$^{\circ}$). We also detect the presence of a softer thermal component, that we associate with the emission from an accretion disk truncated at a distance from the pulsar of 50-115 R$_g$. From these estimates, we derive the magneto-spheric radius for disk accretion to be $\sim$ 0.2 times the classical Alfv\'en radius for radial accretion.

Derivation of Hole Sensitivity Formula for Topology Optimization in Magnetostatic System Using Virtual Hole Concept and Shape Sensitivity

In this paper, hole sensitivity formula is analytically derived for topology optimization in magnetostatic system. With the concept of virtual hole, the hole sensitivity is obtained using continuum shape sensitivity. To demonstrate the validity of the hole sensitivity, topology optimization of synchronous reluctance motor is tested. Since the hole sensitivity formula is represented only with electromagnetic field, it is easily implemented for topology optimization.

Continuum shape sensitivity analysis and what-if study for two-dimensional multi-scale crack propagation problems using bridging scale decomposition

This paper presents a shape sensitivity analysis and what-if study for two-dimensional multi-scale crack propagation problems using bridging scale decomposition. The sensitivity equations are derived in a continuum setting using direct differentiation method based on a continuum variational formulation of the bridging scale. Due to the fact that the crack propagation speed in an atomistic simulation is discrete in design, and cannot be formulated as a continuous function of shape design variables, we propose a hybrid method that combines analytical sensitivity analysis with finite difference approach. The finite difference part of the sensitivity analysis is only intended for calculating the sensitivity of crack growth speed based on the analytically obtained sensitivity coefficients of structural responses. The theoretical development on sensitivity formulation in this paper extends the application of the method to irregular-shaped finite elements and general design velocity fields. Furthermore, we evaluate and compare several performance measures that quantify crack propagation speed based on crack tip locations for sensitivity analysis and ultimately for structural optimization. A two-dimensional beam example is used to verify the accuracy of the proposed sensitivity approach. It is also demonstrated through a what-if study that with an adequate performance measure, the impact of macroscopic shape changes on microscopic crack propagation speed can be accurately predicted.

Local continuum shape sensitivity with spatial gradient reconstruction for nonlinear analysis

Gradient-based optimization for large-scale, multidisciplinary design problems requires accurate and efficient sensitivity analysis to compute design derivatives. Presented here is a nonintrusive analytic sensitivity method, that is relatively easy to implement. Furthermore, it can be as accurate as conventional analytic sensitivity methods, which are intrusive and tend to be difficult, if not infeasible, to implement. The nonintrusive local continuum shape sensitivity method with spatial gradient reconstruction (SGR) is formulated for nonlinear systems. This is an extension of the formulation previously published for linear systems. SGR, a numerical technique used to approximate spatial derivatives, can be leveraged to implement the sensitivity method in a nonintrusive manner. The method is used to compute design derivatives for a variety of applications, including nonlinear static beam bending, nonlinear transient gust response of a 2-D beam structure, and nonlinear static bending of rectangular plates. To demonstrate that the method is nonintrusive, all analyses are conducted using black box solvers. One limiting requirement of the method is that it requires the converged Jacobian or tangent stiffness matrix as output from the analysis tool. For each example the design derivatives of the structural displacement response are verified with finite difference calculations.

Local continuum shape sensitivity with spatial gradient reconstruction

Utilizing gradient-based optimization for large scale, multidisciplinary design problems requires accurate and efficient sensitivity or design derivative analysis. In general, numerical sensitivity methods, such as the finite difference method, are easy to implement but can be computationally expensive and inaccurate. In contrast, analytic sensitivity methods, such as the discrete and continuum methods, are highly accurate but can be very difficult, if not infeasible, to implement. A popular compromise is the semi-analytic method, but it too can be highly inaccurate when computing shape design derivatives. Presented here is an alternative method, which is easy to implement and can be as accurate as conventional analytic sensitivity methods. In this paper a general local continuum shape sensitivity method with spatial gradient reconstruction (SGR) is formulated. It is demonstrated that SGR, a numerical technique, can be used to solve the continuous sensitivity equations (CSEs) in a non-intrusive manner. The method is used to compute design derivatives for a variety of applications, including linear static beam bending, linear transient gust analysis of a 2-D beam structure, linear static bending of rectangular plates, and linear static bending of a beam-stiffened plate. Analysis is conducted with Nastran, and both displacement and stress design derivative solutions are presented. For each example the design derivatives are validated with either analytic or finite difference solutions.

The virtual observatory service TheoSSA: Establishing a database of synthetic stellar flux standards

In the framework of the Virtual Observatory (VO), the German Astrophysical Virtual Observatory (GAVO) developed the registered service TheoSSA (Theoretical Stellar Spectra Access). It provides easy access to stellar spectral energy distributions (SEDs) and is intended to ingest SEDs calculated by any model-atmosphere code, generally for all effective temperature, surface gravities, and elemental compositions. We will establish a database of SEDs of flux standards that are easily accessible via TheoSSA's web interface. The OB-type subdwarf Feige 110 is a standard star for flux calibration. State-of-the-art non-local thermodynamic equilibrium (NLTE) stellar-atmosphere models that consider opacities of species up to trans-iron elements will be used to provide a reliable synthetic spectrum to compare with observations. In case of Feige 110, we demonstrate that the model reproduces not only its overall continuum shape from the far-ultraviolet (FUV) to the optical wavelength range but also the numerous metal lines exhibited in its FUV spectrum. We present a state-of-the-art spectral analysis of Feige 110. We determined $T_\mathrm{eff} = 47\,250 \pm 2000\,\mathrm{K}$, $\log g = 6.00 \pm 0.20$ and the abundances of He, N, P, S, Ti, V, Cr, Mn, Fe, Co, Ni, Zn, and Ge. Ti, V, Mn, Co, Zn, and Ge were identified for the first time in this star. Upper abundance limits were derived for C, O, Si, Ca, and Sc. The TheoSSA database of theoretical SEDs of stellar flux standards guarantees that the flux calibration of astronomical data and cross-calibration between different instruments can be based on models and SEDs calculated with state-of-the-art model-atmosphere codes.

A Global Heat Compliance Measure Based Topology Optimization for the Transient Heat Conduction Problem

For topology optimization with transient loads, heat compliance varies with transient heat analysis. The peak value of the transient heat compliance should be minimized. Thus, this article proposes a global heat compliance measure to handle this kind of topology optimization for the transient heat conduction problem. The optimization model is then constructed by the global heat compliance measure. The finite-element, equivalent static loads, and continuum shape based sensitivity analyses are derived using the adjoint variable method. Through case studies, the effectiveness of the proposed global heat compliance measure for the transient heat conduction topology optimization is validated.

C II] AND12CO(1-0) EMISSION MAPS IN HLSJ091828.6+514223: A STRONGLY LENSED INTERACTING SYSTEM ATz= 5.24

We present Submillimeter Array (SMA) [CII] 158um and Jansky Very Large Array (JVLA) $^{12}$CO(1-0) line emission maps for the bright, lensed, submillimeter source at $z=5.2430$ behind Abell 773: HLSJ091828.6+514223 (HLS0918). We combine these measurements with previously reported line profiles, including multiple $^{12}$CO rotational transitions, [CI], water and [NII], providing some of the best constraints on the properties of the interstellar medium (ISM) in a galaxy at $z>5$. HLS0918 has a total far-infrared (FIR) luminosity L_FIR(8-1000um) = (1.6$\pm$0.1)x10^14 L_sun/mu, where the total magnification mu_total = 8.9$\pm$1.9, via a new lens model from the [CII] and continuum maps. Despite a HyLIRG luminosity, the FIR continuum shape resembles that of a local LIRG. We simultaneously fit all of the observed spectral line profiles, finding four components which correspond cleanly to discrete spatial structures identified in the maps. The two most redshifted spectral components occupy the nucleus of a massive galaxy, with a source plane separation <1 kpc. The reddest dominates the continuum map (de-magnified L_FIR = (1.1$\pm$0.2)x10^13 L_sun), and excites strong water emission in both nuclear components via a powerful FIR radiation field from the intense star formation. A third star-forming component is most likely a region of a merging companion (dV ~ 500 km/s) exhibiting generally similar gas properties. The bluest component originates from a spatially distinct region, and photo-dissociation region (PDR) analysis suggests that it is lower density, cooler and forming stars less vigorously than the other components. Strikingly, it has very strong [NII] emission which may suggest an ionized, molecular outflow. This comprehensive view of gas properties and morphology in HLS0918 previews the science possible for a large sample of high-redshift galaxies once ALMA attains full sensitivity.

The virtual observatory service TheoSSA: Establishing a database of synthetic stellar flux standards

H-rich, DA-type white dwarfs are particularly suited as primary standard stars for flux calibration. State-of-the-art NLTE models consider opacities of species up to trans-iron elements and provide reliable synthetic stellar-atmosphere spectra to compare with observation. We establish a database of theoretical spectra of stellar flux standards that are easily accessible via a web interface. In the framework of the Virtual Observatory, the German Astrophysical Virtual Observatory developed the registered service TheoSSA. It provides easy access to stellar spectral energy distributions (SEDs) and is intended to ingest SEDs calculated by any model-atmosphere code. In case of the DA white dwarf G 191-B2B, we demonstrate that the model reproduces not only its overall continuum shape but also the numerous metal lines exhibited in its ultraviolet spectrum. TheoSSA is in operation and contains presently a variety of SEDs for DA white dwarfs. It will be extended in the near future and can host SEDs of all primary and secondary flux standards. The spectral analysis of G 191-B2B has shown that our hydrostatic models reproduce the observations best at an effective temperature of 60000 +/- 2000K and a surface gravity of log g = 7.60 +/- 0.05. We newly identified Fe VI, Ni VI, and Zn IV lines. For the first time, we determined the photospheric zinc abundance with a logarithmic mass fraction of -4.89 (7.5 times solar). The abundances of He (upper limit), C, N, O, Al, Si, O, P, S, Fe, Ni, Ge, and Sn were precisely determined. Upper abundance limits of 10% solar were derived for Ti, Cr, Mn, and Co. The TheoSSA database of theoretical SEDs of stellar flux standards guarantees that the flux calibration of all astronomical data and cross-calibration between different instruments can be based on the same models and SEDs calculated with different model-atmosphere codes and are easy to compare.

A generic π* shape resonance observed in energy-dependent photoelectron angular distributions from two-colour, resonant multiphoton ionization of difluorobenzene isomers

We present new evidence for the existence of a near threshold π* shape resonance as a common feature in the photoionization of each isomer of difluorobenzene. Experimentally, this is revealed by significant changes in the anisotropy of the photoelectron angular distributions (PADs) following the ionization of the optically aligned S1 state of these molecules at varying photon energies. Continuum multiple scattering Xα calculations reproduce this behaviour well, and allow the visualisation of the continuum shape resonances. The resonances are unusually narrow in energy (<1 eV), but nevertheless appear to extend right down to the ionization thresholds--exactly the low energy range typically accessed in laser-based resonance enhanced multiphoton ionization (REMPI) schemes. The anticipation of such pronounced energy dependence in the PADs and cross-sections sought for other molecules, and an ability to accurately predict such features, should be important for the reliable application and interpretation of experiments involving REMPI probing of those molecules.

Near-infrared spectroscopy of SN 2009ip's 2012 brightening reveals a dusty pre-supernova environment

Abstract We present low-resolution near-infrared (IR) 0.8–2.5 μm spectra of Supernova (SN) 2009ip, taken immediately before, during and just after its rapid brightening in late September/October 2012. The first epoch shows the same general spectral characteristics as the later epochs (smooth continuum, narrow H and He i emission lines), but the IR continuum shape is substantially redder than the later epochs. The epoch 1 continuum can be approximated by reddening the peak-luminosity (epoch 3) spectrum by E(B − V) = 1.0 mag, but the blue colour seen in visual-wavelength spectra at the same time indicates that strong wavelength-dependent extinction by circumstellar dust is not the correct explanation. Instead, we favour the hypothesis that the redder colour before the brightening arises from excess emission from hot ∼2000 K circumstellar dust. The minimum radius ( ≳ 120 au) deduced from the dust temperature and observed luminosity of the transient, combined with the observed expansion speed in the precursor outbursts of SN 2009ip, is consistent with an ejection at least 1.1 yr earlier. The mass of hot dust indicated by the IR excess is ∼4 × 10−7 M⊙, although this is only a lower limit since the near-IR data do not constrain the mass of cooler dust. Thus, the observed pre-SN outbursts of this object were able to efficiently form dust into which the SN ejecta and radiation now propagate. This is consistent with the notion that the same pre-SN eruptions that generally give rise to SNe IIn also give rise to the dust needed for their commonly observed IR echoes. We also discuss some aspects of the IR line profiles, including He i λ10 830.

A NEAR-INFRARED SPECTROSCOPIC STUDY OF YOUNG FIELD ULTRACOOL DWARFS

We present a near-infrared (0.9-2.4 microns) spectroscopic study of 73 field ultracool dwarfs having spectroscopic and/or kinematic evidence of youth (~10-300 Myr). Our sample is composed of 48 low-resolution (R~100) spectra and 41 moderate-resolution spectra (R>~750-2000). First, we establish a method for spectral typing M5-L7 dwarfs at near-IR wavelengths that is independent of gravity. We find that both visual and index-based classification in the near-IR provide consistent spectral types with optical spectral types, though with a small systematic offset in the case of visual classification at J and K band. Second, we examine features in the spectra of ~10 Myr ultracool dwarfs to define a set of gravity-sensitive indices based on FeH, VO, K, Na and H-band continuum shape. We then create an index-based method for classifying the gravities of M6-L5 dwarfs that provides consistent results with gravity classifications from optical spectroscopy. Our index-based classification can distinguish between young and dusty objects. Guided by the resulting classifications, we propose a set of low-gravity spectral standards for the near-IR. Finally, we estimate the ages corresponding to our gravity classifications.