The SDSS-V Black Hole Mapper Reverberation Mapping Project: C iv Broad Absorption Line Acceleration in the Quasar SBS 1408+544

We present the results of an investigation of a highly variable C iv broad absorption line (BAL) feature in spectra of the quasar SBS 1408+544 (z = 2.337) that shows a significant shift in velocity over time. This source was observed as a part of the Sloan Digital Sky Survey (SDSS) Reverberation Mapping project and the SDSS-V Black Hole Mapper Reverberation Mapping project, and has been included in two previous studies, both of which identified significant variability in a high-velocity C iv BAL on timescales of just a few days in the quasar rest frame. Using ∼130 spectra acquired over 8 yr of spectroscopic monitoring with SDSS, we have determined that this BAL is not only varying in strength, but is also systematically shifting to higher velocities. Using cross-correlation methods, we measure the velocity shifts (and corresponding acceleration) of the BAL over a wide range of timescales, measuring an overall velocity shift of Δv=−683−84+89 km s−1 over the 8 yr monitoring period. This corresponds to an average rest-frame acceleration of a = 1.04 −0.13+0.14 cm s−2, though the magnitude of the acceleration on shorter timescales is not constant throughout. We place our measurements in the context of BAL-acceleration models and examine various possible causes of the observed velocity shift.

1. INTRODUCTION Broad absorption lines (BALs) in the spectra of quasars are thought to originate in winds that are launched from quasar accretion disks (e.g., Murray et al. 1995;Proga 2000;Higginbottom et al. 2014;Naddaf et al. 2023).BALs are defined as absorption features with velocity widths wider than 2000 km s −1 (Weymann et al. 1991;Hall et al. 2013) and are found in roughly 10-15% of all optically selected quasars (Gibson et al. 2009;Allen et al. 2011).The winds/outflows that produce BALs may play an important role in the evolution of galaxies -if these outflows are sufficiently energetic, the gas can produce significant feedback that interferes with star formation within the host galaxy and/or further growth of the supermassive black hole (e.g., Di Matteo et al. 2005;Moll et al. 2007;King 2010).Characterizing these BALs and their environments to constrain models for how they are produced, how they evolve, and how they affect their galaxies, is thus important for our understanding of galaxy evolution.
BALs are seen across a wide variety of ionization species (Turnshek 1984;Arav et al. 2001), one of the most common of which is C iv, which appears as a doublet at rest-frame wavelengths of 1548.20 Å and 1550.77Å. C iv is one of the higher-ionization BAL species that appears the most often in quasar spectra, and is high enough in abundance that it is often saturated.It is thus often probed using Si iv as a tracer, as Si iv is much lower abundance and thus usually unsaturated and can be used for density estimations (e.g., Arav et al. 2018).Xu et al. (2019) use photoionization modeling to determine the physical conditions of BAL outflows using Si iv and find that high-ionization BALs such as C iv can probe regions of gas with electron densities ranging from 10 3 -10 5.5 cm −3 , column densities of log N H = 20 -22.5, and log U H between −2 and 0 (see their Figure 4); in addition, their photoionization models return a range of temperatures for these outflows ranging from 15,000 K to 20,000 K.Many of the studies that have used photoionization modeling to determine the distance of these outflows find that the outflows are at distances of > 500 pc from the central source (see Arav et al. 2018 and references therein).
BALs are variable on rest-frame timescales ranging from days to years (e.g., Lundgren et al. 2007;Gibson et al. 2008;Capellupo et al. 2012;Filiz Ak et al. 2013;Vivek et al. 2014;Grier et al. 2015Grier et al. , 2016)).The most commonly observed modes of variability are changes in strength and/or profile shape.This observed variability allows us to place constraints on the geometry, distance, and dynamics of the outflows themselves, possibly revealing information on the energetics of the outflows to inform models of quasar/host-galaxy feedback (e.g., Arav et al. 2013).In addition to changes in equivalent width and profile shape, monolithic velocity shifts, or "acceleration", of BAL outflows have also been reported.The presence (or lack) of observed acceleration in BALs provides an important test for models describing the production and evolution of BAL outflows, as some of these models predict visible acceleration (e.g., Higginbottom et al. 2014).
Over the past two decades, there have been a handful of reports of observed velocity shifts (implying BAL acceleration or deceleration) in several different studies (Vilkoviskij & Irwin 2001;Rupke et al. 2002;Gabel et al. 2003;Hall et al. 2007;Joshi et al. 2014;Grier et al. 2016;Joshi et al. 2019;Lu & Lin 2019, 2020;Xu et al. 2020;Byun et al. 2022).However, velocity shifts caused by an increase or decrease in speed of outflowing gas are difficult to identify.First, because BALs are quite variable in line profile, it can be difficult to differentiate between shifts caused by velocity-dependent line-profile variability (due to changes in the ionization state of the absorbing gas or changes in the column density coverage, for example) and shifts caused by an actual change in the speed of the outflow.In addition, previously measured velocity shifts have mostly been small in magnitude over short timescales, with typical velocity shifts of only a few hundred km s −1 over a few years in the quasar rest frame.We thus expect to require long time baselines to observe a significant velocity shift.
The Sloan Digital Sky Survey Reverberation Mapping Project (SDSS-RM; Shen et al. 2015) and the SDSS-V Black Hole Mapper Reverberation Mapping Program (BHM-RM; Kollmeier et al. 2019) provide us with an excellent opportunity to explore BAL variability, and acceleration, in quasars.The SDSS-RM program observed a single field of ∼850 quasars from 2014-2020; these observations began as a part of the SDSS-III and SDSS-IV surveys (Eisenstein et al. 2011;Blanton et al. 2017).Roughly 380 of these quasars continue to be monitored as a part of the SDSS-V BHM-RM program, further extending the time baseline of these observations.While the primary goal of the SDSS-RM and BHM-RM monitoring programs is to measure black-hole masses using the technique of reverberation mapping, there are roughly 90 quasars in this sample that show BAL features.With more than 100 spectral observations of these quasars over 8 years (and counting), this survey is very well suited for studies of BAL variability, allowing us to explore it on both short (few-day; e.g., Hemler et al. 2019) and long (several-year) timescales.
One of the sources observed by the SDSS-RM and BHM-RM programs is the quasar SBS 1408+544 (SDSS J141007.72+541203.6,hereafter referred to as RM 613).This quasar has a redshift of z = 2.337 ± 0.003 (Shen et al. 2019), an apparent i-band magnitude m i = 18.1 (Alam et al. 2015), and an absolute magnitude M i = − 27.69.This source was somewhat serendipitously discovered to show strong variability in the equivalent width of its C iv BAL on very short timescales - Grier et al. (2015) studied the first 32 observations of this quasar during the first year of monitoring (each observation was taken on a different night, hereafter referred to as an epoch) and found that the equivalent width of the BAL was changing significantly on timescales down to just 1.20 days in the quasar rest-frame.Hemler et al. (2019) also included RM 613 in their sample study of C iv BAL variability in ∼30 SDSS-RM BAL quasars, which included four years of SDSS-RM monitoring.They confirmed the short-term variability that was originally reported by Grier et al. (2015), but in the subsequent three years of observations, this BAL weakened somewhat and did not show additional dramatic variability.We have since acquired roughly 70 additional spectra of the SDSS-RM field over four additional years of monitoring by SDSS, so we are following up on all of the sources that were found to have significant short-timescale variability by Hemler et al. (2019) to investigate whether they continued to show strong variability.We here revisit the quasar RM 613 as a part of this followup effort.
In Section 2, we discuss the observations, the preparation of the spectra for analysis, and our continuumnormalization procedure.Section 3 describes the measurements made and the tests we performed, and Section 4 includes a discussion of these results, their implications, and their relevance to models of BALs.

