On the origin of stars with and without planets. Tc trends and clues to Galactic evolution

Abstract: We explore a sample of 148 solar-like stars to search for a possible correlation between the slopes of the abundance trends versus condensation temperature (known as the Tc slope) with stellar parameters and Galactic orbital parameters in order to understand the nature of the peculiar chemical signatures of these stars and the possible connection with planet formation. We find that the Tc slope significantly correlates (at more than 4sigma) with the stellar age and the stellar surface gravity. We also find tentative evidence that the Tc slope correlates with the mean galactocentric distance of the stars (Rmean), suggesting that those stars that originated in the inner Galaxy have fewer refractory elements relative to the volatiles. While the average Tc slope for planet-hosting solar analogs is steeper than that of their counterparts without planets, this difference probably reflects the difference in their age and Rmean. We conclude that the age and probably the Galactic birth place are determinant to establish the star's chemical properties. Old stars (and stars with inner disk origin) have a lower refractory-to-volatile ratio.

\section{Introduction}

Despite the huge progress in developing instrumentation and observational techniques during the past decade,
the study of extrasolar planets' properties via direct observations is still a very difficult task, and
the precise study and characterization of known extroplanets cannot be dissociated from the study of planet host stars.

The connection between stellar and planetary properties has been widely explored. In particular, the very first correlation
observed in this field of research, that of giant-planet  - metallicity \citep[e.g.,][]{Gonzalez-97, Santos-01, Santos-04, Fischer-05, Sousa-11}, was one of
the most important constraints on planet formation theories \citep[e.g.,][]{Mordasini-09}. Afterwards it was shown that
not only does the presence of planets correlate with metallicity (usually abundance of iron), but planet-host stars also show a systematic
enhancement of other elements \citep[e.g., $\alpha$-elements - ][]{Haywood-08b, Adibekyan-12a, Adibekyan-12b}. More recently, several studies
have revealed that
stellar metallicity also plays an important role in the architecture of planetary orbits \citep[e.g.,][]{Dawson-13, Beauge-13, Adibekyan-13a}.

Naturally, this connection is bidirectional: not only do stellar properties play an important role in planet formation
and evolution, but the planet formation can also have an impact on stellar properties. As an example, \cite{Israelian-09} and \cite{DelgadoMena-13}
find strong evidence that solar analogs with planets show enhanced Li depletion compared to their non-host counterparts.
The studies mentioned above are not the only ones.  Recently, \cite{Melendez-09} have claimed that the Sun shows a deficiency in refractory
elements with respect to other solar twins probably because they were trapped in the terrestrial planets in our solar system.
The same conclusion has also been reached by \cite{Ramirez-09}  who analyzed 64 solar twins and analogs.  However, recent results by
\cite[][- GH10,13]{Jonay-10, Jonay-13} strongly challenge the relation between the presence of planets and the abundance peculiarities of the stars.
Other works also have examined this possible connection between the chemical peculiarities and formation of planets
\citep[e.g.,][]{Smith-01, Ecuvillon-06, Gonzalez-10, Gonzalez-11, Jonay-11, Schuler-11, Ramirez-14}.

In this letter, we explore the origin of the possible trend observed between [X/H] (or [X/Fe]) and condensation
temperature ($T_{c}$) using a sample of 148 solar-like stars.

\section{Data}

Our initial sample is a combination of two samples of solar analogs (95 stars) and ``hot'' analogs (61 stars) taken from
GH10,13. We have
cross-matched this sample with the Geneva-Copenhagen Survey sample \citep[GCS-][]{Nordstrom-04}, for which \cite{Casagrande-11}
provides the Galactic orbital parameters, the space velocity components, and the ages of 148 of the stars considered in our study%
\footnote{Throughout the paper, BASTI expectation ages are used as suggested by \cite{Casagrande-11}.%
}. Fifty-seven of these stars are planet hosts, while for the remaining 91 no planetary companion has been detected up to now.

The stellar atmospheric parameters and the slopes of the $\Delta$[X/Fe]$_{SUN-star}$ versus  T$_{c}$ were derived using very high-quality
HARPS spectra%
\footnote{Zero slope means solar chemical composition, and a positive slope corresponds to a smaller refractory-to-volatile ratio compared to the Sun.%
}. Twenty-five elements from C (Z = 6) to Eu (Z = 63) have been used for this analysis.
These slopes are corrected for the Galactic chemical evolution trends as discussed in GH10,13.

The stars in the sample have effective temperatures 5604 \emph{K} $\leq$ T$_{eff}$ $\leq$ 6374 \emph{K},
metallicites -0.29 $\leq$ {[}Fe/H{]} $\leq$ 0.38 dex, and surface gravities
4.14 $\leq$ $\log g$ $\leq$ 4.63 dex. Throughout the letter we defined solar analogs as stars with; T$_{eff}$ = 5777$\pm$200 K;
logg = 4.44$\pm$0.20 dex; [Fe/H] = 0.0$\pm$0.2 dex. Solar twins are defined as with; T$_{eff}$ = 5777$\pm$100 K;
logg = 4.44$\pm$0.10 dex; [Fe/H] = 0.0$\pm$0.1 dex. Fifteen out of 58 solar analogs in this sample are known to be orbited by planets and
three out of 15 solar twins are planet hosts.

