One-zone Model Research Articles

Experimental and Numerical Evaluation of an HCCI Engine Fueled with Biogas for Power Generation under Sub-Atmospheric Conditions

Energy transition to renewable sources and more efficient technologies is needed for sustainable development. Although this transition is expected to take a longer time in developing countries, strategies that have been widely explored by the international academic community, such as advanced combustion modes and microgeneration, could be implemented more easily. However, the implementation of these well-known strategies in developing countries requires in-depth research because of the specific technical, environmental, social, and economic conditions. The present research relies on the use of biogas-fueled HCCI engines for power generation under sub-atmospheric conditions provided by high altitudes above sea level in Colombia. A small air-cooled commercial Diesel engine was modified to run in HCCI combustion mode by controlling the air–biogas mixture temperature using an electric heater at a high speed of 1800 revolutions per minute. An experimental setup was implemented to measure and control the most important experimental variables, such as engine speed, biogas flow rate, intake temperature, crank angle degree, intake pressure, NOx emissions, and in-cylinder pressure. High intake temperature requirements of around 320 ∘C were needed to achieve stable HCCI combustion; the maximum net indicated mean effective pressure (IMEPn) was around 1.5 bar, and the highest net indicated efficiency was close to 32%. Higher intake pressures and the addition of ozone to the intake mixture were numerically studied as ways to reduce the intake temperature requirements for stable HCCI combustion and improve engine performance. These strategies were studied using a one-zone model along with detailed chemical kinetics, and the model was adjusted using the experimental results. The simulation results showed that the addition of 500 ppm of ozone could reduce the intake temperature requirements by around 50 ∘C. The experimental and numerical results achieved in this research are important for the design and implementation of HCCI engines running biogas for microgeneration systems in developing countries which exhibit more difficult conditions for HCCI combustion implementation.

General Physical Properties of Fermi Blazars

We study the general physical properties of Fermi blazars using the Fermi fourth source catalog data (4FGL-DR2). The quasi-simultaneous multiwavelength data of Fermi blazars are fitted by using the one-zone leptonic model to obtain some physical parameters, such as jet power, magnetic field, and Doppler factor. We study the distributions of the derived physical parameters as a function of black hole mass and accretion disk luminosity. The main results are as follows. (1) For a standard thin accretion disk, the jet kinetic power of most flat-spectrum radio quasars can be explained by the Blandford–Payne (BP) mechanism. However, the jet kinetic power of most BL Lacertae objects (BL Lacs) cannot be explained by either the Blandford–Znajek mechanism or the BP mechanism. The BL Lacs may have advection-dominated accretion flows surrounding their massive black holes. (2) After excluding the redshift, there is a moderately strong correlation between the jet kinetic power and jet radiation power and the accretion disk luminosity for Fermi blazars. These results confirm a close connection between jet and accretion. The jet kinetic power is slightly larger than the accretion disk luminosity for Fermi blazars. (3) There is a significant correlation between jet kinetic power and gamma-ray luminosity and radio luminosity for Fermi blazars, which suggests that gamma-ray luminosity and radio luminosity can be used to indicate the jet kinetic power.

EOS-dependent millihertz quasi-periodic oscillation in low-mass X-ray binary

ABSTRACT We studied the frequency and critical mass accretion rate of millihertz quasi-periodic oscillations (mHz QPOs) using a one-zone X-ray burst model. The surface gravity is specified by two kinds of equation of states: neutron star (NS) and strange star (SS). The base flux, Qb, is set in the range of 0–2 MeV nucleon−1. It is found that the frequency of mHz QPO is positively correlated to the surface gravity but negatively to the base heating. The helium mass fraction has a significant influence on the oscillation frequency and luminosity. The observed 7–9 mHz QPOs can be either explained by a heavy NS/light SS with a small base flux or a heavy SS with a large base flux. As base flux increases, the critical mass accretion rate for marginally stable burning is found to be lower. Meanwhile, the impact of metallicity on the properties of mHz QPOs was investigated using one-zone model. It shows that both the frequency and critical mass accretion rate decrease as metallicity increases. An accreted NS/SS with a higher base flux and metallicity, combined with a lower surface gravity and helium mass fraction, could be responsible for the observed critical mass accretion rate ($\dot{m}\simeq 0.3\dot{m}_{\rm Edd}$). The accreted fuel would be in stable burning if base flux is over than ∼2 MeV nucleon−1. This finding suggests that the accreting NSs/SSs in low-mass X-ray binaries showing no type I X-ray bursts possibly have a strong base heating.

