Context: White dwarf - Brown dwarf short period binaries (P$_{\rm orb}$ $\lesssim$ 2 hours) are some of the most extreme irradiated atmospheric environments known. These systems offer an opportunity to explore theoretical and modelling efforts of irradiated atmospheres different to typical hot Jupiter systems. Aims: We aim to investigate the three dimensional atmospheric structural and dynamical properties of the Brown dwarf WD0137-349B. Methods: We use the three dimensional GCM model Exo-FMS, with a dual-band grey radiative-transfer scheme to model the atmosphere of WD0137-349B. The results of the GCM model are post-processed using the three dimensional Monte Carlo radiative-transfer model \textsc{cmcrt}. Results: Our results suggest inefficient day-night energy transport and a large day-night temperature contrast for WD0137-349B. Multiple flow patterns are present, shifting energy asymmetrically eastward or westward depending on their zonal direction and latitude. Regions of overturning are produced on the western terminator. We are able to reproduce the start of the system near-IR emission excess at $\gtrsim$ 1.95 $\mu$m as observed by the GNIRS instrument. Our model over predicts the IR phase curve fluxes by factors of $\approx$1-3, but generally fits the shape of the phase curves well. Chemical kinetic modelling using \textsc{vulcan} suggests a highly ionised region at high altitudes can form on the dayside of the Brown dwarf. Conclusions: We present a first attempt at simulating the atmosphere of a short period White dwarf - Brown dwarf binary in a 3D setting. Further studies into the radiative and photochemical heating from the UV irradiation is required to more accurately capture the energy balance inside the Brown dwarf atmosphere. Cloud formation may also play an important role in shaping the emission spectra of the Brown dwarf.

Observations of brown dwarfs (BDs), free-floating planetary-mass objects, and directly imaged extrasolar giant planets (EGPs) exhibit rich evidence of large-scale weather. Understanding the mechanisms driving the vigorous atmospheric circulation of BDs and directly imaged EGPs and its effects on their observed lightcurve variability and spectral properties is a pressing need. Our previous work has suggested a strong role of radiative cloud feedback on inducing a spontaneous time evolution in a simple one-dimensional framework. Yet the radiative cloud feedback in a three-dimensional (3D) dynamical framework remains unexplored for conditions relevant to these objects. Here we present a series of atmospheric circulation models that self-consistently couple dynamics with idealized cloud formation and its radiative effects. We demonstrate that vigorous atmospheric circulation can be triggered and self-maintained by cloud radiative feedback. The circulation is dominated by cloud-forming and clear-sky vortices that evolve over timescales from several to tens of hours. The typical horizontal lengthscales of dominant vortices are strongly constrained by the rotation, showing a linear dependence on the inverse of rotation rate with stronger rotation leading to thinner clouds. Domain-mean outgoing radiative flux exhibits variability over timescales of tens of hours due to the statistical evolution of storms. The circulation driven by cloud radiative feedback represents a natural mechanism generating significant surface inhomogeneity as well as irregular flux variability. Our results also have important implications for near-IR colors of dusty BDs and EGPs, including the scatter in the near-IR color-magnitude diagram and the viewing-geometry dependent near-IR colors.

In this paper, we estimate the cross-correlation power spectra between the Planck 2018 cosmic microwave background (CMB) temperature anisotropy map and the unresolved $\gamma$-ray background (UGRB) from the 9-years Fermi-Large Area Telescope (LAT) data. In this analysis, we use up to nine energy bins over a wide energy range of [0.631, 1000] GeV from the Fermi-LAT UGRB data. Firstly, we find that the Fermi data with the energy ranges [1.202, 2.290] GeV and [17.38, 36.31] GeV show the positive evidence for the Integrated Sachs-Wolfe (ISW) effect at $1.8\sigma$ confidence level, and the significance would be increased to $2.7\sigma$ when using these two energy bins together. Secondly, we apply the single power-law model to normalize the amplitude and use all the nine Fermi energy bins to measure the significance of the ISW effect, we obtained ${\rm A_{amp}}=0.95 \pm 0.53$ ($68\%$ C.L.). For the robustness test, we implement a null hypothesis by randomizing the Fermi mock maps of nine energy bins and obtain the non-detection of ISW effect, which confirms that the ISW signal comes from the Fermi-LAT diffuse $\gamma$-ray data and is consistent with the standard $\Lambda$CDM model prediction essentially. We use a cross-correlation coefficient to show the relation between different energy bins. Furthermore, we vary the cut ranges $|b|$ of galactic plane on the mask of Fermi map and carefully check the consequent influence on the ISW signal detection.