We would like to point out that stellar parameters derived in GH10,13 and in \cite{Casagrande-11} are in perfect agreement
($\Delta$T$_{eff}$ = -8$\pm$50 K, $\Delta$ $\log g$ = 0.06$\pm$0.07 dex, and $\Delta$[Fe/H] = -0.02$\pm$0.06 dex), which means
that the ages derived in \cite{Casagrande-11} are indeed coherent with the rest of the observables.

\begin{table}
\centering
\caption{Significance levels of the obtained correlations.}
\label{ks}
\begin{tabular}{lcccc}
\hline
\hline
Parameter & Cor. coeff & Significance & Cor. coeff & Significance \tabularnewline
& \multicolumn{2}{c}{  All } & \multicolumn{2}{c}{Solar analogs}\tabularnewline
\hline
Age  & 0.35 & \textbf{4.37} & 0.60 & \textbf{4.57}\tabularnewline
$\log g$  & -0.50 & \textbf{6.01} & -0.73 & \textbf{5.57}\tabularnewline
$[$Fe/H$]$  & -0.14 & 1.73 & 0.00 & 0.05\tabularnewline
T$_{eff}$ & -0.05 & 0.59 & 0.03 & 0.17\tabularnewline

\hline
\end{tabular}
\end{table}

\section{Correlations with Tc slope}

We searched for possible correlations between the T$_{c}$ slope and, in turn, atmospheric parameters and also Galactic orbital parameters and age,
in order to understand which is/are the main factor(s) responsible for the abundance trends with T$_{c}$.

\subsection{ T$_{c}$ slope against stellar parameters and age}

In Fig.~\ref{fig_slope_logg} we show the relation between the T$_{c}$ slopes and $\log g$  both for the full sample and for the
solar analogs exclusively. One can note that the T$_{c}$ slopes strongly correlate with the $\log g$. To evaluate the significance of the correlation
we performed a simple bootstrapped Monte Carlo test. For more details about the test we refer the reader to \cite{Figueira-13}
and \cite{Adibekyan-13b}.
The correlation coefficients and the significance levels (z-scores) of the correlations are presented in Table 1. The same table
shows that T$_{c}$ slopes do not significantly correlate with other stellar parameters. The T$_{c}$ slope versus {[}Fe/H{]} and T$_{eff}$ plots
are shown in Fig.~\ref{fig_slope_atmos_param}.

The observed significant relation between the surface gravity and T$_{c}$ slope means that the chemical abundance trends depend
either on the evolutionary stage of the star or on its age%
\footnote{Since a solar type star evolves from ZAMS, its $\log g$  slowly increases with age along the main sequence.%
}. In Fig.~\ref{fig_slope_age} we plot the T$_{c}$ slope against
the stellar age. The plot and Table 1 clearly show that the correlation with age is quite significant and confirms the result obtained for
the surface gravity: old stars are more “depleted“ in refractory elements (lower refractory-to-volatile ratios)
than their younger counterparts.

For FGK dwarf stars in the main sequence one does not expect significant changes in their atmospheric chemical abundances with age. This means
that the observed correlation between the T$_{c}$ slope and age probably reflects the chemical evolution in the Galaxy.
We note that this is the simplest assumption we can make based on our limited current knowledge of stellar evolution,
and we caution the reader that might be other effects that could severely affect the composition of stars  as a function of age.
For example, in the solar wind, elements with low first ionization potentials (FIP) are about four times more abundant than in the
photosphere of the sun \citep[e.g.,][]{Geiss-98, Raymond-99}. However, this so-called FIP-effect, which will depend on age and
has significant effect on coronal abundances \citep{Wood-06}, is not expected to affect the photospheric abundances.

In this context it is worth mentioning the work of \cite{Melendez-09}, where the authors discussed several possible explanations
of the solar ``peculiar'' abundances when compared to the solar twins. Most of the discussed possible effects (e.g.,
supernova pollution, early dust separation, etc.) however are not expected to have a dependence on age and can indeed be
responsible for peculiarities in chemical composition for individual cases (stars).
In fact the median age of their comparison solar-twin stars was of 4.1 Gyr, a very similar to that of the Sun.
The authors also discussed the Galactic evolution effects, taking into account the possibility that the Sun may have
migrated from an inner Galactic orbit \citep{Wielen-96}, but ended up not considering it as a probable explanation.

\subsection{ T$_{c}$ slope and Galactic orbital parameters}