DATA AND DATA PREPARATION 2.1. Spectral Data
The spectra used in this study are from the SDSS-RM project (e.g., Shen et al. 2015), which was carried out as a part of the SDSS-III and SDSS-IV programs (Eisenstein et al. 2011;Blanton et al. 2017) from 2014-2020.In addition, we obtained 41 spectra of this source as a part of the SDSS-V BHM-RM program (Kollmeier et al. 2019), which began observations in 2021.There is also an additional spectrum of this source taken as a part of the SDSS-III survey in 2013 May (before the SDSS-RM program began) that was released to the public as a part of the SDSS Data Release 12 (DR12; Alam et al. 2015); we include this spectrum in our study as well.In total, we have 132 spectra of this quasar taken with the 2.5-m SDSS telescope at Apache Point Observatory with the BOSS spectrograph (Gunn et al. 2006;Dawson et al. 2013;Smee et al. 2013) spanning roughly nine years (eight years of dedicated monitoring, plus the early SDSS spectrum taken a year earlier).The BOSS spectrograph covers a wavelength range of roughly 3650-10,400 Å with a spectral resolution of R ∼ 2000, which results in a velocity resolution of ∼69 km s −1 in the C iv region of the spectrum.The spectra from the first two years of monitoring (2013 and 2014) were processed with the standard SDSS-III pipeline (version 5 7 1) and the remaining years of data were processed with the updated SDSS-IV eBOSS pipeline (version 5 10 1). Figure 1 shows the mean spectrum of RM 613, created from these spectra.
With the recent upgrade to the SDSS-V robotic fiber positioning system (FPS), our original field of 849 quasars had to be reduced to approximately 380 sources for continuing observations.RM 613 was dropped from the SDSS-V monitoring because its position in the SDSS-RM field is outside the FPS field of view, so we have not been able to continue monitoring this particular source with SDSS-V.However, we were able to obtain a spectrum of RM 613 on 30 May 2023 with the Hobby-Eberly Telescope (Ramsey et al. 1998;Hill et al. 2021) using the Low Resolution Spectrograph 2 (LRS2; Chonis et al. 2016).Figure 2 shows that the C iv BAL in question has weakened significantly at this point.This spectrum was processed separately from the SDSS spectra; see Appendix A for details.
We also note that an additional spectrum of this source was taken in 1991 by Chavushyan et al. (1995).Grier et al. (2015) inspected this spectrum and noted that the high-velocity C iv BAL was not detectable in the spectrum at the time, though the SNR was too low for its absence to be established definitively.We do not include this spectrum in our analysis, as the signal is too low for useful measurements to be obtained.
Visual inspection of the mean spectrum (Figure 1) confirms that this quasar hosts a prominent highvelocity C iv BAL feature that is detached from the C iv emission line.Following Grier et al. (2015), we will henceforth refer to this BAL as "Trough A".In addition to this higher-velocity C iv BAL, we see two narrower C iv absorption systems that are superimposed onto the C iv emission line (hereafter "Trough B" for the middle feature and "Trough C" for the narrower feature that lies nearly at line center; see Figure 1).We see no BALs in lower-ionization transitions (e.g., Mg ii or Al iii), which means this source is considered a "high-ionization" BAL quasar.There is a hint of Si iv absorption at similar velocities as C iv but the absorption is too weak to formally meet the definition of a BAL.While it is visible in the mean spectrum, this absorption is too shallow for us to well measure the properties of this feature in individual spectra.In addition to the difficulties presented by the shallowness of the line, its blue edge is contaminated by Oi/ Si ii emission, and there is a strong, narrow absorption line located on the red edge.Both of these features interfere with attempts to make measurements of the Si iv absorption even in mean/stacked spectra.We also see some N v absorption just redward of Lyα, but it is too contaminated by Lyα absorption for us to include in our study.

Data Preparation
Prior to any further processing, we visually inspected all spectra and noticed a few epochs with significant excess noise.We quantify the SNR of each spectrum at rest-frame 1700 Å by measuring the median SNR of pixels between rest-frame 1650-1750 Å (SNR 1700 ), which corresponds to 5506-5840 Å in the observed frame for RM 613.Based on visual inspection, we determined that spectra with SNR 1700 < 3 were unlikely to provide any useful constraints, so we excluded all epochs with SNR 1700 below this threshold.This excluded three epochs: 7, 53, and 115.This brings the total number of usable BOSS spectra down to 129 (including the early BOSS spectrum).
We first cropped off all pixels at wavelengths less than 3650 Å and greater than 10300 Å (in the observed frame), as both the blue and red edges of the spectra showed significant noise and telluric contamination.We then searched for pixels that may have problems, as  flagged by the SDSS pipeline using the bitmasks provided with the spectra.We linearly interpolated over any pixels that were flagged by the SDSS "ANDMASK" as having issues.The uncertainties on the interpolated pixels were multiplied by a factor of 10 to represent the increased uncertainties due to interpolation.We follow previous work (e.g.Gibson et al. 2009;Grier et al. 2015) and correct for Galactic extinction and reddening in the spectra, adopting a R V = 3.1 Milky Way extinction model (Cardelli et al. 1989) and R V values following Schlafly & Finkbeiner (2011).We then shifted the spectra to the quasar rest frame using a redshift of z = 2.337 (Grier et al. 2015).All further discussion/analysis of these spectra will refer to the spectra in the quasar rest frame.