A model for GN-z11: top-heavy stellar initial mass functions in forming galactic nuclei and ultra-compact dwarfs

ABSTRACT Recent JWST observations of the z = 10.6 galaxy GN-z11 have revealed a very high gas-phase nitrogen abundance (higher than four times the solar value), a very small half-light radius (≈60 pc), and a large stellar mass (Ms ≈ 109 M⊙) for its size. We consider that this object is a forming galactic nucleus or ultra-compact dwarf galaxy rather than a proto globular cluster, and thereby investigate the chemical abundance pattern using one-zone chemical evolution models. The principal results of the models are as follows. The observed log (N/O) > −0.24, log (C/O) > −0.78, and 12 + log (O/H) ≈ 7.8 can be self-consistently reproduced by the models both with very short star formation time-scales (<107 yr) and with top-heavy stellar initial mass functions. The adopted assumption of no chemical enrichment by massive (m > 25 M⊙) core collapse supernovae (CCSNe) is also important for the reproduction of high gas-phase log (N/O), because such CCSNe can decrease high log (N/O) of gas polluted by OB and Wolf–Rayet stars. GN-z11 can have a significant fraction (>0.5) of nitrogen-rich ([N/Fe] > 0.5) stars, which implies a possible link between nitrogen-rich stellar populations of the inner Galaxy and giant elliptical galaxies and high-z objects with high gas-phase log (N/O) like GN-z11.

Investigating the Multiband Nonthermal Radiative Properties of HESS J1420−607

HESS J1420−607 is a γ-ray emitting source associated with the pulsar wind nebula (PWN) powered by the energetic pulsar PSR J1420−6048. Based on 14 yr of data obtained with the Fermi Large Area Telescope, we re-analyzed its GeV γ-ray radiative properties, resulting in detailed spectra obtained within the band 10–200 GeV. Moreover, we use a one-zone time-dependent model for the multiband nonthermal emission from pulsar wind nebulae to investigate the radiative properties of the nebula associated with HESS J1420−607. Assuming that the electrons/positrons are injected into the PWN with a broken power law spectrum with indexes of ∼1.6 and ∼2.7, as well as a break Lorentz factor of ∼5 × 106, the results indicate that the multi-wavelength spectral energy distribution is consistent with the detected fluxes in both X-rays and γ-rays. The results support that the γ-rays detected can be produced via inverse Compton scattering of the energetic electrons/positrons within the PWN.

Analytical Constraints on the Radius and Bulk Lorentz Factor in the Lepto-Hadronic One-Zone Model of BL Lacs

In this work, we study the parameter space of neutrino-emitting BL Lacs under the framework of the one-zone lepto-hadronic model. We show that constraints on the model come from various aspects of observations such as the variability timescale of blazar flares, gamma ray opacity and the spectral energy distribution of electromagnetic emission, as well as the inferred neutrino emissivity of the blazar. We apply our method to two potential neutrino sources, i.e., TXS 0506+056 and PKS 0735+178, which are BL Lacs. Then, we explore and summarize the allowed range of parameters such as the bulk Lorentz factor and the blob radius under different distributions of injected protons. We find that the parameter space that is available to explain the BL Lac–neutrino association is sensitive to the proton distribution, and usually, an injected proton luminosity significantly exceeding the Eddington luminosity is required for both sources. Our results suggest that the simple lepto-hadronic one-zone model may not be a reasonable interpretation for BL Lac–neutrino associations.