Several studies have shown that the mean of the apo- and pericentric distances (R$_{mean}$) are good indicators of the stellar birthplace
\citep[e.g.,][]{Grenon-87, Edvardsson-93, Nordstrom-99,
Rocha-Pinto-04, Bensby-13}. \cite{Haywood-08a} shows that orbital parameters (R$_{mean}$) of the metal-poor and metal-rich thin disk stars are
significantly different from those of the main thin-disk population, and an outer and inner galactic disk origin was suggested for them
both, respectively.
A word of caution should be added here, if the radial migration proposed by \citet{Sellwood-02}%
\footnote{The term ``churning'' was introduced in \cite{Schonrich-09} for this type of migration.%
} is efficient, then using R$_{mean}$ as a proxy for the birthplace of a star could be dubious. However, \cite{Haywood-13} cast doubt upon
the efficiency of churning for contaminating the solar neighborhood.

In Fig.~\ref{fig_slope_rmean} we study the dependence of the T$_{c}$ slopes on R$_{mean}$, which we use as a proxy of the birth radii.
Even if the correlation between the T$_{c}$ slopes and R$_{mean}$ is not strong and is not significant (at 1-$\sigma$ level), one can see
that most of the stars with the lowest mean galactocentric distances have steeper average T$_{c}$ slopes ($\approx0.071\pm0.065$) than the average of the stars
with R$_{mean}$ = 8$\pm$1 kpc (T$_{c}$ $\approx-0.086\pm0.035$). The significance of this difference is 2.1$\sigma$ if one applies two sample \textit{t}-test.
For the solar analog sample, the significance of this difference is less significant (at a level of 1.5$\sigma$).

In Fig.~\ref{fig_slope_twins} we show the relation between the T$_{c}$ slopes and both the stellar age and the R$_{mean}$ for
solar twins in our sample. As one can see, all the correlations obtained for the full sample and sample of solar analogs are also valid for
the solar twins. Solar twins follow the general trends discussed in this section and in Sect. 3.1.

\begin{figure}
\begin{center}
\begin{tabular}{c}
\includegraphics[angle=0,width=0.8\linewidth]{slope_age_basti.ps}
\end{tabular}
\end{center}
\vspace{-0.6cm}
\caption{T$_{c}$ slopes versus ages for the full sample (\textit{top}) and for the solar analogs (\textit{bottom}).
Gray solid lines provide linear fits to the data points.}
\label{fig_slope_age}
\end{figure}

\cite{Rocha-Pinto-04}, among other authors, have already shown that R$_{mean}$ correlates with age%
\footnote{\cite{Wielen-96} have shown that this relation  can only be considered in a statistical sense.%
}in such a way that young objects (mostly younger than one Gyr)  all have R$_{mean}$ $\approx$ R$_{\odot}$, while older stars present a
higher proportion of objects coming from different galactocentric radii.
In Fig.~\ref{fig_rmean_age} we plot R$_{mean}$ versus stellar age. One can note that there is a weak correlation between the two parameters
(with a correlation coefficient of $\approx$ 0.16 at a significance level of  $\approx$ 2$\sigma$),
even if not as an apparent way as in Fig. 7 of \cite{Rocha-Pinto-04}, probably because we do not have very young stars in our sample.
However, in Fig.~\ref{fig_rmean_age} one can note that old stars (older than 6-7 Gyr) from the inner Galaxy ( R$_{mean}$ $<$ 7)
have predominantly positive T$_{c}$ slopes, while those in the solar circle have negative slopes.
In addition, the same figure shows that young stars with inner R$_{mean}$ values show an overabundance of
negative T$_{c}$ slopes, while their older counterparts mostly have positive slopes.

To test our findings observationally, i.e., if the T$_{c}$ slopes are different at different galactocentric distances (\textit{R}),
we used the Galactic abundance gradients derived in \cite{Lemasle-08, Lemasle-13}.
Using the linear fits (provided by these authors)
of the gradient for each elements, we derived [X/Fe] abundance ratios at four galactocentric radii (6, 8, and 10 kpc). Then
we plotted [X/Fe] against T$_{c}$ (see Fig.~\ref{fig_lemasle}), and after performing a linear fit, we calculated the T$_{c}$ slopes.
From the figure it is evident that at smaller galactocentric distances, the T$_{c}$ slope is steeper%
\footnote{We note that here one should changed the sign of the T$_{c}$ slopes to be consistent with the T$_{c}$ slopes derived in
GH10,13.%
} which qualitatively (not quantitatively) agrees well with what we obtained in the T$_{c}$ versus R$_{mean}$ plot.
We note that in  \cite{Lemasle-08, Lemasle-13}, the authors used young Galactic Cepheids to derive the gradients. This means that
the gradients should be considered as current ones, which in turn means that they are not affected by any type of migration.
Although not quantitative, this comparison shows that at a fixed age (the present day) the
T$_{c}$ correlates with the galactocentric radius.
This is one more hint that R$_{mean}$ could be a relevant parameter for "single-age" populations
(as already noted for the other end of the age domain explored in the paper, i.e. old stars, discussed in Fig.~\ref{fig_rmean_age}).