Continuum and Emission-Line Fits
To isolate the variability of the BAL from the variability of the rest of the quasar, we fit a continuum model to each spectrum.We follow previous work (e.g., Filiz Ak et al. 2012) and model the quasar continuum as a reddened power law with an SMC-like reddening coefficient (Pei 1992).We fit the continuum using a nonlinear least-squares algorithm, selecting four "line-free regions" -wavelength ranges that are largely uncontaminated by strong emission and absorption features (1270 Å -1290 Å, 1700 Å -1750 Å, 1950 Å -2050 Å, and 2200 Å -2300 Å).An example spectrum with its continuum fit is shown in Figure 3.We calculated the uncertainties in the continuum fits using "flux randomization" Monte Carlo iterations, where we altered the flux of each individual pixel by a random Gaussian deviate based on the size of its uncertainty.We then fit the continuum to the new altered spectrum, and repeated this process 100 times to determine the standard deviation of the model continuum flux at each pixel, which we adopt as the uncertainty in the continuum fit.
The primary BAL of interest in our work, Trough A, is at a sufficiently high velocity that it appears detached from the C iv emission line.However, there are two additional absorption features that lie on top of the C iv emission line; in order to examine their behavior, we need to isolate the C iv emission line from the BALs.To do this, we follow Hemler et al. (2019) and fit an emission-line profile to each individual spectrum and divide it out with the continuum.We chose to use a Voigt profile (e.g., Gibson et al. 2009), which was found by Hemler et al. (2019) to be a good fit in this particular source.We used an iterative fitting technique to exclude wavelength bins where the flux deviated from the fit by more than 3σ and manually excluded the regions that showed significant absorption.As expected, the C iv emission line model was not so wide that it affects the high-velocity Trough A BAL; as such, our primary analysis, which focuses on Trough A, does not rely at all on our emission-line fits.
We also investigated the possibility that the Si iv λ1393 emission line, with line center at 1393.755 Å, may interfere with our measurements.We fit a Voigt profile to the Si iv line using the same procedure as with C iv , but found that the red wing of the Si iv line ended significantly blueward from the blue wing of Trough A; there is no overlap between the Si iv emission line and the BAL feature at any epoch.As with the C iv line, the Si iv line thus does not affect any of our measurements of Trough A. The Si iv fits are shown in Figure 3 for demonstrative purposes, but are not included in the majority of our analysis.
We added the continuum and C iv emission-line models together and divided the original spectra by this combined continuum+emission-line fit to obtain a set of "normalized" spectra.The uncertainties from the continuum fit were propagated along with the spectral uncertainties to determine the final uncertainties on the normalized spectra.All subsequent measurements and analysis were performed on the normalized spectra unless otherwise noted.Figure 4 shows the mean normalized spectrum, focused on the C iv region of the spectrum.Because we only fit the continuum and C iv emission line, other emission features are still visible at the edges of this range (we see Si iv around 1400 Å and lowlevel contributions from a wide variety of species redward of C iv); however, none of these features has any effect on our measurements of Trough A, as the local continuum around Trough A is well fit by the continuum.In addition to Troughs A, B, and C, there are a number of narrow absorption features present in the spectrum as well, some of which overlap on top of Trough A and B. These have been identified as intervening C iv and Si iv systems at a variety of redshifts.

BAL Measurements
Our first goal was to examine the variability of the C iv BAL Trough A with the additional four years of monitoring with SDSS that were not included by Hemler et al. (2019).To do this, we first characterize the properties of Trough A by measuring its rest-frame equivalent width (EW, in units of Å), as defined by where f n represents the continuum-normalized flux.We also measure the absorbed-flux-weighted velocity centroid in units of km s −1 , and the mean fractional depth (d, in units of normalized flux).We first smoothed the spectra using a boxcar smoothing algorithm over five pixels; this smoothed spectrum was used to determine the upper and lower wavelength limits of the BALs, defined as the location where the flux in the BAL reaches 90% of the continuum flux (corresponding to a normalized flux density of 0.9 as per standard BAL definitions/conventions; see, e.g., Gibson et al. 2009) on either side of the BAL. Figure 5 shows these boundaries over the course of the monitoring period.We then measured the properties of the BAL  within these outer limits, using Monte Carlo randomizations to measure the uncertainties in these parameters.These measurements are shown in Figure 6. Figure 6 indicates how the BAL varied over the monitoring period using these three measurements.We again confirm the strong short-timescale variability observed by Grier et al. (2015) during the first year of observations.This was followed by a period of about 4 years where the BAL remained in a weaker state before it again strengthened.The mean depth shows the same variability as the EW.However, inspection of v cent over time shows a very strong trend toward higher (more negative) velocities over the monitoring period.This could be indicative of acceleration of the BAL.Examination of the upper and lower limits of the BAL as a function of time (Figure 5) confirms that the upper and lower bounds of the BAL are also slowly shifting blueward over time.
Figure 7 shows a "trailed spectrogram" that includes all 129 of the SDSS spectra included in our study as a way of visualizing the normalized flux over time.For visualization purposes, we used the continuum-normalized spectra without the C iv emission-line fit to produce this figure.Examination of this plot confirms our conclusions based on the measurements of the BAL.First, the blueward shift of the BAL over the course of the observations (as suggested by our v cent vs. time plot) is readily apparent.The trailed spectrogram also indicates the relative weakening of the absorption during the middle 3-4 years of monitoring, followed by a regaining of strength in the final three years of monitoring.A similar trend in strength is observed in the Trough B absorption feature (appearing at roughly 1535 Å), though Trough C, the narrow doublet feature that lies close to line center, does not appear to vary significantly in strength.Faint, thin, vertical blue lines that are visible in Figure 7 (particularly those just blueward of Trough A) represent the narrow intervening absorption systems that we identified (Section 2.3).
Despite the significant changes in EW, Grier et al. (2015) found no significant line-profile variability on short timescales during the first year of monitoring.However, our Figure 7 indicates a change in the shape of the BAL over the longer time period.To examine how the shape of the BAL changed over the longer time baseline spanned by our data, we combined each year of monitoring, producing a mean spectrum for each individual year, and inspected the shape of the line in these higher-SNR spectra.Figure 8 shows that the average line profile remained similar through the first several years of monitoring; however, over the last 2-3 years of monitoring, the trough changed from having two distinct sub-troughs to a more uniform, single-trough feature.The slow blueward shift in velocity of this BAL is still quite apparent in Figure 8.Interestingly, our recent HET spectrum (see Appendix A) shows that the trough had returned to its two-pronged structure by mid-2023, with Trough A presenting as two distinct troughs where the normalized flux returns to a value of 1 between them (formally, this feature would be identified as two mini-BALs in this spectrum rather than a BAL).
Previous studies of BAL acceleration typically only include 2-3 spectra, though some studies have had as many as five (Byun et al. 2022).The more spectra, the more clear it is whether or not the observed shift in velocity is "real" or not -i.e., whether the shift may be due to velocity-dependent variability across the trough, or is actually a shift in velocity that could be caused by acceleration of the outflow itself.With 129 epochs spanning 9 years, it is clear that this BAL is undergoing an actual shift in velocity (underneath any additional variability in shape and strength) rather than undergoing velocity-dependent variations that mimic a velocity shift.