Revisiting the Chandra Observation on the Region of PSR J1J1809–193: Indication of the Existence of an X-Ray Halo and Implication for the Origin of HESS J1809–193

HESS J1809–193 is an extended TeV γ-ray source and the origin of its γ-ray emission remains ambiguous. The pulsar wind nebula (PWN) of PSR J1J1809–193 lying inside the extended γ-ray emission is a possible candidate. Powered by the central pulsar, ultrarelativistic electrons in the PWN can produce radio to X-ray emission through synchrotron and γ-ray emission by inverse Compton (IC) scattering. To check whether this PWN is the counterpart of HESS J1809–193, we analyzed the Chandra X-ray radial intensity profile and spectral index profile of this PWN. We then adopted a one-zone isotropic diffusion model to fit the keV and TeV data. We found diffuse nonthermal X-ray emission extending beyond the PWN, which is likely an X-ray halo radiated by escaping electron/positron pairs from the PWN. A relatively strong magnetic field of ∼20 μG is required to explain the spatial evolution of the X-ray spectrum (i.e., the significant softening of the spectrum with increasing distance from the pulsar), which, however, would suppress the IC radiation of pairs. Our result implies that a hadronic component may be needed to explain HESS J1809–193.

Modelling the variable emission states of γ-ray-emitting narrow-line Seyfert 1 galaxies

ABSTRACT γ-ray-emitting narrow-line Seyfert 1 galaxies (γ-NLS1) constitute an intriguing small population of active galactic nuclei with γ-ray emission resembling low-power flat-spectrum radio quasars (FSRQ), but with differing physical properties. They are jetted, γ/radio-loud Seyfert galaxies, with relatively low black hole masses, accreting at exceptionally high, near-Eddington rates. Certain of these sources exhibit highly variable emission states on relatively short time-scales, the physical origin of which remains elusive. In this work, varying emission states of two bona fide NLS1s, 1H 0323+342 and PMN J0948+0022, and one little-studied FSRQ/intermediate object, B2 0954+25A, are examined. For each source, we analysed quasi-simultaneous multiwavelength data for different states of γ-ray activity and present the results of their broad-band emission modelling, taking into account all available physical constraints to limit the range of the model parameters. Two different scenarios are discussed, in the framework of a one-zone leptonic model, where the high-energy emission is due to the inverse Compton scattering of the disc and broad line region (BLR) or torus photons by relativistic electrons within the jet. The transition from low to high state is well described by variations of the jet parameters, leaving the external photon fields unchanged. The parameterization favours an emission scenario with particle injection on a stationary shock inside the jet. When considering all physical constraints, the disc and BLR scenario is preferred for all three sources. We use the multi-epoch modelling to characterize total jet powers and discuss the intrinsic nature of γ-NLS1 galaxies and FSRQs.

Probing gamma-ray burst afterglows with the Cherenkov Telescope Array

ABSTRACT Detection of delayed sub-TeV photons from gamma-ray bursts (GRBs) by MAGIC and HESS has proven the promising future of GRB afterglow studies with the Cherenkov Telescope Array (CTA), the next-generation gamma-ray observatory. With the unprecedented sensitivity of CTA, afterglow detection rates are expected to increase dramatically. In this paper, we explore the multidimensional afterglow parameter space to see the detectability of sub-TeV photons by CTA. We use a one-zone electron synchrotron and synchrotron self-Compton model to obtain the spectral energy distribution. We consider bursts going off in a medium of homogenous density. The blast wave is assumed to be radiatively inefficient and evolving adiabatically. Considering that the electron acceleration is not efficient if the acceleration time-scale exceeds the radiative cooling time-scale, we find that the sub-TeV emission is always due to the self-Compton process. We find that jets with high kinetic energy or large bulk Lorentz factor decelerating into a dense ambient medium offer better detection prospects for CTA. For relatively lower values of the downstream magnetic field, electrons are slow-cooling, and the emitted radiation is positively correlated with the magnetic field. For larger magnetic fields, the electron population enters the fast-cooling phase where the radiated flux is inversely proportional to the magnetic field. We apply our results in the context of bright TeV afterglows detected in recent years. Our results indicate that cosmological short GRBs have only moderate prospects of detection by CTA while local neutron star merger counterparts can be detected if the jet is launched towards the observer.