Tidally locked gas giants are typically in several-day orbits, implying a modest role for rotation in the atmospheric circulation. Nevertheless, there exist a class of gas-giant, highly irradiated objects---brown dwarfs orbiting white dwarfs in extremely tight orbits---whose orbital and hence rotation periods are as short as 1-2 hours. Phase curves and other observations have already been obtained for this class of objects, raising fundamental questions about the role of increasing planetary rotation rate in controlling the circulation. So far, most modeling studies have investigated rotation periods exceeding a day, as appropriate for typical hot Jupiters. Here we investigate atmospheric circulation of tidally locked atmospheres with decreasing rotation periods down to 2.5 hours. With decreasing rotation period, the width of the equatorial eastward jet decreases, consistent with the narrowing of the equatorial waveguide due to a decrease of the equatorial deformation radius. The eastward-shifted equatorial hot spot offset decreases accordingly, and the off-equatorial westward-shifted hot areas become increasingly distinctive. At high latitudes, winds become weaker and more rotationally dominated. The day-night temperature contrast becomes larger due to the stronger influence of rotation. Our simulated atmospheres exhibit variability, presumably caused by instabilities and wave interactions. Unlike typical hot Jupiter models, thermal phase curves of rapidly rotating models show a near alignment of peak flux to secondary eclipse. This result helps to explain why, unlike hot Jupiters, many brown dwarfs orbiting white dwarfs exhibit IR flux peaks aligned with secondary eclipse. Our results have important implications for understanding fast-rotating, tidally locked atmospheres.

We report the identification of a positive cross-correlation signal between the unresolved $\gamma$-ray emission, measured by the \emph{Fermi} Large Area Telescope, and four different galaxy cluster catalogues. The selected catalogues peak at low-redshift and span different frequency bands, including infrared, optical and X-rays. The signal-to-noise ratio of the detected cross-correlation amounts to 3.5 in the most significant case. We investigate and comment about its possible origin, in terms of compact $\gamma$-ray emission from AGNs inside clusters or diffuse emission from the intracluster medium. The analysis has been performed by introducing an accurate estimation of the cross-correlation power-spectrum covariance matrix, built with mock realisations of the gamma and galaxy cluster maps. Different methods to produce the mock realizations starting from the data maps have been investigated and compared, identifying suitable techniques which can be generalized to other cross-correlation studies.

Phase curve observations provide an opportunity to study the full energy budgets of exoplanets by quantifying the amount of heat redistributed from their daysides to their nightsides. Theories explaining the properties of phase curves for hot Jupiters have focused on the balance between radiation and dynamics as the primary parameter controlling heat redistribution. However, recent phase curves have shown deviations from the trends that emerge from this theory, which has led to work on additional processes that may affect hot Jupiter energy budgets. One such process, molecular hydrogen dissociation and recombination, can enhance energy redistribution on ultra-hot Jupiters with temperatures above $\sim2000$ K. In order to study the impact of H$_{2}$ dissociation on ultra-hot Jupiters, we present a phase curve of KELT-9b observed with the Spitzer Space Telescope at 4.5 $\mu$m. KELT-9b is the hottest known transiting planet, with a 4.5-$\mu$m dayside brightness temperature of $4566^{+140}_{-136}$ K and a nightside temperature of $2556^{+101}_{-97}$ K. We observe a phase curve amplitude of $0.609 \pm 0.020$ and a hot spot offset of $18.7^{+2.1}_{-2.3}$ degrees. The observed amplitude is too small to be explained by a simple balance between radiation and advection. General circulation models (GCMs) and an energy balance model that include the effects of H$_{2}$ dissociation and recombination provide a better match to the data. The GCMs, however, predict a maximum hot spot offset of $5$ degrees, which disagrees with our observations at $>5\sigma$ confidence. This discrepancy may be due to magnetic effects in the planet's highly ionized atmosphere.