BAL Velocity-Shift Measurements
To quantify the velocity shift in our C iv BAL, we adopt two different measurements.First, we measure the change in v cent between spectra, ∆v cent .This is a straightforward way to measure the change in velocity of the trough over time and quantify the acceleration of the BAL. Figure 9 shows the change in v cent between all 8256 possible pairs of spectral epochs (each of the 129 spectra paired with every other spectrum following it).We see a clear trend of an increasing velocity shift with an increase in ∆t, supporting the idea that the BAL is in fact accelerating.We do note a slight "bifurcation" in the ∆v cent measurements shown in Figure 9 that arises due to different pairs of epochs that have similar ∆t (for example, measurements between day 1 and day 400, and between day 400 and day 800, will both have the same x-value on the plot).This suggests that the magnitude of the velocity shift changes throughout the monitoring period.
While v cent provides a way to characterize a velocity shift, in principle it is also sensitive to changes in the line profile.This means that velocity-dependent variability within the trough can cause a change in v cent even if the trough itself remains within the same boundaries.While we do not see significant line-profile variability on short timescales in the BAL of RM 613, the shape of the BAL does appear to change during the latter half of our monitoring period.By relying solely on v cent to quantify the acceleration, we risk the possibility that our measurements are somewhat skewed by this shape change in later years.
To help minimize the possibility that our measurements of acceleration are affected by line-profile changes, we adopt a second method of quantifying the velocity shift between two epochs: the Interpolated Cross Correlation Function (ICCF), as adapted for searches for velocity shifts in BALs by Grier et al. (2016).The ICCF procedure (e.g., Peterson et al. 1998) was originally developed to measure time delays between light curves in reverberation mapping data -however, the procedure can be adapted to search for velocity shifts in spectra as well.Using the CCF to measure the velocity shift may be less affected by changes in the line strength and shape than v cent because the ICCF will measure a lower correlation coefficient when this variability is present.
We first isolate the C iv BAL in our normalized spectra by cropping the spectra down to include only the BAL itself and about 2000 km s −1 of "padding" on either side of the BAL (to allow for shifts).We implement the ICCF analysis via the PyCCF software (Sun et al. 2018), which works as follows: The code first measures the correlation coefficient r between two spectra, and then applies a velocity shift (in increments defined by the user) and interpolates the data so that the shifted spectra lie on the same wavelength grid.The coefficient r is then remeasured after applying this shift, and the spectra are shifted again -r is measured after applying all possible velocity shifts within a given range to build up the cross correlation function (CCF).The peak value of the CCF represents the velocity shift at which the two spectra are the most highly correlated.As is commonly done in reverberation-mapping studies, we measure the centroid of the CCF about the peak by including all points with values greater than 0.8r peak to characterize the velocity shift between two spectra (∆v CCF ).
To measure the uncertainties in ∆v CCF , we follow Peterson et al. (1998,2004) and employ Monte Carlo simulations as we did for the continuum fits.We perform 1,000 iterations: The spectra are altered by random Gaussian deviates scaled by their uncertainties, and the CCF is recalculated for each iteration.We then adopt the median of the cross correlation centroid distribution (CCCD) as our best velocity-shift measurement (we also keep track of the peak distribution, or CCPD), and the 1σ uncertainties are calculated corresponding to the 68.3% percentile of the CCCD.We generally measure uncertainties on the order of 1-2 pixels (each SDSS wavelength bin corresponds to ∼69 kms −1 ), depending on the SNR of the individual spectra involved.Figure 10 shows an example pair of spectra and the resulting CCF and CCCD.
As with ∆v cent , we measured the CCF between all 8256 pairs of spectra to explore all possible timescales.Figure 9 shows these measurements as a function of time.We again see significant velocity shifts that increase with ∆t; our detection of acceleration is thus strengthened by the observations of effectively the same acceleration trends using these two independent methods.However, the CCF measured smaller shifts overall than ∆v cent -we suspect that the differences are due to the changing shape of the BAL at the end of the campaign.Because ∆v CCF is less sensitive to these changes in the shape than ∆v cent , it measures a smaller shift between the spectra in pairs where there has a been a substantial shape change between observations.Additionally, the "bifurcation" seen in the ∆v cent measurements in Figure 9 disappears when the CCF is used, further suggesting that the increase in the velocity shifts measured by ∆v cent between epochs in the latter half of the observations is significantly affected by the change in shape of the BAL.

Measurements of the Mean Spectra
With ∼130 spectral epochs and more than 8000 possible pairs of spectra, choosing individual epochs to best quantify the acceleration across the 8 years of monitoring is somewhat challenging.To simplify things, we used the mean spectra from observations in each individual year (Figure 8), which largely have similar noise properties as one another, and made measurements of the BALs using these mean spectra for further use.We adopt the median MJD among each individual observing season as the effective MJD for each mean spectrum; this effectively gives us eight "epochs" to work with.We then measured the EW, mean depth, and v cent from the eight mean spectra and calculated ∆v cent and ∆v(CCF) between sequential pairs of mean spectra.These measurements are provided in Table 1.In addition, we provide measurements from our single HET spectrum, which was taken roughly two years after the last SDSS spectrum (see Appendix A), to characterize the activity of the BAL beyond the SDSS monitoring period.
We also calculate the measured acceleration between each pair of mean spectra using both ∆v cent and ∆v CCF .In some cases, depending on which method is used, we measure acceleration that is consistent with zero to within the 1σ uncertainties, indicating that on one-year timescales, velocity shifts are still often too small to detect at high confidence.However, we measure acceleration with a range of values from 0.04 cm s −2 to 3.90 cm s −2 using ∆v cent (hereafter we will refer to the acceleration calculated from ∆v cent as a cent ), and ranging from 0.05 cm s −2 to 2.05 cm s −2 when using ∆v CCF to calculate the acceleration (hereafter, a CCF ).We suspect that some of the differences between a cent and a CCF measurements are a result of the shape change that occurred in the final few years of monitoring; we thus see a more significant increase in a cent over that period, but suspect that not all of the measured velocity shift can be attributed to the actual acceleration of the gas.Previous studies (e.g., Grier et al. 2016, Byun et al. 2022, Xu et al. 2020) that report possible acceleration candidates have also measured a wide range of acceleration values (ranging from −0.25 to 1.5 cm s −2 ), supporting the idea that the acceleration can be variable in magnitude.
Examination of Figure 6 indicates a slow decrease in v cent over the monitoring period; however, the relationship between v cent and time does not appear to be linear.Figure 11 again shows v cent as a function of time, but this time we have fit a linear relation to the data and we see that the trend is not well fit (the reduced χ 2 of a linear fit to the data is 20.9) -there is an overall curvature to the trend, particularly during the latter half of the campaign.We instead fit a 2nd order polynomial to the data and we find this fit in much better agreement than the linear fit (the reduced χ 2 of the 2nd order polynomial is 10.2), though both are poor fits in general, as there is   ), and all acceleration measurements were all made between sequential epochs.All measurements were made using the mean spectra from each observing season, aside from the HET spectrum, which is the only available spectrum from that year.
significant short-timescale variability on top of the longterm trend.This would suggest that the magnitude of the acceleration is increasing over time.However, this interpretation is complicated by the fact that we also see a substantial change in the shape of the BAL during the latter half of the campaign (e.g., Figure 8).For the first half of the campaign, v cent and v CCF measure changes that are roughly on par with one another; however, in the latter half of the campaign, where the shape changes become more substantial, v CCF measures overall smaller shifts than v cent .This suggests that the non-linearity in the relationship between v cent and time is due to the changing shape of the BAL over this time rather than actual acceleration of the gas.
In addition to measuring a between sequential years, we also measure the acceleration between Years 1 and 8 to obtain the average acceleration over the entire monitoring period.We measure ∆v cent = −901 ± 35 km s −1 and ∆v CCF = −683 +88 −84 km s −1 , corresponding to a cent = −1.37±0.05cm s −2 and a CCF = −1.04+0.14 −0.13 cm s −2 .These values are the same order of magnitude as measurements from previous work, which examine BALs in quasars that have broadly similar properties.