Toward nebular spectral modeling of magnetar-powered supernovae

Context. Many energetic supernovae (SNe) are thought to be powered by the rotational energy of a highly magnetized, rapidly rotating neutron star. The emission from the associated luminous pulsar wind nebula (PWN) can photoionize the SN ejecta, leading to a nebular spectrum of the ejecta with signatures that might reveal the PWN. SN 2012au is hypothesized to be one such SN. Aims. We investigate the impact of different ejecta and PWN parameters on the SN nebular spectrum, and test whether any photoionization models are consistent with SN 2012au. We study how constraints from the nebular phase can be linked into modeling of the diffusion phase and the radio emission of the magnetar. Methods. We present a suite of late-time (1–6 yr) spectral simulations of SN ejecta powered by an inner PWN. Over a large grid of one-zone models, we study the behavior of the physical state and line emission of the SN as the PWN luminosity (LPWN), the injected spectral energy distribution (SED) temperature (TPWN), the ejecta mass (Mej), and the composition (pure O or realistic) vary. We discuss the resulting emission in the context of the observed behavior of SN 2012au, a strong candidate for a PWN-powered SN. We used optical light-curve models and broadband PWN models to predict possible radio emission from SN 2012au. Results. The SN nebular spectrum varies as TPWN varies because the ejecta become less ionized as TPWN increases. Ejecta models with low mass and high PWN power obtain runaway ionization for O I, and in extreme cases, also O II, causing a sharp decrease in their ion fraction over a small change in the parameter space. Certain models can reproduce the oxygen line luminosities of SN 2012au reasonably well at individual epochs, but we find no model that fits over the whole time evolution. This is likely due to uncertainties and simplifications in the model setup. Using our derived constraints from the nebular phase, we predict that the magnetar powering SN 2012au had an initial rotation period ~15 ms, and it is expected to be a strong radio source (F > 100 μJy) for decades.

General Physical Properties of Gamma-Ray-emitting Radio Galaxies

We study the radio galaxies with known redshifts detected by the Fermi satellite after 10 yr of data (4FGL-DR2). We use a one-zone leptonic model to fit the quasi-simultaneous multiwavelength data of these radio galaxies and study the distributions of the derived physical parameters as a function of black hole mass and accretion disk luminosity. The main results are as follows. (1) We find that the jet kinetic power of most radio galaxies can be explained by the hybrid jet model based on ADAFs surrounding Kerr black holes. (2) After excluding the redshift, there is a significant correlation between the radiation jet power and the accretion disk luminosity, while the jet kinetic power is weakly correlated with the accretion disk luminosity. (3) We also find a significant correlation between inverse Compton luminosity and synchrotron luminosity. The slope of the correlation for radio galaxies is consistent with the synchrotron self-Compton (SSC) process. The result may suggest that the high-energy components of radio galaxies are dominated by the SSC process.

Multiwavelength monitoring of the nucleus in PBC J2333.9−2343: the giant radio galaxy with a blazar-like core

ABSTRACT PBC J2333.9−2343 is a giant radio galaxy at z = 0.047 with a bright central core associated to a blazar nucleus. If the nuclear blazar jet is a new phase of the jet activity, then the small orientation angle suggests a dramatic change of the jet direction. We present observations obtained between 2018 September and 2019 January (cadence larger than three days) with Effeslberg, SMARTS-1.3m, ZTF, ATLAS, Swift, and Fermi-LAT, and between 2019 April and 2019 July (daily cadence) with SMARTS-1.3 m and ATLAS. Large (>2 ×) flux increases are observed on time-scales shorter than a month, which are interpreted as flaring events. The cross correlation between the SMARTS-1.3 m monitoring in the NIR and optical shows that these data do not show significant time lag within the measured errors. A comparison of the optical variability properties between non-blazars and blazars AGN shows that PBC J2333.9−2343 has properties more comparable to the latter. The SED of the nucleus shows two peaks, that were fitted with a one-zone leptonic model. Our data and modelling show that the high energy peak is dominated by External Compton from the dusty torus with mild contribution from Inverse Compton from the jet. The derived jet angle of 3 deg is also typical of a blazar. Therefore, we confirm the presence of a blazar-like core in the centre of this giant radio galaxy, likely a Flat Spectrum Radio Quasar with peculiar properties.