The tentative correlation between T$_{c}$ and the mean galactocentric distance of the stars, and the observational
results of \cite{Lemasle-08}, suggests that the chemical
composition of the ``birth place'' is a factor partially responsible for the observed  T$_{c}$ trends.
The observed (weak) correlations suggest that at a fixed time the
steeper  T$_{c}$ slopes (lower refractory-to-volatile ratios) are associated to stars that probably
have originated in the inner Galaxy.  This result is in line with the recent chemical evolution models of the Milky Way invoking
stellar migration \citep[e.g.,][]{Minchev-13}. Older stars are exposed to perturbations and scattering for a longer period of time
and have also had more time to migrate,
therefore some fraction of them originated in the inner Galaxy, while almost all the young stars have birth radii within the solar
radius \citep[][]{Wang-13, Minchev-13}.

\begin{figure}
\begin{center}
\begin{tabular}{c}
\includegraphics[angle=0,width=0.8\linewidth]{slope_rmean.ps}
\end{tabular}
\end{center}
\vspace{-0.6cm}
\caption{T$_{c}$ slopes versus R$_{mean}$ for the full sample (\textit{top}) and for the solar analogs (\textit{bottom}).
Gray solid lines provide linear fits to the data points.}
\label{fig_slope_rmean}
\end{figure}

\subsection{ T$_{c}$ slope and planets}

In this letter we do not develop a long discussion about the signatures of planet formation on the abundance trends since it has already been
discussed  in detail by GH10,13. Using the same sample they have already shown that there is no statistically significant
difference in T$_{c}$ slopes for planet-hosting stars (in particular for rocky planet hosts) and stars without any detected planetary companion.
In fact the average T$_{c}$ slope of the planet hosting stars (-0.08$\pm$0.38) is even slightly smaller (the difference is not significant)
than that of the non-host stars (-0.05$\pm$0.39) from the full sample.

Following our definition of solar analogs, we found that the average of the T$_{c}$ slope for planet hosting solar analogs is
greater (0.012$\pm$0.31) than that of their non-host counterparts (-0.16$\pm$0.34). The Kolmogorov-Smirnov (K-S) statistics
predict the $\approx$ 0.21  probability ($P_{KS}$) that these two
subsamples came from the same underlying distribution for T$_{c}$ slope. At the same time, the same statistics predict a
$P_{KS}$ $\approx$ 0.20 probability that they stem from the same underlying age distributions.
The latter can be seen in Fig.~\ref{fig_slope_age}: most of the planet-hosting
stars tend to be relatively old ($\gtrsim$ 5 Gyr).
Moreover,  planet host and non-host samples show a different distribution of  R$_{mean}$ -- $P_{KS}$ $\approx$ 0.007.
As can also be seen in Fig. ~\ref{fig_slope_rmean}, 10 planet hosts out of 15 (66\%) have R$_{mean}$ smaller than 7.5 kpc (where
slopes are usually high) and only 16 out of 43 stars without detected planets (37\%) have similarly low R$_{mean}$ values.
Clearly the two subsamples are not consistent with respect to the mean galactocentric distance and age.
Interestingly, \cite{Haywood-09} has already shown that (giant) planet host stars tend to have smaller R$_{mean}$ and probably
originate in the inner disk, which follow the same direction as our findings.

We note that the mean value of the T$_{c}$ slopes of the super-Earth hosts (M$_{planet}$ $<$ 0.05 M$_{J}$) is very close to the
average value of the Jovian hosts (M$_{planet}$ $>$ 0.1 M$_{J}$), while the average T$_{c}$ slope of the three Neptunian hosts%
\footnote{The Neptunian hosts also host Jovian and super-Earth like planets.}%
(0.05 M$_{J}$ $\geq$  M$_{planet}$ $\geq$ 0.1 M$_{J}$) is greater. For a more detailed comparison  of the T$_{c}$ slopes for planet hosts
at different planetary mass regimes, we refer the reader to GH10,13.

These results suggest that the difference in T$_{c}$ slopes observed for solar analogs with and without planets is probably due to the
differences in their ``birth places'' and times.

\section{Summary and conclusion}

We used a sample of 148 solar-type stars from GH10,13 to explore the main factors responsible for the abundance
trends with condensation temperature. For these stars the stellar atmospheric parameters and the T$_{c}$ slopes were taken from the
above-mentioned studies, while the stellar ages,  Galactic orbital parameters, and velocity components are from \cite{Casagrande-11}.

Our study reveals a strong correlation between stellar ages and T$_{c}$ slope: old stars show steeper slope i.e., less
refractory elements relative to volatiles. The same result reflects the statistically significant correlation between T$_{c}$ slopes
and surface gravity: more evolved (old) stars have a lower refractory-to-volatile ratio.

Moving one step further, we found tentative evidence that the T$_{c}$ slopes also correlate  with the mean galactocentric distance
of the stars; this suggest that stars which probably have origin in the inner Galaxy (small R$_{mean}$) have steeper slopes.
The result fits well in the recent evolution picture of the Milky Way, showing that some fraction of old stars in
the solar neighborhood might have its origin in the inner disk \citep[e.g.,][]{Minchev-13}.

Briefly exploring the possible reasons why one can see a difference in T$_{c}$ slopes for planet-hosting solar analogs and solar analogs
without detected planets, we found that in the current sample these two subsamples have different distributions of age and R$_{mean}$
that correlate with the T$_{c}$ slope.
These differences might explain the differences in T$_{c}$ slope distribution, suggesting that there are no
signatures of planet formation in the observed abundance trends with the condensation temperature.