Evolution of Trough A "Sub-Troughs"
Figures 7 and 8 clearly indicate a change in overall shape of Trough A; at the beginning of the monitoring period, we see two separate "sub-troughs" within the BAL, whereas the last year of monitoring shows a more broad, single-component trough.In an attempt to disentangle the behavior of these individual sub-troughs within Trough A from the observed velocity shift of the BAL, we model the BAL in each mean spectrum using two Gaussian profiles.We include only pixels that fall within the formal BAL limits in each mean spectrumi.e., our fits are restricted to the region where the normalized flux density lies below 0.9 (see Appendix B to see the actual model fits).We then examine the behavior of the two Gaussian components over time.Figure 12 shows the measured parameters of the Gaussian component representing each sub-trough (hereafter referred to as the red sub-trough and blue sub-trough) as a function of time.
We see that both sub-troughs are shifting blueward, though at slightly different rates: The red sub-trough begins with a line center of roughly 1470.5 Å and shifts blueward to a line center of 1467 Å (a shift of roughly 3.5 Å over the observing period), and the blue subtrough begins with a line center of roughly 1460.5 Å and shifts to 1457.5 Å, corresponding to a total shift of roughly 3 Å.The amplitude of the Gaussian fit of the blue sub-trough follows the overall EW trend during this period, starting off strong, followed by a significant decrease in amplitude in Years 2-6 and then a substantial increase in Years 7 and 8.The amplitude of the Gaussian representing the red sub-trough remains substantially weaker than that of the blue sub-trough throughout the campaign, and does not increase in amplitude during Years 7 and 8; by Year 7, the model is best fit by one major Gaussian (the blue sub-trough) that has significantly widened, with the second Gaussian contributing minimally.The measured widths of the blue Gaussian component increases over the monitoring period, whereas the width of the red component remains lower throughout.
The behavior of the two Gaussian components in our model fits suggests that both of the sub-troughs are shifting; however, the right sub-trough seems to be shifting blueward at a slightly higher rate than the left sub-trough: Between Years 1 and 8, the center of the red component shifted in velocity by 754 km s −1 (from −15405 km s −1 to −16159 km s −1 ), whereas the blue component moved by 625 km s −1 (from −17485 km s −1 to −18110 km s −1 ).This, combined with the increase in EW of the blue sub-trough of the BAL during the last two years of monitoring, causes the two sub-troughs to appear to merge into a single, broader trough while they shift.However, our HET spectrum taken two years after the SDSS observations indicates that the two components still remain distinct from one another, suggesting that they continue to behave as separate entities.

Measurements of C iv Troughs B and C
In addition to our measurements of Trough A, we also investigate the two narrower C iv absorption features present in RM 613, Troughs B and C (see Figure 4).Trough B has a velocity width of ∼ 1311 km s −1 , formally classifying it as a "mini-BAL" (e.g., Hall et al. 2013), and a centroid velocity of about −4800 km s −1 .Trough C is also formally classified as a mini-BAL, with a velocity width of ∼ 1380 km s −1 and a centroid velocity of about −1000 km s −1 .Grier et al. (2015) examined the behavior of both of these troughs during the first year of SDSS-RM monitoring and found that Trough B showed similar behavior in strength (EW and mean depth) as BAL Trough A, but Trough C did not show any significant variability in strength.The coordination between Troughs A and B implied that the observed variability is due to a change in ionization state of the gas.
Figure 13 shows the EW as a function of time for all three troughs for the eight-year monitoring period.We again find that the EW and depth of Trough B tracks that of Trough A; while there is somewhat more noise (likely due to the BAL being superimposed on the alsovariable C iv emission line), we see similar behavior in the strength of the BAL.However, we do not see any velocity shift in Trough B; the upper and lower boundaries of Trough B remain constant to within the uncertainties during the entire monitoring period.Trough C does not show any of the same variability trends in strength as the other two troughs -the EW shows very little variability

Possible Causes of Velocity Shifts
We below consider a few possible models to describe the behavior observed in RM 613 and evaluate the suitability of these models in producing this behavior.

Changes in ionization state
We first consider the possibility that changes in the ionization state of the gas or the covering factor have caused the observed velocity shift.While this possibility is difficult to rule out in previous cases where we had only 1-2 spectra to examine, the sheer number of observations we have of this source suggests that this is unlikely in the case of RM 613.It would take a remarkable coincidence of circumstances for the BAL to slowly and consistently vary in just such a way as to appear to show a monolithic increase in velocity over such a long period of time.For example, our modelling of Trough A as two Gaussians might have yielded two Gaussians at fixed velocities, with the red one weakening and the blue one strengthening to give the impression of acceleration.Instead, the best-fit model found increasingly negative velocities for both Gaussians.
In addition, inspection of the continuum light curve over the monitoring period (Figure 14; Shen et al. 2023) shows no unusual variability in the continuum flux of RM 613 during this time.However, the continuum flux is slightly elevated in Years 7 and 8 of the campaign, which also corresponds with the period in which we see the most dramatic changes in shape and the highest rates of acceleration.We examined the light curves of both the C iv and He ii 1640 Å emission lines, which sometimes trace the ultraviolet continuum more closely than the optical continuum (e.g., De Rosa et al. 2015), and we similarly see no unusual or extreme variability in either.The slope of the continuum increases in the negative direction as the quasar gets brighter (i.e., as the quasar continuum luminosity increases, the quasar gets bluer), again showing no unusual behavior throughout the monitoring period.

Geometric Effects
We next consider the possibility that geometric effects may produce the observed velocity shift; for example, if the outflow is launched from a rotating disk, its continued rotation may also cause an observed velocity shift (e.g., Hall et al. 2013, Grier et al. 2016) even if the speed of the outflow remains constant.However, assuming the rotating wind is silhouetted against the accretion disk, we would only see the wind that is appearing to decelerate rather than accelerating: the wind rotates into view with its maximum blueshift, gradually shifts to moving perpendicular to our line of sight (no shift), then becomes redshifted until it rotates out of view with its maximum redshift.Thus the rotation of a wind silhouetted against the accretion disk cannot produce our observed increase in velocity.
There is another possible geometric effect that may cause a velocity shift: a change in flow angle at con-stant speed.If the velocity change is due to a change in flow direction, for example due to a change in the azimuthal and/or polar angle of an outflow from an accretion disk where the outflow crosses our link of sight, we can determine the change in angle required to produce the observed velocity shift as follows.If θ LOS (t) is the angle between the flow and the line of sight and |v true | is the speed of the outflow, then: (3) With t = t 1 for the first observing season, we can solve for v true and obtain the line of sight velocity at any other time t when the outflow has the same v true but makes an angle θ LOS (t) with our line of sight: With our measured v LOS (t 1 ) = 16800 km s −1 , the minimum value of θ LOS (t 1 ) is 19 degrees (which requires θ LOS (t 8 ) = 0 degrees in the 8th observing season).If θ LOS (t 1 ) = 45 degrees, θ LOS (t 8 ) = 41.8 degrees.The maximum plausible θ LOS (t 1 ) = 73.7,as that yields v true = 60, 000 km s −1 , matching the fastest known UV outflows (Rogerson et al. 2015); in that case, θ LOS (t 8 ) = 72.8degrees.Thus, depending on the initial angle of the flow to the line of sight, a variation from 1 -19 degrees in that angle can explain the observed acceleration as a change in flow direction at fixed speed.This model cannot be ruled out, but does not explain all of the observed changes in v cent , which shows variations that are poorly described by a linear change (Figure 11).It is possible, however, that this effect, combined with the overall shape change that occurs in the BAL toward the latter half of the observations, may explain our observations.