Broadband X-Ray Spectroscopy of the Pulsar Wind Nebula in HESS J1640-465

We present updated measurements of the X-ray properties of the pulsar wind nebula associated with the TeV γ-ray source HESS J1640-465 derived from Chandra and Nuclear Spectroscopic Telescope Array data. We report a high N H value along the line of sight, consistent with previous work, which led us to incorporate the effects of dust scattering in our spectral analysis. Due to uncertainties in the dust scattering, we report a range of values for the PWN properties (photon index and unabsorbed flux). In addition, we fit the broadband spectrum of this source and found evidence for spectral softening and decreasing unasborbed flux as we go to higher photon energies. We then used a one-zone time-dependent evolutionary model to reproduce the dynamical and multiwavelength spectral properties of our source. Our model suggests a short spin-down timescale, a relatively higher than average magnetized pulsar wind, a strong pulsar wind nebula magnetic field and maximum electron energy up to PeV, suggesting HESS J1640-465 could be a PeVatron candidate.

Offsets between X-Ray and Radio Components in X-Ray Jets: The AtlasX

The X-ray emission mechanism of powerful extragalactic jets—which has important implications for their environmental impacts—is poorly understood. The X-ray/radio positional offsets in the individual features of jets provide important clues. Extending previous work in Reddy et al., we present a detailed comparison between X-ray maps, deconvolved using the Low-count Image Reconstruction and Analysis tool, and radio maps of 164 components from 77 Chandra-detected X-ray jets. We detect 94 offsets (57%), with 58 new detections. In FR II–type jet knots, the X-rays peak and decay before the radio in about half the cases, disagreeing with the predictions of one-zone models. While a similar number of knots lack statistically significant offsets, we argue that projection and distance effects result in offsets below the detection level. Similar deprojected offsets imply that X-rays could be more compact than radio for most knots, and we qualitatively reproduce this finding with a “moving-knot” model. The bulk Lorentz factor (Γ) derived for knots under this model is consistent with previous radio-based estimates, suggesting that kiloparsec-scale jets are only mildly relativistic. An analysis of the X-ray/radio flux ratio distributions does not support the commonly invoked mechanism of X-ray production from inverse Compton scattering of the cosmic microwave background, but does show a marginally significant trend of declining flux ratio as a function of the distance from the core. Our results imply the need for multi-zone models to explain the X-ray emission from powerful jets. We provide an interactive list of our X-ray jet sample at http://astro.umbc.edu/Atlas-X.

Disentangling the nature of the prototype radio weak BL Lac

Context. The γ-ray emitting source WISE J141046.00+740511.2 has been associated with a Fermi-LAT detection by crossmatching with Swift/XRT data. It has shown all the canonical observational characteristics of a BL Lac source, including a power-law, featureless optical spectrum. However, it was only recently detected at radio frequencies and its radio flux is significantly low. Aims. Given that a radio detection is fundamental to associate lower-energy counterparts to Fermi-LAT sources, we aim to unambiguously classify this source by performing a multiwavelength analysis based on contemporaneous data. Methods. By using multifrequency observations at the Jansky Very Large Array, Giant Metrewave Radio Telescope, Gran Telescopio Canarias, Gemini, William Herschel Telescope and Liverpool observatories, together with Fermi-LAT and Swift data, we carried out two kinds of analyses. On one hand, we studied several known parameters that account for the radio loudness or weakness characterization and their application to blazars (in general) and to our source (in particular). And, on the other hand, we built and analyzed the observed spectral energy distribution (SED) of this source to try to explain its peculiar characteristics. Results. The multiwavelength analysis indicates that WISE J141046.00+740511.2 is a blazar of the high-frequency peaked (HBL) type that emits highly polarized light and that is likely located at a low redshift. In addition, the one-zone model parameters that best fit its SED are those of an extreme HBL (EHBL); this blazar type has been extensively predicted in theory to be lacking in the radio emission that is otherwise typical of canonical γ-ray blazars. Conclusions. We confirm that WISE J141046.00+740511.2 is indeed a highly polarized BL Lac of the HBL type. Further studies will be conducted to explain the atypical low radio flux detected for this source.