We may conclude that the T$_{c}$ slope depends on the age of a star at a fixed galactocentric radius.
At the same time, stars with smaller galactocentric radii show stepper  T$_{c}$ slopes at a fixed time (age).
In other words, the age and galactic birth place may determine the chemical pattern of the stars.

\begin{acknowledgements}

{We gratefully thank Misha Haywood and Ivan Minchev for fruitful comments and discussion.
We gratefully acknowledge the (second)
anonymous referee for a very exhaustive and honest review of the manuscript that led to many constructive comments and suggestions.
This work was supported by the European Research Council/European Community under the FP7 through Starting Grant agreement
number 239953. V.Zh.A., S.G.S., and E.D.M are supported by grants SFRH/BPD/70574/2010,
SFRH/BPD/47611/2008, and SFRH/BPD/76606/2011 from the FCT (Portugal), respectively.
NCS also acknowledges support in the form of a Investigador FCT contract funded by FCT/MCTES (Portugal) and POPH/FSE (EC).
G.I., S.B.L, and J.I.G.H. acknowledge financial support from the Spanish Ministry project MINECO AYA2011-29060,
and J.I.G.H. also received support from the Spanish Ministry of Economy
and Competitiveness (MINECO) under the 2011 Severo Ochoa Program MINECO SEV-2011-0187.
PF is supported by the FCT and POPH/FSE (EC) through an Investigador FCT contract with application reference IF/01037/2013 and
and POPH/FSE (EC) by FEDER funding through the program "Programa Operacional de Factores de Competitividade - COMPETE.
}
\end{acknowledgements}

Solar and Stellar Astrophysics (astro-ph.SR); Earth and Planetary Astrophysics (astro-ph.EP)


Sound Event Detection in Urban Audio With Single and Multi-Rate PCEN

Abstract: Recent literature has demonstrated that the use of per-channel energy normalization (PCEN), has significant performance improvements over traditional log-scaled mel-frequency spectrograms in acoustic sound event detection (SED) in a multi-class setting with overlapping events. However, the configuration of PCEN's parameters is sensitive to the recording environment, the characteristics of the class of events of interest, and the presence of multiple overlapping events. This leads to improvements on a class-by-class basis, but poor cross-class performance. In this article, we experiment using PCEN spectrograms as an alternative method for SED in urban audio using the UrbanSED dataset, demonstrating per-class improvements based on parameter configuration. Furthermore, we address cross-class performance with PCEN using a novel method, Multi-Rate PCEN (MRPCEN). We demonstrate cross-class SED performance with MRPCEN, demonstrating improvements to cross-class performance compared to traditional single-rate PCEN.

\section{Introduction}
\label{sec:intro}

Noise suppression is a critical step in acoustic signal detection, particularly so in the case of practical sound event detection (SED) in field recordings.
Popular approaches to this task use convolutional operators, mimicking the methods used in computer vision \cite{8434391}.
When using audio as an input to these methods, the images are typically time-frequency representations (spectrograms).
Traditionally, log-scaling is the primary approach for noise reduction, and is the standard approach when spectrograms are the feature of interest.
For single-source audio in clean acoustic environments, this is sufficient for SED.

However, real world environments typically don't have clean, separated events of a single class; this is particularly true in the case of urban audio.
Field recordings often have multiple sound sources of varying acoustic qualities, leading to varying cross-class performance.
Furthermore, the introduction of auditory deformations can lead to rapid performance degradation \cite{DBLP:journals/corr/WangGHLS16}.

Former research has proposed using per-channel energy normalization (PCEN) as a time-frequency representation to mitigate the effects of background noise, demonstrating its use as an input to convolutional methods in SED.
PCEN has proven beneficial in various single-source tasks, including bioacoustic event detection \cite{Lostanlen_2019}, keyword spotting \cite{DBLP:journals/corr/WangGHLS16}, and vocal detection in music \cite{schlueter2018_ismir}.
The per-channel background suppression characteristics make PCEN an attractive choice in SED in urban environments, as background noise in urban environments tends to be Brownian, rather than Gaussian, an assumption made when log-scaling spectrograms \cite{8170052}.

One constraint of PCEN is its dependence on parameter configurations in relation to the specific acoustic properties of the sound event of interest \cite{8514023}.
This makes PCEN poorly suited for multi-class classification, as a given PCEN parameter configuration will likely suit only one particular class of sound events while performing poorly across other classes.
In this article, we demonstrate the effectiveness of PCEN for SED in urban environments.
Furthermore, we propose a new approach to acoustic event detection that combines the foreground separation characteristic of PCEN while preserving multi-class performance: Multi-Rate PCEN (MRPCEN).
By computing a multi-layered PCEN spectrogram at different parameter configurations, we gain the advantages of PCEN, particularly robustness to varying acoustic conditions and auditory deformations without the performance loss associated with cross-class performance.