Gas Dynamics
Another possible scenario is that the gas may be accelerating due to hydrodynamical effects such as overpressure, buoyancy, or entraining 1 .For example, Waters et al. (2021) discuss how bubbles of hot gas can form in a disk atmosphere and rise out of it, accelerating some cooler gas out of the atmosphere in the process.If the absorbing gas structure (traveling at ∼17000 km s −1 ) moves into a region of lower pressure, the structure will expand.Such expansion would only happen at the sound velocity (∼10 km s −1 ) or, at most, the turbulent velocity.A very large turbulent velocity of 900 km s −1 would be required to produce the observed velocity shift in our observations, which seems unlikely.This model would also predict that the acceleration is temporary, lasting only until the structure comes into pressure equilibrium again.Future observations of RM 613 will be revealing in this regard.
Alternatively, if the absorbing structure moves into a region where the surrounding gas is denser than the absorbing gas, the absorbing gas will accelerate away from the quasar due to a buoyancy force.However, denser surrounding gas should be lower-ionization than the absorbing gas and should itself produce absorption, so a pure buoyancy scenario seems implausible.
The visualizations of simulated disk wind outflows presented by Proga et al. (2012) illustrate how geometric and hydrodynamical effects can produce velocity shifts.In their Figure 1, between panels 3, 4 and 5 the blue edge of an absorption trough increases in outward velocity due to the motion of gas of different outward velocities across accretion-disk regions of different surface brightnesses, mimicking acceleration of a single flow structure.However, the red edge of the simulated trough does not exhibit a significant shift, which may indicate that a case like ours with a shift in velocity at both trough edges is more likely to be acceleration by radiation pressure, as discussed below.