Recent observations of ultra-hot Jupiters with dayside temperatures in excess of $2500~\mathrm{K}$ have found evidence for new physical processes at play in their atmospheres. In this work, we investigate the effects of the dissociation of molecular hydrogen and recombination of atomic hydrogen on the atmospheric circulation of ultra-hot Jupiters. To do so, we incorporate these effects into a general circulation model (GCM) for hot Jupiter atmospheres, and run a large suite of models varying the incident stellar flux, rotation period, and strength of frictional drag. We find that including hydrogen dissociation and recombination reduces the fractional day-to-night temperature contrast of ultra-hot Jupiter atmospheres and causes the speed of the equatorial jet to decrease in simulations with fixed rotation. This is because the large energy input required for hydrogen dissociation cools the dayside of the planet, and the energy released due to hydrogen recombination warms the nightside. The resulting decrease in the day-to-night temperature contrast reduces the day-to-night pressure gradient that drives the circulation, resulting in weaker wind speeds. The results from our GCM experiments qualitatively agree with previous theory which found that the fractional day-night temperature contrast of ultra-hot Jupiters should decrease with increasing equilibrium temperature due to hydrogen dissociation and recombination. Lastly, we compute full-phase light curves from our suite of GCM experiments, finding that the reduced day-to-night temperature contrast in ultra-hot Jupiter atmospheres causes a smaller phase curve amplitude. The reduction in phase curve amplitude due to hydrogen dissociation and recombination could explain the relatively small phase curve amplitudes of observed ultra-hot Jupiters.

Brown dwarfs and directly imaged giant planets exhibit significant evidence for active atmospheric circulation, which induces a large-scale patchiness in the cloud structure that evolves significantly over time, as evidenced by infrared light curves and Doppler maps. These observations raise critical questions about the fundamental nature of the circulation, its time variability, and the overall relationship to the circulation on Jupiter and Saturn. Jupiter and Saturn themselves exhibit numerous robust zonal (east-west) jet streams at the cloud level; moreover, both planets exhibit long-term stratospheric oscillations involving perturbations of zonal wind and temperature that propagate downward over time on timescales of ~4 years (Jupiter) and ~15 years (Saturn). These oscillations, dubbed the Quasi Quadrennial Oscillation (QQO) for Jupiter and the Semi-Annual Oscillation (SAO) on Saturn, are thought to be analogous to the Quasi-Biennial Oscillation (QBO) on Earth, which is driven by upward propagation of equatorial waves from the troposphere. To investigate these issues, we here present global, three-dimensional, high-resolution numerical simulations of the flow in the stratified atmosphere--overlying the convective interior--of brown dwarfs and Jupiter-like planets. The effect of interior convection is parameterized by inducing small-scale, randomly varying perturbations in the radiative-convective boundary at the base of the model. In the simulations, the convective perturbations generate atmospheric waves and turbulence that interact with the rotation to produce numerous zonal jets. Moreover, the equatorial stratosphere exhibits stacked eastward and westward jets that migrate downward over time, exactly as occurs in the terrestrial QBO, Jovian QQO, and Saturnian SAO. This is the first demonstration of a QBO-like phenomenon in 3D numerical simulations of a giant planet.

The Gamma-Ray Integrated Detectors (GRID) is a space mission concept dedicated to monitoring the transient gamma-ray sky in the energy range from 10 keV to 2 MeV using scintillation detectors onboard CubeSats in low Earth orbits. The primary targets of GRID are the gamma-ray bursts (GRBs) in the local universe. The scientific goal of GRID is, in synergy with ground-based gravitational wave (GW) detectors such as LIGO and VIRGO, to accumulate a sample of GRBs associated with the merger of two compact stars and study jets and related physics of those objects. It also involves observing and studying other gamma-ray transients such as long GRBs, soft gamma-ray repeaters, terrestrial gamma-ray flashes, and solar flares. With multiple CubeSats in various orbits, GRID is unaffected by the Earth occultation and serves as a full-time and all-sky monitor. Assuming a horizon of 200 Mpc for ground-based GW detectors, we expect to see a few associated GW-GRB events per year. With about 10 CubeSats in operation, GRID is capable of localizing a faint GRB like 170817A with a 90% error radius of about 10 degrees, through triangulation and flux modulation. GRID is proposed and developed by students, with considerable contribution from undergraduate students, and will remain operated as a student project in the future. The current GRID collaboration involves more than 20 institutes and keeps growing. On August 29th, the first GRID detector onboard a CubeSat was launched into a Sun-synchronous orbit and is currently under test.