\section{Audio Representation}
While log-scaling mel-spectrograms is a simple and computationally efficient method of range compression, it has limited effects in certain conditions, particularly those in which background noise is non-Gaussian noise.
PCEN has been used as a pre-processing method for time-frequency representations that reduces the effects of noise on convolutional neural networks by Gaussianizing the distribution of magnitudes across mel-frequency spectrogram coefficients \cite{DBLP:journals/corr/WangGHLS16}.
\begin{figure}
    \centering
    \includegraphics[width=\columnwidth]{spec_demo.png}
    \caption{A comparison of standard SED audio featurizations compared to Multi-Rate PCEN. Note the unevenly distributed noise in the log-mel spectrogram, with more noise in the lower end of the spectrum (due to Brown noise found in urban environments), which is removed in the PCEN spectrogram.}
    \label{fig:pcen_spec_comparison}
\end{figure}

\subsection{PCEN}
PCEN is a sequence of audio processing steps on a spectrogram $\textbf{E}(t,f)$ using adaptive gain control scaled by the response an autogregressive filter $\phi_T$, followed by dynamic range compression. It takes the form:

\begin{equation}
    \textbf{PCEN}(t,f) = \left(\frac{\textbf{E}(t,f)}{(\epsilon + (\mathbf{E} * \phi_T)(t,f))^\alpha} + \delta \right)^r - \delta^r
\end{equation}

\noindent The gain $\alpha$ scales the smoothing of the spectrogram.
$\delta$ and $r$ are the bias and power of the dynamic range compression that PCEN provides.
We can see that the filter response $\mathbf{M}(t,f) = (\mathbf{E}*\phi_T)(t,f)$ scales our PCEN magnitude by:

\begin{equation}
    \mathbf{M}(t,f) = (\mathbf{E}*\phi_T)(t,f) = s \mathbf{E}(t,f) + (1-s)\mathbf{M}(t-\tau,f)
\end{equation}

\noindent where $0<s<1$ is the weight of 1\textsuperscript{st} order autogregressive filter (AR(1)), and $\tau$ is the discretized step defined by the hop size of the input spectrogram ($\epsilon$ is an offset term).
The resulting filter is a low-pass filter of 0dB gain and a cuttoff frequency of $\omega_c = \frac{2 \pi \tau}{T} = \arccos{(1-\frac{s^2}{2(1-s)})}$ at 3 dB, and sidelobe falloff of
10 dB per decade near $\omega_c$ \cite{8514023}.

We consider this filter response $\mathbf{M}(t,f)$ an approximation to the relative magnitude of stationary background noise at each frequency band $f$.
The effect of scaling our output PCEN spectrogram by the reciprocal of this response is amplifying/magnifying the response of foreground events and suppressing background events on a per-channel basis.
In regimes where the background noise isn't Gaussian, such as the case with Brownian noise in urban environments, this decorrelates noise across different frequency bands.

A critical tuning parameter of this smoothing of the autoregressive filter is the rate parameter $T$, which defines the cutoff frequency by $\omega_c = \frac{2 \pi \tau}{T}$.
Setting $T$ too low will reduce the noise reduction effects of PCEN, but setting it too high will suppress the sound event of interest, especially if the sound event of interest is stationary.
Prior practical recommendations for the rate parameter are defined by the stationarity of the sound, the frequency range, and the chirp rate.
For single sources in a consistent acoustic environment, tuning the rate parameter $T$ to these specifications proves sufficient for audio detection and classification tasks \cite{8514023}.

\subsection{Multi-rate PCEN}
In a field setting, with several audio classes with varying acoustic properties and variable recording conditions, a singular ideal value for $T$ is much harder to identify.
For example, in the case of urban audio, the combination of short, fast decay sounds (such as gunshots or dog barks) compared to longer ambient sounds (e.g. sirens, air conditioners) lead to widely varying preferred values of $T$.
To capture information across varying regimes of our rate parameter, $T$, we take inspiration from 3-channel RGB images in computer vision.
These three separate but correlated feature maps can be passed to a convolutional neural network (CNN) which can use this multi-frequency image to make predictions that a black-and-white image could not.

We replicate this multi-regime approach, but instead of frequency responses, we vary our rate parameter $T$ in each layer, as shown in figure \ref{fig:pcen_spec_comparison}.
Each $i$\textsubscript{th} layer of the image has a differing level of gain control applied based on the cutoff threshold from $T_i$.
By using multiple logarithmically-scaled values of $T_i$, we can produce a multi-layered image that captures information of sounds at varying decay lengths.
This ensures that sound events that may be suppressed using one value of $T_i$ will be preserved in another.
The resulting multi-channel image can be used as an input to a CNN much like a multi-channel color image is used in machine vision.
By incorporating information at varying degrees of gain control, our model preserves the robustness of PCEN without degrading multi-class performance.

\section{Experimental Design}