Optical/γ-ray blazar flare correlations: understanding the high-energy emission process using ASAS-SN and Fermi light curves

ABSTRACT Using blazar light curves from the optical All-Sky Automated Survey for Supernovae (ASAS-SN) and the γ-ray Fermi-LAT telescope, we performed the most extensive statistical correlation study between both bands, using a sample of 1180 blazars. This is almost an order of magnitude larger than other recent studies. Blazars represent more than 98 per cent of the AGNs detected by Fermi-LAT and are the brightest γ-ray sources in the extragalactic sky. They are essential for studying the physical properties of astrophysical jets from central black holes. However, their γ-ray flare mechanism is not fully understood. Multiwavelength correlations help constrain the dominant mechanisms of blazar variability. We search for temporal relationships between optical and γ-ray bands. Using a Bayesian Block Decomposition, we detect 1414 optical and 510 γ-ray flares, we find a strong correlation between both bands. Among all the flares, we find 321 correlated flares from 133 blazars, and derive an average rest-frame time delay of only 1.1$_{-8.5}^{+7.1}$ d, with no difference between the flat-spectrum radio quasars, BL Lacertae-like objects or low, intermediate, and high-synchrotron peaked blazar classes. Our time-delay limit rules out the hadronic proton-synchrotron model as the driver for non-orphan flares and suggests a leptonic single-zone model. Limiting our search to well-defined light curves and removing 976 potential but unclear ‘orphan’ flares, we find 191 (13 per cent) and 115 (22 per cent) clear ‘orphan’ optical and γ-ray flares. The presence of ‘orphan’ flares in both bands challenges the standard one-zone blazar flare leptonic model and suggests multizone synchrotron sites or a hadronic model for some blazars.

Investigating the radiative properties of LHAASO J1908 + 0621

ABSTRACT LHAASO J1908+0621 has recently been detected as a source emitting γ-rays with energies above 100 TeV, and multiband observations show that a break around 1 TeV appears in the γ-ray spectrum. We have reanalysed the GeV γ-ray properties of the 100-TeV source using 14 years of data recorded by the Fermi Large Area Telescope (Fermi-LAT). The spectrum in the energy range range 30–500 GeV has an index of 1.50 ± 0.26, which is much smaller than that detected in TeV γ-rays. Additionally, the radiation properties of this source are investigated based on a one-zone time-dependent model. In the model, LHAASO J1908+0621 is associated with a pulsar wind nebula (PWN) powered by the pulsar PSR J1907 + 0602. High-energy particles composed of electrons and positrons are injected into the nebula. Multiband non-thermal emission is produced via synchrotron radiation and inverse Compton scattering (ICS). Taking the effect of radiative energy losses and adiabatic cooling into account, the spectral energy distribution from a model with a broken power law for the distribution of the injected particles can explain the fluxes detected in the γ-ray bands. The results support the idea that LHAASO J1908 + 0621 originates from the PWN powered by PSR J1907 + 0602, and γ-rays with energy above 100 TeV are produced by electrons/positrons in the nebula via ICS.

A multimessenger study of the blazar PKS 0735+178: a new major neutrino source candidate