The study of the cross-correlation angular power spectrum between gravitational tracers and electromagnetic signals can be a powerful tool to constrain Dark Matter (DM) microscopic properties. In this work we correlate \Fermi\ diffuse \g-ray maps with catalogues of galaxy clusters. To emphasize the sensitivity to a DM signal, we select clusters at low-redshift $0<z<0.2$ and with large-halo mass $M_{500}>10^{13}M_\odot$. The analysis is performed with four catalogues in different wavebands, including infrared, optical and X-rays. No evidence for a DM signal is identified. On the other hand, we derive competitive bounds: the thermal cross-section is excluded at 95\% C.L. for DM masses below 20 GeV and annihilation in the $\tau^+-\tau^-$ channel.

Growing observational evidence has suggested active meteorology in atmospheres of brown dwarfs (BDs) and directly imaged extrasolar giant planets (EGPs). In particular, a number of surveys have shown that near-IR brightness variability is common among L and T dwarfs. Despite initial understandings of atmospheric dynamics which is the major cause of the variability by previous studies, the detailed mechanism of variability remains elusive, and we need to seek a natural, self-consistent mechanism. Clouds are important in shaping the thermal structure and spectral properties of these atmospheres via large opacity, and we expect the same for inducing atmospheric variability. In this work, using a time-dependent one-dimensional model that incorporates a self-consistent coupling between the thermal structure, convective mixing, cloud radiative heating/cooling and condensation/evaporation of clouds, we show that radiative cloud feedback can drive spontaneous atmospheric variability in both temperature and cloud structure in conditions appropriate for BDs and directly imaged EGPs. The typical periods of variability are one to tens of hours with typical amplitude of the variability up to hundreds of Kelvins in effective temperature. The existence of variability is robust over a wide range of parameter space, but the detailed evolution of variability is sensitive to model parameters. Our novel, self-consistent mechanism has important implications for the observed flux variability of BDs and directly imaged EGPs, especially those evolve in short timescales. It is also a promising mechanism for cloud breaking, which has been proposed to explain the L/T transition of brown dwarfs.

The L/T transition is an important evolutionary phase in brown dwarf atmospheres, providing us with a unique opportunity to explore the effects of clouds, convection, winds, gravity and metallicity across a very narrow temperature range. Understanding these physical processes is critical for understanding ultracool atmospheres. In the next decade, we will answer three key questions regarding L/T transition atmospheres: 1. What is the physical mechanism behind the L/T transition? 2. What is the spatial extent of atmospheric structures at the L/T transition? 3. How do gravity and metallicity affect the L/T transition? The theory and methods developed for brown dwarfs will be used in the 2030s and beyond for solar-system age giant exoplanets and eventually habitable zone earth analogues. Developing these techniques now are crucial.

The Probe of Inflation and Cosmic Origins (PICO) is a probe-class mission concept currently under study by NASA. PICO will probe the physics of the Big Bang and the energy scale of inflation, constrain the sum of neutrino masses, measure the growth of structures in the universe, and constrain its reionization history by making full sky maps of the cosmic microwave background with sensitivity 80 times higher than the Planck space mission. With bands at 21-799 GHz and arcmin resolution at the highest frequencies, PICO will make polarization maps of Galactic synchrotron and dust emission to observe the role of magnetic fields in Milky Way's evolution and star formation. We discuss PICO's optical system, focal plane, and give current best case noise estimates. The optical design is a two-reflector optimized open-Dragone design with a cold aperture stop. It gives a diffraction limited field of view (DLFOV) with throughput of 910 square cm sr at 21 GHz. The large 82 square degree DLFOV hosts 12,996 transition edge sensor bolometers distributed in 21 frequency bands and maintained at 0.1 K. We use focal plane technologies that are currently implemented on operating CMB instruments including three-color multi-chroic pixels and multiplexed readouts. To our knowledge, this is the first use of an open-Dragone design for mm-wave astrophysical observations, and the only monolithic CMB instrument to have such a broad frequency coverage. With current best case estimate polarization depth of 0.65 microK(CMB}-arcmin over the entire sky, PICO is the most sensitive CMB instrument designed to date.

Secondary eclipse observations of ultra-hot Jupiters have found evidence that hydrogen is dissociated on their daysides. Additionally, full-phase light curve observations of ultra-hot Jupiters show a smaller day-night emitted flux contrast than that expected from previous theory. Recently, it was proposed by Bell & Cowan (2018) that the heat intake to dissociate hydrogen and heat release due to recombination of dissociated hydrogen can affect the atmospheric circulation of ultra-hot Jupiters. In this work, we add cooling/heating due to dissociation/recombination into the analytic theory of Komacek & Showman (2016) and Zhang & Showman (2017) for the dayside-nightside temperature contrasts of hot Jupiters. We find that at high values of incident stellar flux, the day-night temperature contrast of ultra-hot Jupiters may decrease with increasing incident stellar flux due to dissociation/recombination, the opposite of that expected without including the effects of dissociation/recombination. We propose that a combination of a greater number of full-phase light curve observations of ultra-hot Jupiters and future General Circulation Models that include the effects of dissociation/recombination could determine in detail how the atmospheric circulation of ultra-hot Jupiters differs from that of cooler planets.