\subsection{UrbanSED}
We used the UrbanSED dataset to evaluate the performance on this method for the task of sound event detection\cite{8170052}.
UrbanSED contains 10,000 synthetic soundscapes of 10 distinct sound classes, with each class having approximately 1000 instances per class, drawn from the UrbanSound8K dataset.
Each soundscape contains 1-9 time/class labeled foreground events with additive background Brownian noise.
The dataset is pre-split in a 6-2-2 training, validation, and evaluation subsets.
Because UrbanSED is synthetic, we can ensure no spurious unlabeled audio is included in the soundscapes, making it a standard benchmark for state-of-the-art SED models.

\subsection{Data augmentation}
Because one of the primary predicted benefits of PCEN is improved robustness to audio deformations, we augmented our dataset with several reverberant duplicates of UrbanSED.

In this implementation, reverb is modeled as a convolution of a source signal and an impulse response of a given acoustic environment.
Augmenting a given audio clip with reverb is done by computing the convolution of the impulse and the audio file.
This effectively computes the response of the audio clip in the acoustic environment associated with the impulse.
We used 6 distinct reverb responses, three impulse responses recorded in different acoustic environments, and three synthetic impulse responses of white noise with an exponential decay envelope $e^{-t/\tau_c}$.

The three real impulse responses were recorded in a bedroom\footnote{``My Bedroom'', https://freesound.org/people/Uzbazur/sounds/382907/}, an alleyway\footnote{``alley.wav'', https://freesound.org/people/NoiseCollector/sounds/126804/}, and a tunnel\footnote{``tunnel\_2013'', https://freesound.org/people/recordinghopkins/sounds/175358/}, each with increasing decay times.
The three synthetic impulse responses had decay time constants $\tau_c$ of $0.1$s, $0.3$s, and $0.5$s.
In addition to reverb augmentation, we also duplicated our dataset by pitch shifting each sample by \{$\pm 1, \pm 2$\} semitones.
We applied pitch shifts and convolutional to reverb to our dataset using using \textit{MuDA}, a library for musical data augmentation \cite{muda2015}.

\subsection{Audio featurization}
Training audio was processed using librosa 0.7.2 \cite{brian_mcfee-proc-scipy-2015} which generated a PCEN spectrogram for each audio sample and rate parameter $T$.
Each spectrogram had the following parameters: sampling rate $44.1$kHz, window size $1024$ samples, hop length $512$ samples, and $128$ mel-frequency bands.
Our resulting frame rate was $44.1$kHz$/512 \approx 86$Hz.
At a 10 second length sample length, each spectrogram had a horizontal length of 862 samples per band.

Our rate parameters were 10 logarithmically spaced values, ranging from $2^0=1$ to $2^9=512$, corresponding to averaging over windows ranging from $10$ms to $6$s.
The adaptive gain control bias was set to $\epsilon=10^{-6}$, the gain was set to $\alpha = 0.98$, the dynamic range compression bias was set to $\delta=2$, and the compression power $r=0.5$ for all rate configurations; these parameters were chosen based on prior work in audio classification \cite{Lostanlen_2019}.
Each audio sample had 10 PCEN spectrograms, one per rate parameter.
Individual layers/subsets were selected in training and evaluation.
Following the channels-last conventions of Keras, the resulting PCEN spectrograms have the shape of (128,862,10), and the resulting log-mel spectrograms had a shape (128,862,1).

\subsection{Network architecture}
The network architecture is inspired by the $L^3$ audio sub-network \cite{DBLP:journals/corr/ArandjelovicZ17} for discriminative audio-video correspondence embeddings.
This architecture has demonstrated success in its use for classification and predictions at a fine time-resolution with mel-spectrogram inputs\cite{8434391}.
We follow the implementation of this architecture found in \cite{8434391}, with the input layer adjusted to accept multi-layer images to ensure that all methods were being evaluated on near-identical architecture.
Additional parameters in the input layer may lead to lower performance due to additional depth requirements, but this is assumed to be negligible for this application.

Models were built and trained using Keras 2.3.1 \cite{chollet2015keras}, with Tensorflow 2.2.0 \cite{tensorflow2015-whitepaper}.
The model was trained using the Adam optimizer \cite{kingma2014method}, using UrbanSED's pre-folded training and validation sets.
The loss function was binary cross-entropy on a per-class per-frame bases.
The validation metric was accuracy.
If no improvement in the accuracy was seen in 10 epochs, the learning rate would be reduced.
Early stopping was implemented if no improvement was seen after 30 epochs.
Evaluation was handled via the sed\_eval package \cite{sed_eval_art}, which computed segment-based classification metrics, including precision, recall, and error rate.
We used the F1-score, the harmonic mean of precision and recall, as our main metric of effectiveness.
These metrics were computer both per-class and overall.
For reproducability, the implementation and experimental framework, including data preprocessing, model training, and evaluation, is publicly available on github\footnote{https://github.com/ChrisIck/pcen-t-varying}.

\subsection{Data preparation}
We trained and evaluated models across multiple datasets to test robustness and stability to audio deformation.
We primarily used the \textit{dry} dataset as our baseline dataset, which was UrbanSED with no reverb augmentation, to ensure we were achieving state of the art performance seen in \cite{8434391}.
The \textit{realreverb} set is built from the reverb-augmented audio using the 3 real recorded impulse responses, and is the primary set we evaluate, as this is closest to real-world conditions.
The \textit{simreverb} set, using the remaining 3 synthetic impulse responses, demonstrated the strongest separation between PCEN and log-mel models, but the results don't generalize as well as when using real impulse responses.
All models were trained on data augmented by $\{\pm 1, \pm 2\}$ pitch-shifts, but were validated and evaluated on non-shifted data.

\subsection{Rate parameters}
We trained and evaluated models on fixed sets of rate parameters.
We experimented with single rate parameter models and varying sets of multi-layered models, each with rate constants in $T_k = 2^k$ for $k \in \{0,1,...,9\}$.
We tested a total of 44 unique rate parameter configurations for our PCEN models.
This included single-rate parameter models and $n$-layer models including rate parameters from $T_i$ to $T_{i+n-1}$.
As a baseline, we also computed a model using a traditional log-scaled mel-spectrogram, as seen in previous literature for this application and dataset \cite{8434391}.

\subsection{Training and evaluation}
Each set of models was trained and evaluated on both datasets to see how stable each model was to reverberation conditions that were both contained in the respective training dataset, as well as in conditions distinct from those in the training set.
Evaluation was done on bootstrapped subsamples of each evalution set, sampling 100 evaluation examples with replacement, 100 times per model and dataset.
Evaluations were completed both in overall micro-averaged performance metrics across classes, as well as on a per-class basis.

\section{Results}
\begin{figure}
    \centering
    \includegraphics[width=\columnwidth]{singlerate_compare_rr.pdf}
    \caption{F1-score of single-rate PCEN models at various $T$ values (blue dots), compared to 8-layer MRPCEN (dashed green) and log-mel (dotted red) models evaluated on the 100 bootstrapped samples of the \textit{realreverb} dataset.}
    \label{fig:singlerate_compare}
\end{figure}

\subsection{Per-class performance}
We demonstrate PCEN's robustness to noise by comparing performance of models trained and evaluated on the \textit{realreverb} dataset.
We see in figure \ref{fig:singlerate_compare} that single rate PCEN demonstrates performance at or above log-scaled spectrograms, depending on its choice of rate parameter $T$ and the class of the target event.
Certain classes, like the jackhammer and gunshot classes, prefer higher $T$ values in the neighborhood of $T = 2^7$ where mid-low frequency is filtered out with the adaptive gain aspect of PCEN.
Similarly, for higher and full-band frequency stationary events, such as air conditioners and sirens, lower values of $T$ in the neighborhood of $T \in [2^1,2^3]$ perform better, as the low-frequency noise doesn't interfere with the relevant frequency bands.
However, it is crucial to note that due to the diversity of sonic characteristics of these different classes, there does not exist a single ideal value of $T$ that produces the strongest performance overall.

\subsection{Cross-class performance}
MRPCEN successfully outperforms most single-rate models in cross-class performance, as seen in figure \ref{fig:singlerate_compare}.
The plot shows per-class performance of each model evaluated on the \textit{realreverb} dataset, compared to a log-mel based model and an 8-layer MRPCEN.
In comparison, MRPCEN model preserves information at multiple rates, it consistently performs at or above the majority of the single rate models.
Furthermore, MRPCEN performs at or above the level of the standard log-mel spectrogram with the exception of the ``car\_horn'' and ``jackhammer'' classes, providing it with the highest overall performance metric on this dataset (compared to the strongest performing single-rate PCEN model $T=2^1$) in table \ref{fig:overall_f1_table}.

\begin{table}
    \centering
        \begin{tabular}{l|p{1.7cm}|p{1.7cm}|p{1.7cm}}
        Model & Overall F1 (\textit{realreverb}) & Overall F1 (\textit{simreverb}) & Overall F1 (All Data)\\
        \hline
        Logmel & 0.334 & 0.175 & 0.274\\
        PCEN & 0.345 & 0.167 & 0.268\\
        MRPCEN & \textbf{0.356} & \textbf{0.204} & \textbf{0.285}
    \end{tabular}
    \caption{Overall Micro-averaged F1 Scores on the \textit{realreverb} and \textit{simreverb} datasets}
    \label{fig:overall_f1_table}
\end{table}

\section{Conclusion}
The results here show PCEN as a viable alternative to log-mel spectrograms, showing equivalent or improved performance depending on rate parameter choice, which in turn depends on the acoustic characteristics of the target sound event and acoustic environment.
We can also see that MRPCEN provides cross-class performance improvements over single-rate PCEN models.
MRPCEN permits simultaneous prediction of differing classes with distinct acoustic characteristics in a single model.
In a field setting, this will lead to less per-class knowledge and tuning expertise needed to effectively deploy a model that can perform well in multi-class applications.

\vfill\pagebreak
\section{Acknowledgements}
This work was supported in part through by NSF award 1955357, and in part by the NYU IT High Performance Computing resources, services, and staff expertise.

\bibliographystyle{IEEEbib}

Audio and Speech Processing (eess.AS); Machine Learning (cs.LG); Sound (cs.SD)