ABSTRACT The blazar PKS 0735+178 is possibly associated with multiple neutrino events observed by the IceCube, Baikal, Baksan, and KM3NeT neutrino telescopes while it was flaring in the γ-ray, X-ray, ultraviolet, and optical bands. We present a detailed study of this peculiar blazar to investigate the temporal and spectral changes in the multiwavelength emission when the neutrino events were observed. The analysis of Swift-XRT snapshots reveal a flux variability of more than a factor 2 in about 5 × 103 s during the observation on 2021 December 17. In the γ-ray band, the source was in its historical highest flux level at the time of the arrival of the neutrinos. The observational comparison between PKS 0735+178 and other neutrino source candidates, such as TXS 0506+056, PKS 1424+240, and GB6 J1542+6129, shows that all these sources share similar spectral energy distributions, very high radio and γ-ray powers, and parsec scale jet properties. Moreover, we present strong supporting evidence for PKS 0735+178 to be, like all the others, a masquerading BL Lac. We perform comprehensive modelling of the multiwavelength emission from PKS 0735+178 within one-zone lepto-hadronic models considering both internal and external photon fields and estimate the expected accompanying neutrino flux. The most optimistic scenario invokes a jet with luminosity close to the Eddington value and the interactions of ∼ PeV protons with an external UV photon field. This scenario predicts ∼0.067 muon and anti-muon neutrinos over the observed 3-week flare. Our results are consistent with the detection of one very high-energy neutrino like IceCube-211208A.

Chemical evolution of elliptical galaxies I: supernovae and AGN feedback

ABSTRACT We study the formation and evolution of elliptical galaxies and how they suppress star formation and maintain it quenched. A one-zone chemical model which follows in detail the time evolution of gas mass and its chemical abundances during the active and passive evolution is adopted. The model includes both gas infall and outflow as well as detailed stellar nucleosynthesis. Elliptical galaxies with different infall masses, following a down-sizing in star formation scenario, are considered. In the chemical evolution simulation, we include a novel calculation of the feedback processes. We include heating by stellar wind, core-collapse supernovae (SNe), Type Ia SNe (usually not highlighted in galaxy formation simulations), and active galactic nucleus (AGN) feedback. The AGN feedback is a novelty in this kind of models and is computed by considering a Bondi-Eddington limited accretion onto the central supermassive black hole. We successfully reproduce several observational features, such as the [α/Fe] ratios increasing with galaxy mass, mass-metallicity, MBH–σ and MBH–M* relations. Moreover, we show that stellar feedback and in particular Type Ia SNe, has a main role in maintaining quenched the star formation after the occurrence of the main galactic wind, especially in low-mass ellipticals. For larger systems, the contribution from AGN to thermal energy of gas appears to be necessary. However, the effect of the AGN on the development of the main galactic wind is negligible, unless an unreasonable high-AGN efficiency or an extremely low-stellar feedback are assumed. We emphasize the important role played by Type Ia SNe in the energy budget of early-type galaxies.

Globular cluster formation with multiple stellar populations: a single-binary composite scenario

ABSTRACT We discuss a GC formation scenario in which the first generation (1G) of single asymptotic giant branch (AGB) stars and intermediate-mass close binaries (IMCBs) eject gas, from which the second generation (2G) of stars can be formed. The two key parameters in the scenario are the fractions of binary stars (fb) and the slopes (α) of the stellar initial mass functions (IMFs) for 1G stars. Principle results derived by analytic and one-zone models of GC formation are as follows. The mass fraction of 2G stars (f2g) can be higher than ≈0.4 for α < 1.8 and is not so dependent on fb. The ratio of the initial mass of a GC to the present-day mass (Mgc) ranges from 2 to 7 depending on α for 0.5 ≤ fb ≤ 0.9. The differences in [Na/Fe] between 1G and 2G stars can be as large as 0.7 for a wide range of model parameters. The Li abundances of 2G stars can be as high as those of 1G even if the pristine gas from IMCBs is assumed to be Li-free. Formation histories of 2G stars show at least two peaks owing to two peaks in the total ejection rate of gas from IMCB populations. The observed correlation between f2g and Mgc can be due to α depending on Mgc. The hypothetical long duration of 2G formation (≈108 yr) is possible, because massive star formation can be suppressed through frequent dynamical interaction between 1G stars and gas clouds.