Growing observations of brown dwarfs have provided evidence for strong atmospheric circulation on these objects. Directly imaged planets share similar observations, and can be viewed as low-gravity versions of brown dwarfs. Vigorous condensate cycles of chemical species in their atmospheres are inferred by observations and theoretical studies, and latent heating associated with condensation is expected to be important in shaping atmospheric circulation and influencing cloud patchiness. We present a qualitative description of the mechanisms by which condensational latent heating influence the circulation, and then illustrate them using an idealized general circulation model that includes a condensation cycle of silicates with latent heating and molecular weight effect due to rainout of condensate. Simulations with conditions appropriate for typical T dwarfs exhibit the development of localized storms and east-west jets. The storms are spatially inhomogeneous, evolving on timescale of hours to days and extending vertically from the condensation level to the tropopause. The fractional area of the brown dwarf covered by active storms is small. Based on a simple analytic model, we quantitatively explain the area fraction of moist plumes, and show its dependence on radiative timescale and convective available potential energy. We predict that, if latent heating dominates cloud formation processes, the fractional coverage area by clouds decreases as the spectral type goes through the L/T transition from high to lower effective temperature. This is a natural consequence of the variation of radiative timescale and convective available potential energy with spectral type.

The full-phase infrared light curves of low-eccentricity hot Jupiters show a trend of increasing fractional dayside-nightside brightness temperature difference with increasing incident stellar flux, both averaged across the infrared and in each individual wavelength band. The analytic theory of Komacek & Showman (2016) shows that this trend is due to the decreasing ability with increasing incident stellar flux of waves to propagate from day to night and erase temperature differences. Here, we compare the predictions of this theory to observations, showing that it explains well the shape of the trend of increasing dayside-nightside temperature difference with increasing equilibrium temperature. Applied to individual planets, the theory matches well with observations at high equilibrium temperatures but, for a fixed photosphere pressure of $100 \ \mathrm{mbar}$, systematically under-predicts the dayside-nightside brightness temperature differences at equilibrium temperatures less than $2000 \ \mathrm{K}$. We interpret this as due to as the effects of a process that moves the infrared photospheres of these cooler hot Jupiters to lower pressures. We also utilize general circulation modeling with double-grey radiative transfer to explore how the circulation changes with equilibrium temperature and drag strengths. As expected from our theory, the dayside-nightside temperature differences from our numerical simulations increase with increasing incident stellar flux and drag strengths. We calculate model phase curves using our general circulation models, from which we compare the broadband infrared offset from the substellar point and dayside-nightside brightness temperature differences against observations, finding that strong drag or additional effects (e.g. clouds and/or supersolar metallicities) are necessary to explain many observed phase curves.

We report the discovery of a new planetary system around the K giant $\eta$ Cet (HIP 5364, HD 6805) based on 118 high-precision optical radial velocities taken at Lick Observatory since July 2000. Since October 2011 an additional nine near-infrared Doppler measurements have been taken using the ESO CRIRES spectrograph (VLT, UT1). The visible data set shows two clear periodicities. Although we cannot completely rule out that the shorter period is due to rotational modulation of stellar features, the infrared data show the same variations as in the optical, which strongly supports that the variations are caused by two planets. Assuming the mass of $\eta$ Cet to be 1.7 $M_\odot$, the best edge-on coplanar dynamical fit to the data is consistent with two massive planets ($m_b\sin i$ = 2.6 $\pm$ 0.2 $M_{\mathrm{Jup}}$, $m_c\sin i$ = 3.3 $\pm$ 0.2 $M_{\mathrm{Jup}}$), with periods of $P_b$ = 407 $\pm$ 3 days and $P_c$ = 740 $\pm$ 5 days and eccentricities of $e_b$ = 0.12 $\pm$ 0.05 and $e_c$ = 0.08 $\pm$ 0.03. We tested a wide variety of edge-on coplanar and inclined planetary configurations for stability, which agree with the derived radial velocities. We find that in certain coplanar orbital configurations with moderate $e_b$ eccentricity, the planets can be effectively trapped in an anti-aligned 2:1 mean motion resonance. A much larger non-resonant stable region exists in low-eccentricity parameter space, although it appears to be much farther from the best fit than the 2:1 resonant region. In all other cases, the system is categorized as unstable or chaotic. Another conclusion from the coplanar inclined dynamical test is that the planets can be at most a factor of $\sim$ 1.4 more massive than their suggested minimum masses. This stability constraint on the inclination excludes the possibility of two brown dwarfs, and strongly favors a planetary system.