Acceleration due to Radiation Pressure
We next consider that the observed acceleration could be due to an increase in speed of the outflow by incident ionizing radiation (e.g., Murray et al. 1995).This could be a case of our line of sight intersecting an outflow in a location where it is being accelerated to its terminal outflow velocity v ∞ , or a case of a previously coasting outflow being newly accelerated.
In § 4.1 of Grier et al. (2016), we provided an equation for the gas velocity and acceleration as a function of distance r from a black hole for a Murray et al. (1995) disk wind launched from a radius r = r L .The velocity is given by v(r) = v ∞ (1 − r L /r) 1.15 and the ratio of the acceleration to the velocity satisfies a(r)/v(r) = 1.15v(r)rL /r 2 .If we adopt an average acceleration of a = 1.21 ± 0.12 cm s −2 at a velocity of v = 1.6807 × 10 9 cm s −1 and assume a plausible value for r L , we can find values of r and v ∞ for which this model will match the observations.Following Grier et al. (2016) and assuming M BH ∼ 2 × 10 9 M ⊙ and r L = 3.6 × 10 17 cm = 0.12 pc = 610 R Sch , we find r = 2.73r L = 9.83 × 10 17 cm = 0.33 pc, and v ∞ = 28, 400 km s −1 .These numbers are only approximate, but they do indicate that to explain our observations with this model of gas accelerating to a fixed v ∞ , the gas needs to be located within a parsec of the central engine.For example, doubling r requires quadrupling r L (so that r = 1.365rL ) and more than doubling v ∞ to an implausible 76, 600 km s −1 , and cutting r L in half yields r = 3.86r L , or 0.23 pc, and v ∞ = 23, 700 km s −1 .
If we consider that the acceleration may have increased over time (Figure 11), then in this model either the gas is at a radius r < 1.65r L (since this model's acceleration is only increasing at such small radii) or the terminal velocity of the wind is changing with time.The terminal velocity is approximately (Laor & Brandt 2002).Here, v K (r L ) is the Keplerian circular speed at the launch radius and the effective Eddington ratio Υ eff = L inc /L Edd is the ratio of the quasar's luminosity incident on the wind gas (after accounting for any absorption interior to the wind) to its Eddington luminosity, and F M is the force multiplier acting on the wind gas.As discussed in § 2 of Murray et al. (1995), the force multiplier is larger where ∂vr ∂r is large (i.e., in regions of large acceleration) and where the ionization is low.
Naddaf et al. ( 2023) presents a different model in which radiation pressure on dust is the primary driver for the BAL outflow.From Fig. 4 of that paper, the launch radius for their model wind is about 10-20 times larger than for our model above, ranging from (6 − 1.2) × 10 3 R Sch .Dust could provide this wind with a larger F M , enabling it to reach the observed outflow velocities from a larger r L .This wind would have a similar v(r) profile to our model, with a length scale for reaching terminal velocity similar to the launch radius.Thus, in this dust-driven wind model, the BAL gas observed to be accelerating in this quasar would need to be located within 10-20 pc of the central engine.
As mentioned at the start of this section, it is also possible that the gas was previously coasting but has recently begun to accelerate.While we do not observe any dramatic/extreme variability in the optical continuum flux, small increases in the optical continuum flux may trace a larger increase in shorter-wavelength flux (such as that in the extreme ultraviolet) which could accelerate the outflow gas without over-ionizing it.A simple increase in the ionizing continuum level seems unlikely to always result in acceleration, since continuum variability is seen in BAL quasars without acceleration and deceleration being common (e.g., Filiz Ak et al. 2013;Grier et al. 2016).An unusually large or sustained ionizing continuum increase still might explain the acceleration, and we note that the velocity change in RM 613 was observed to be largest after MJD 58500, when the near-UV continuum flux and the EWs of Troughs A and B were all increasing.However, Wang et al. (2015) show that increases in the near-UV continuum are accompanied by increases in BAL trough EWs only ∼25% of the time (their Figure 7b), from which they conclude that most carbon atoms in the gas producing BAL troughs are in states more highly ionized than C iv.
In summary, it could be that the gas in the BAL trough observed in this object is located at a relatively small radius where it is still being accelerated by the quasar's ultraviolet radiation (≲ 1 pc in the Murray et al. 1995 model;≲ 10 − 20 pc in the Naddaf et al. 2023 model), even absent any increase in that radiation field.Alternatively, if the gas was previously coasting, we suggest that an unusual change in the ionizing SED consisting of an increase in ionizing flux and/or a softening of the ionizing spectrum could in principle explain both the acceleration and the increase in the C iv absorption EW, by producing more photons that scatter off and accelerate ions in the BAL gas but fewer photons that ionize C iv.We leave detailed study of these possibilities to future work.

On the Efficacy of Searches for BAL Acceleration
We are in the unique position of having more than a hundred epochs of spectra for our source; this makes the identification of acceleration fairly straightforward, as we can rule out velocity-dependent variability that mimics a velocity shift simply by visual inspection.However, in most cases, reports of velocity shifts/acceleration in BAL quasars have only a handful of epochs, frequently separated by several years.To date there has only been one systematic search for BAL acceleration in a large sample of quasars (Grier et al. 2016;hereafter G16), and this study included sources with only three epochs of observations.G16 provide a suggested procedure for identifying and quantifying the acceleration of a BAL based on the cross-correlation of spectra from two epochs and also requiring that the overall shape of the BAL remain constant between the observed epochs.This second criterion was included by necessity -with only 2-3 epochs, it is plausible that the cause of observed velocity shifts is changes in the ionization state of the gas or the transverse motion of gas into the line of sight, rather than the actual acceleration of the outflow.
G16 report a detection rate of acceleration in ∼ 1% − 10% of quasars in their study.However, this may not be the true incidence of acceleration in BALs.Due to observational constraints, such as the limited number of spectra with which most acceleration candidates are identified, G16 may not have been able to identify all cases in which the BAL actually accelerated.With ∼130 spectra providing us with a solid detection in RM 613, we can estimate, via simulations, the likelihood that we would have detected the acceleration in RM 613 if we only had 3 epochs of spectra over this time period.
As a way of determining this likelihood, we perform a series of simulations using our data from RM 613.We divide the spectra into three groups: Group 1 contains spectra from years 1-3, Group 2 spans years 4-6, and Group 3 spans years 7 and 8.We then randomly draw a spectrum from each set to produce a set of three epochs; division into the three groups is meant to enforce a time separation between spectra that is on par with that between spectra used by G16.We then carry out the tests defined by G16 to determine whether or not those three epochs are indicative of acceleration: We first measure the CCF between spectra 1 and 2, and between spectra 2 and 3, to determine if a velocity shift is present between these pairs.We then perform a χ2 test between these pairs both before and after having applied the measured velocity shift from the CCF.If the observed velocity shift is more than 3σ inconsistent with zero and the original spectra are a "bad" match, but the shifted spectra are a "good" match (as determined by the p-value measured in the χ 2 test), the pair of spectra is considered to have passed the criteria for acceleration.If both pairs (epochs 1-2 and 2-3) pass these tests, G16 determined that they are likely to show real acceleration signatures and are considered "acceleration candidates".
We perform this test with our RM 613 data 1000 times, randomly drawing epochs within each group and carrying out the G16 test for acceleration between the three randomly selected spectra.The results were quite striking: Out of 1000 draws, only 1 produced an "acceleration candidate" that passed the tests defined by G16.This low success rate has two potential causes: First, there is a much shorter timescale between the first pair of observations used in our study compared to those of G16 (the median ∆t rest between epochs 1 and 2 in our random draws was 322 days in the quasar rest frame, whereas the median ∆t rest between epochs 1 and 2 in the quasar sample examined by G16 is 1146 days). 2 The shorter timescales involved mean that the velocity shifts will be smaller, so we more rarely detect the velocity shift between epochs 1 and 2 at > 3σ significance in our simulations simply because there was not enough time for a large velocity shift to occur.Secondly, as noted above, the strength of the BAL in RM 613 is quite variable throughout the campaign and there is significant shape variability at the end of the monitoring period, causing epochs that are well separated in time to fail the required χ 2 tests.
Based on these simulations, we conclude that studies searching for acceleration that only include a handful of epochs, such as that of G16, will miss "real" acceleration much of the time, particularly if those epochs are not separated in time by at least a few years in the quasar rest frame.The actual rate of acceleration in BALs is likely much higher than that reported by G16 -the data, however, have thus far been generally insufficient for detection.High-cadence studies of BAL quasars over long timescales will be required to determine the actual rate of acceleration.5. SUMMARY Over the years, roughly a dozen cases of significant velocity shifts in BALs have been reported (see Section 1 and references therein).In those cases, only a handful of spectra were used, causing difficulties with disentangling line-profile variability with actual shifts in velocity.We have investigated the variability of a highvelocity C iv BAL trough with ∼130 spectroscopic observations spanning more than eight years, allowing us to explore variability and acceleration on both short and long timescales.Our major findings are as follows: 1. We have observed a long-term velocity shift in the C iv BAL of the quasar RM 613 over the monitoring period; the density of our observations indicates that this observed shift is not due solely to velocity-dependent variability in the BAL profile that might mimic an acceleration signature (Fig- ures 7 and 8).
2. We observe a very slow change in the shape of the BAL that becomes noticeable in the second half of the observing period.We do not see significant changes in shape on shorter timescales; however, over the span of our observations, the BAL changed from a two-pronged feature to a single smooth trough (Figures 7 and 8).
3. Combining all spectra taken within an individual observing season, we measure the acceleration of the trough between each observing season (Table 1), as well as over the entire monitoring period.Using the CCF, we measure the average velocity shift over the eight-year period to be −683 +89 −84 km s −1 , which corresponds to an average acceleration of a = −1.04+0.14 −0.13 cm s −2 .The magnitude of this acceleration varies throughout the monitoring period (Figures 9 and 11).
4. We consider a variety of possible causes of the observed velocity shift and determine that our observations are most likely to explained by geometric effects (Section 4.1.2) or the acceleration of the outflow by radiation pressure (Section 4.1.4).If our line of sight intersects the outflow in its acceleration region, that places the observed gas within 1-20 pc of the black hole, depending on the wind model adopted.
5. Our simulations suggest that studies such as G16 that search for velocity shifts in small samples are likely to miss actual acceleration much of the time due to short time baselines and the prevalence of variability in the shape and strength of BALs (Section 4.2).
Observations of BAL velocity shifts, the implied cause of which is the acceleration of outflow material, are difficult for a myriad of reasons.We used this unique dataset to investigate both the the short and longterm variability of this particular BAL simultaneously, which allowed us to determine that this BAL showed both line-profile variability and a distinct shift in velocity.Future studies focused on BAL acceleration will require long time-baseline, high-cadence observations to determine the frequency of such behavior as well as to disentangle signatures of line-profile variability from velocity shifts due to acceleration of the outflow.The SDSS-V BHM program will include observations on BAL quasars on several different timescales.The multi-epoch spectroscopy part of BHM aims to obtain a handful of observations of hundreds to thousands of quasars, though with a cadence similar to that from G16 we expect the yield of BAL-acceleration candidates from that sample to be low.However, we expect that there are roughly 200 BAL quasars in the sample that is currently being observed by the BHM-RM program, which will have a high observational cadence similar to that in this study -these data will be particularly instructive in searches for BAL acceleration.Additional surveys using multi-object spectrographs will aid further in the construction of a statistically significant sample of BAL acceleration candidates.B. MODEL FITS TO TROUGH A SUB-TROUGHS Section 3.4 presents a discussion of model fits to the two distinct sub-troughs of Trough A. We adopt a simple model, using a single Gaussian component for each sub-trough.Figure 16 shows these model fits, demonstrating the evolution of the two distinct sub-troughs over time.

Figure 1 .
Figure 1.The mean spectrum of RM 613, calculated from the 129 SDSS spectral epochs used in our study (see Section 2.2).The left panel shows the entire wavelength range covered by the SDSS spectra; the right panel is zoomed in on the C iv region of interest.The top panels show the flux density in units of 10 −17 erg s −1 cm −2 Å−1 and the bottom panels show the spectral uncertainties in the same units.

Figure 2 .
Figure 2. The spectrum of RM 613 taken on 30 May 2023 by the Hobby-Eberly Telescope (HET).The top panels show the flux density in units of 10 −15 erg s −1 cm −2 Å−1 and the bottom panels show the spectral uncertainties in the same units.The red shaded region indicates the C iv region of interest.The HET spectrum has been rebinned to the same wavelengths as the SDSS spectra using the SpectRes python package for display purposes.

Figure 3 .
Figure 3.An example SDSS spectrum of RM 613 and its corresponding continuum fit.The spectrum itself is shown in black, the continuum fit -including the C iv and Si iv emission-line fits -in red, and the green shaded regions highlight the four line-free regions used in the continuum fit.The vertical dashed green line shows the rest-frame wavelength of C iv.The right panel shows the same spectrum zoomed in to the C iv region of interest.The top panels show the flux density in units of 10 −17 erg s −1 cm −2 Å−1 and the bottom panels give the spectral uncertainties in the same units.

Figure 4 .
Figure 4. Mean spectra used in our study, after normalizing by the continuum and C iv emission-line fit.The three panels are zoomed into three regions of interest (the C iv , Si iv, and N v regions).Dotted vertical lines indicate the restframe line center of the three emission lines; the velocity given on the x-axes represents the velocity calculated against each particular emission line.The solid blue horizontal line indicates a normalized flux density of 1.0, and the red dotted horizontal line shows a flux density of 0.9, to aid the eye.The top subpanel for each region shows the mean normalized flux density, and the smaller subpanel shows the uncertainties on the mean normalized spectrum.

Figure 5 .
Figure 5.The redward (top panel) and blueward (middle panel) boundaries of the Trough A BAL as a function of the Modified Julian Date (MJD).Black points represent SDSS spectra, and the magenta triangles represent the measurements made from the HET spectrum.These boundaries are defined as the location at which the flux of the BAL reaches a normalized flux of 0.9 on either side of the BAL feature.Horizontal dashed lines indicate the median value throughout the monitoring, to aid the eye.The bottom panel shows the velocity width (vmax − vmin) of the BAL as a function of time using the determined limits.

Figure 6 .
Figure 6.Measurements of Trough A over time.The top panel shows the equivalent width (in units of Å), the second panel shows the mean depth of the BAL, and the third panel shows the absorbed-flux-weighted centroid velocity (vcent) of the BAL.Black points represent measurements made from each individual SDSS spectrum and the magenta triangles represent the measurements made from the HET spectrum.Red squares show the values measured from the mean spectra of each individual year (see Section 3.3).

Figure 7 .
Figure 7.A trailed spectrogram showing all of our continuum-normalized spectra over time.For visualization purposes, we did not include the emission-line fit in the normalization for this figure.Each row represents a different spectrum, and the color indicates the strength of the normalized flux density in each wavelength bin: absorption features appear as cool/blue vertical stripes, and emission lines appear as warm/red stripes.Vertical black dashed lines indicate the highest and lowest boundaries of the C iv BAL throughout the monitoring period.Horizontal gray lines indicate the breaks between observing seasons.Occasional narrow horizontal "stripes" of noise that appear are spectra that are noticeably lower in SNR than most of the other spectra.It is apparent that Trough A shifts to the left (to higher blueshifted velocities) throughout the monitoring.

Figure 8 .
Figure 8. Mean spectra from each observing year, focused in on the C iv BAL.The bottom (blue) line shows the mean from 2014, and each line moving upward represents the next year in sequence (2015-2021), with the top curve showing the mean in 2021.These mean spectra have been smoothed by five pixels using a boxcar smoothing algorithm, and each spectrum has been shifted upward in flux for visibility purposes.The horizontal blue dotted lines indicate a normalized flux level of 0.9 for each spectrum to aid the eye in identifying the upper and lower limits of the BAL.

Figure 9 .
Figure 9.The measured velocity shift as a function of time between all possible pairs of epochs, in the quasar rest frame.The left panel shows the velocity shift measured from the CCCD (∆vCCF) and the right panel shows the shift as measured by the difference in vcent between epochs (∆vcent).Gray lines indicate measurements made between all individual epoch pairs; black points show measurements made between the Year 1 mean spectrum and all additional mean spectra in subsequent years.

Figure 10 .
Figure 10.CCF, CCCD, CCPD analysis between the mean spectrum of 2014 (Year 1) and mean spectrum of 2021 (Year 8).The upper half of the figure shows the normalized spectra and their errors, along with the superposition of both on one figure.The vertical lines show the upper and lower boundaries of the BAL.The lower half of the figure shows the CCF, CCCD, and CCPD (see Section 3).

Figure 11 .
Figure 11.The velocity centroid vcent as a function of time.Black points represent measurements made from each individual SDSS spectrum and the magenta triangle represents the measurements made from the HET spectrum.Red squares show the values measured from the mean spectra of each individual year (see Section 3.3).A linear fit to the data is shown by the blue line, and a 2nd order polynomial fit is shown as a dashed green line.

Figure 12 .
Figure 12.The model fit parameters of the two Gaussian components of the Trough A model as a function of time.Blue symbols represent the higher-velocity, or bluer, Gaussian sub-trough component, and red represents the lowervelocity, or redder sub-trough component.The top panel shows the center wavelength of each component; the middle panel shows the amplitude of the Gaussian fits, and the bottom panel shows the width σ of the fit.The units of the amplitude are normalized flux density times Å.

Figure 13 .
Figure 13.The rest-frame EW of Troughs A, B, and C as a function of time (MJD).
beyond statistical scatter.Similar to Trough B, Trough C remains stable in velocity throughout the duration of our observations.
4. DISCUSSION While there have been other reports of monolithic velocity shifts in BALs over the years (see Section 1 for a list of such reports), this is the first time we have made a detection in a source with such dense time sampling to allow us to track the behavior of the BAL on shortto-long timescales.This allows us to evaluate several different mechanisms for producing the observed velocity shifts as well as to evaluate the efficacy of previous searches for BAL velocity shifts/acceleration.

Figure 14 .
Figure 14.Continuum flux, as measured from both photometric and spectroscopic observations of the quasar over time (light curves are from Shen et al. 2023).
RW and CJG acknowledge support from NSF grantsAST-2310211 and AST-2309930.Support  for this research was provided by the Office of the Vice Chancellor for Research and Graduate Education at the University of Wisconsin-Madison with funding from the Wisconsin Alumni Research Foundation.PBH acknowledges support from NSERC grant 2023-05068.WNB acknowledges support from NSF grant AST-2106990.YS acknowledges partial support from NSF grants AST-2009947 and AST-2108162.Y.H. was supported as an Eberly Research Fellow by the Eberly College of Science at the Pennsylvania State University.CR acknowledges support from the Fondecyt Iniciacion grant 11190831 and ANID BASAL project FB210003.RJA was supported by FONDECYT grant number 123171 and by the ANID BASAL project FB210003.MLM-A acknowledges financial support from Millenium Nucleus NCN19-058 and NCN2023 002 (TITANs).C.A.N. thanks support from CONAHCyT project Paradigmas y Controversias de la Ciencia 2022-320020 and DGAPA-UNAM grant IN111422.Funding for the Sloan Digital Sky Survey V has been provided by the Alfred P. Sloan Foundation, the Heising-Simons Foundation, the National Science Foundation, and the Participating Institutions.SDSS acknowledges support and resources from the Center for High-Performance Computing at the University of Utah.The SDSS web site is www.sdss.org.SDSS is managed by the Astrophysical Research Consortium for the Participating Institutions of the SDSS Collaboration, including the Carnegie Institution for Science, Chilean National Time Allocation Committee (CNTAC) ratified researchers, the Gotham Participation Group, Harvard University, Heidelberg University, The Johns Hopkins University, L'Ecole polytechnique fédérale de Lausanne (EPFL), Leibniz-Institut für Astrophysik Potsdam (AIP), Max-Planck-Institut für Astronomie (MPIA Heidelberg), Max-Planck-Institut für Extraterrestrische Physik (MPE), Nanjing University, National Astronomical Observatories of China (NAOC),

Figure 15 .
Figure 15.The HET spectrum (black line) after being normalized by the continuum.The measured limits of the BAL are shown as vertical red lines.The horizontal blue dotted line represents a normalized flux of 1.0, and the horizontal red line shows a normalized flux of 0.9.

Figure 16 .
Figure 16.The two-Gaussian model fits to the two sub-troughs within Trough A for the mean spectra of each individual year of monitoring.Each mean spectrum is shown in gray; red shows the redward Gaussian component of the model fit to the data and blue represents the blueward Gaussian component.

Table 1 .
BAL Measurements from Mean Spectra and HET Spectrum