The frequency ratios and of HD 49933 exhibit an increase at high frequencies. This behavior also exists in the ratios of other stars, which is considered to result from the low signal-to-noise ratio and the larger line width at the high-frequency end and could not be predicted by stellar models in previous work. Our calculations show that the behavior not only can be reproduced by stellar models, but can be predicted by asymptotic formulas of the ratios. The frequency ratios of the Sun, too, can be reproduced well by the asymptotic formulas. The increased behavior derives from the fact that the gradient of mean molecular weight at the bottom of the radiative region hinders the propagation of p-modes, while the hindrance does not exist in the convective core. This behavior should exist in the ratios of stars with a large convective core. The characteristic of the ratios at high frequencies provides a strict constraint on stellar models and aids in determining the size of the convective core and the extent of overshooting. Observational constraints point to a star with $M=1.28\pm0.01 M_{\odot}$, $R=1.458\pm0.005 R_{\odot}$, $t=1.83\pm0.1$ Gyr, $r_{cc}=0.16\pm0.02 R_{\odot}$, $\alpha=1.85\pm0.05$, and $\delta_{ov}=0.6\pm0.2$ for HD 49933.

We present six years of new radial-velocity data from the Anglo-Australian and Magellan Telescopes on the HD 73526 2:1 resonant planetary system. We investigate both Keplerian and dynamical (interacting) fits to these data, yielding four possible configurations for the system. The new data now show that both resonance angles are librating, with amplitudes of 40 degrees and 60 degrees, respectively. We then perform long-term dynamical stability tests to differentiate these solutions, which only differ significantly in the masses of the planets. We show that while there is no clearly preferred system inclination, the dynamical fit with i=90 degrees provides the best combination of goodness-of-fit and long-term dynamical stability.

We present an updated analysis of radial velocity data of the HD 82943 planetary system based on 10 years of measurements obtained with the Keck telescope. Previous studies have shown that the HD 82943 system has two planets that are likely in 2:1 mean-motion resonance (MMR), with the orbital periods about 220 and 440 days (Lee et al. 2006). However, alternative fits that are qualitatively different have also been suggested, with two planets in a 1:1 resonance (Gozdziewski & Konacki 2006) or three planets in a Laplace 4:2:1 resonance (Beauge et al. 2008). Here we use \c{hi}2 minimization combined with parameter grid search to investigate the orbital parameters and dynamical states of the qualitatively different types of fits, and we compare the results to those obtained with the differential evolution Markov chain Monte Carlo method. Our results support the coplanar 2:1 MMR configuration for the HD 82943 system, and show no evidence for either the 1:1 or 3-planet Laplace resonance fits. The inclination of the system with respect to the sky plane is well constrained at about 20(+4.9 -5.5) degree, and the system contains two planets with masses of about 4.78 MJ and 4.80 MJ (where MJ is the mass of Jupiter) and orbital periods of about 219 and 442 days for the inner and outer planet, respectively. The best fit is dynamically stable with both eccentricity-type resonant angles {\theta}1 and {\theta}2 librating around 0 degree.

A novel theoretical method is developed to study the polycyclic aromatic hydrocarbon - diffuse interstellar band (PAH-DIB) hypothesis. In this method, a computer program is used to enumerate all PAH molecules with up to a specific number of fused benzene rings. Fast quantum chemical calculations are then employed to calculate the electronic transition energies, oscillator strengths, and rotational constants of these molecules. An electronic database of all PAHs with up to any specific number of benzene rings can be constructed this way. Comparison of the electronic transition energies, oscillator strengths, and rotational band contours of all PAHs in the database with astronomical spectra allows one to identify possible individual PAH carriers of some of the intense narrow DIBs. Using the current database containing up to 10 benzene rings we have selected 8 closed-shell PAHs as possible carriers of the intense lambda6614 DIB.