We introduce SF-LIFE, a large-scale simulated movement dataset designed to accelerate research in transportation, mobility, and machine learning. The dataset contains 3,024,000,000,000 location records capturing complete, noise-free, multi-modality trajectories of 500,000 simulated agents observed at a 1Hz frequency navigating the San Francisco Bay Area network over a 70-day period. The data captures (1) needs-driven daily agendas of individual agents generated by an agent-based simulation of human patterns of life and (2) detailed kinematic trajectories moving agents across the OpenStreetMap representation of San Francisco using data from 40+ transit agencies across 9 counties. SF-LIFE provides unprecedented scale and detail as trajectories are based on real transit infrastructure using San Francisco General Transit Feed Specification (GTFS) data, having agent movements across multiple modalities, including bus, rail, bike, automobile, and walking. For this high-fidelity simulated representation of San Francisco, we provide (1) the full trajectory data annotated with transportation mode labels, (2) reduced-size versions of the trajectory data with reduced temporal frequency, (3) a
Urban Air mobility has gained momentum with recent advancements in the electric vertical take-off and landing (eVTOL) vehicles, offering faster point-to-point air taxi services that could help relieve traffic congestion in chronically overburdened cities. The research assesses the feasibility and systems-of-systems level adoption potential of UAM operations in the San Francisco Bay Area by comparing passenger departure, waiting, travel, and arrival times across key regional nodes, including San Francisco, Oakland, San Jose, and Palo Alto airports, with conventional ground transportation. A multi-agent simulation was developed in MATLAB to evaluate the fleet operations and to model demand arrival using a Poisson process under stochastic passenger flows and turnaround constraints. Results indicate that utilizing UAM during peak demand could reduce total travel times up to eighty percent across the region. The findings of this paper highlight the critical operational factors for fleet schedule optimization. Especially how the fleet size, passengers' request volumes, and turnaround time directly influence waiting time, operating cost, and overall user acceptance.
It is well known that Charles Hermite kept an intense correspondence with many of the word's leading mathematicians of his time. This paper focuses on Hermite's letters to Francisco Gomes Teixeira, a Portuguese mathematician, who exchanged letters with Hermite for more than twenty years.
This study presents the development of two new sedimentary velocity models for the San Francisco Bay Area (SFBA) to improve the near-surface representation of shear-wave velocity ($V_S$) for large-scale, broadband numerical simulations, with the ultimate goal of enhancing the representation of the sedimentary layers in the Bay Area community velocity model. The first velocity model is stationary and is based solely on $V_{S30}$; the second velocity model is spatially varying and has location-specific adjustments. They were developed using a dataset of 200 measured $V_S$ profiles. Both models were formulated within a hierarchical Bayesian framework, using a parameterization that ensures robust scaling. The spatially varying model includes a slope adjustment term modeled as a Gaussian process to capture site-specific effects based on location. Residual analysis shows that both models are unbiased for $V_S$ values up to 1000 m/sec. Along-depth variability models were also developed using within-profile residuals. The proposed models show higher $V_S$ in the San Jose area and Livermore Valley compared to the USGS Bay Area community velocity model by a factor of two or more in some case
Congestion pricing, while adopted by many cities to alleviate traffic congestion, raises concerns about widening socioeconomic disparities due to its disproportionate impact on low-income travelers. We address this concern by proposing a new class of congestion pricing schemes that not only minimize total travel time, but also incorporate an equity objective, reducing disparities in the relative change in travel costs across populations with different incomes, following the implementation of tolls. Our analysis builds on a congestion game model with heterogeneous traveler populations. We present four pricing schemes that account for practical considerations, such as the ability to charge differentiated tolls to various traveler populations and the option to toll all or only a subset of edges in the network. We evaluate our pricing schemes in the calibrated freeway network of the San Francisco Bay Area. We demonstrate that the proposed congestion pricing schemes improve both the total travel time and the equity objective compared to the current pricing scheme. Our results further show that pricing schemes charging differentiated prices to traveler populations with varying value-of-t
The University of California San Francisco Brain Metastases Stereotactic Radiosurgery (UCSF-BMSR) dataset is a public, clinical, multimodal brain MRI dataset consisting of 560 brain MRIs from 412 patients with expert annotations of 5136 brain metastases. Data consists of registered and skull stripped T1 post-contrast, T1 pre-contrast, FLAIR and subtraction (T1 pre-contrast - T1 post-contrast) images and voxelwise segmentations of enhancing brain metastases in NifTI format. The dataset also includes patient demographics, surgical status and primary cancer types. The UCSF-BSMR has been made publicly available in the hopes that researchers will use these data to push the boundaries of AI applications for brain metastases. The dataset is freely available for non-commercial use at https://imagingdatasets.ucsf.edu/dataset/1
In 2020, California required San Francisco to consider equity in access to resources such as housing, transportation, and emergency services as it re-opened its economy post-pandemic. Using a public dataset maintained by the San Francisco Fire Department of every call received related to emergency response from January 2003 to April 2021, we calculated the response times and distances to the closest of 48 fire stations and 14 local emergency rooms. We used logistic regression to determine the probability of meeting the averages of response time, distance from a fire station, and distance to an emergency room based on the median income bracket of a ZIP code based on IRS statement of income data. ZIP codes in the lowest bracket ($25,000-$50,000 annually) consistently had the lowest probability of meeting average response metrics. This was most notable for distances to emergency rooms, where calls from ZIP codes in the lowest income bracket had an 11.5% chance of being within the city's average distance (1 mile) of an emergency room, while the next lowest probability (for the income bracket of $100,000-$200,000 annually) was 75.9%. As San Francisco considers equity as a part of Califo
Climate change in California may require construction of a barrier separating the Pacific Ocean from San Francisco Bay and the Sacramento River-San Joaquin River Delta simply because Southern California is remarkably dependent on freshwater exported from the Delta. We offer a new kind of salt barrier, a macroproject built of impermeable textile materials stretched across the Golden Gate beneath the famous bridge. We anticipate it might eventually substitute for a recently proposed San Francisco In-Stream Tidal Power Plant harnessing a 1.7 m tide at the Bay entrance if future climate conditions Statewide is conducive. First-glance physics underpin our macroproject.
San Francisco Crime Classification is an online competition administered by Kaggle Inc. The competition aims at predicting the future crimes based on a given set of geographical and time-based features. In this paper, I achieved a an accuracy that ranks at top %18, as of May 19th, 2016. I will explore the data, and explain in details the tools I used to achieve that result.
Urban transportation and land use models have used theory and statistical modeling methods to develop model systems that are useful in planning applications. Machine learning methods have been considered too 'black box', lacking interpretability, and their use has been limited within the land use and transportation modeling literature. We present a use case in which predictive accuracy is of primary importance, and compare the use of random forest regression to multiple regression using ordinary least squares, to predict rents per square foot in the San Francisco Bay Area using a large volume of rental listings scraped from the Craigslist website. We find that we are able to obtain useful predictions from both models using almost exclusively local accessibility variables, though the predictive accuracy of the random forest model is substantially higher.
We have observed fallout from the recent Fukushima Dai-ichi reactor accident in samples of rainwater collected in the San Francisco Bay area. Gamma ray spectra measured from these samples show clear evidence of fission products - 131,132I, 132Te, and 134,137Cs. The activity levels we have measured for these isotopes are very low and pose no health risk to the public.
This paper presents a geometric time-domain method for identifying three-phase load equivalents from instantaneous voltage and current measurements at the point of common coupling. Measured waveforms are interpreted as trajectories in Euclidean signal spaces, and load-equivalent parameters are recovered from the geometry of those trajectories. The method extends a previously published single-phase geometric identification formulation to three- and four-wire systems and places special emphasis on the three-wire case, where no neutral voltage is measured and the terminal data must satisfy coupled Kirchhoff constraints. The main advance over the earlier analytical formulation is a sampled-data implementation based on local time windows, normalized matrix equations, harmonic-projection derivative and primitive coordinates, explicit geometric identifiability tests, passivity constraints, and energy/Kirchhoff residuals. The method does not force a model when the measured trajectory lacks enough information; instead, it reports low-rank or ill-conditioned windows as low-confidence evidence. Numerical simulations with clean data, measurement noise, window-length sweeps, and sensor delay sh
Recent work on causal abstraction, in particular graphical approaches focusing on causal structure between clusters of variables, aims to summarize a high-dimensional causal structure in terms of a low-dimensional one. Existing methods for learning such summaries from data assume that both the high- and low-dimensional structures are acyclic, which is helpful for causal effect identification and reasoning but excludes many high-dimensional models and thus limits applicability. We show that in the linear non-Gaussian (LiNG) setting, the high-dimensional acyclicity assumption can be relaxed while still allowing recovery of a low-dimensional causal directed acyclic graph (DAG). We further connect identifiability of this low-dimensional DAG to existing results: LiNG models with cycles are observationally identifiable only up to an equivalence class whose members differ by reversals of directed cycles; our low-dimensional DAG, which is invariant across all members of a given equivalence class, thus forms a natural representative of the class. While existing approaches for learning this observational equivalence class over high-dimensional variables have exponential time complexity, our
Henze and Malikiosis (2017) have shown that the Lonely Runner Conjecture (LRC) can be restated as a convex-geometric question on the so-called LR zonotopes, lattice zonotopes with one more generator than their dimension. This relation naturally suggests a more generel statement, the "shifted" LRC, the zonotopal version of which concerns a classical parameter, the covering radius. Theorems A and B in Malikiosis-Schymura-Santos (2025) use the zonotopal restatements of both the original and the shifted LRC to prove a linearly-exponential bound on the size of the (integer) speeds for which the conjectures need to be checked in order to establish them for each fixed number of runners; in the shifted version their statement and proof rely on a certain assumption on two-dimensional rational vector configurations, the so-called "Lonely Vector Property". In this paper we do two things: We push the analogies between the two versions of LRC and their zonotopal counterparts, in particular highlighting that the proofs of Theorems A and B in Malikiosis-Schymura-Santos are more transparent, and the statments more general, if regarded in terms of two quite general classes of lattice zonotopes: the
This paper develops a branch-level energy-localization framework for three-phase loads. The instantaneous terminal power of an admissible lumped equivalent is decomposed uniquely as Joule dissipation plus magnetic and electric stored-energy rates, branch by branch. Three formal results are established: a Branch-Level Localization Theorem (uniqueness given an admissible topology); a Topology-Indeterminacy Theorem (multiple admissible topologies reproduce identical terminal data with distinct localizations); and a Generalized Energetic Duality Theorem that organizes classical electrical dualities (Norton-Thevenin, series--parallel, L vs C, R vs G) as restrictions to Linear Time Invariant (LTI) sinusoidal regimes of a single time-domain principle in which constant-parameter equivalence is replaced by time-varying parameters. The framework is exercised on six test cases including the de Leon--Cohen open-phase paradox, switched-resistive loads, three-wire delta-versus-wye-virtual indeterminacy, fluctuating-phase loads, and a four-wire nonlinear load with hysteretic, linear, and switched branches. The framework is positioned as complementary to IEEE Std. 1459, CPC, instantaneous p-q, and
We use Diffusing Wave Spectroscopy (DWS) to perform the first direct space- and time-resolved measurement of the dissipation rate~$ε$ at the boundary of a turbulent flow. We have shown in a previous publication that this technique provides maps of the dissipation rate of Newtonian fluids~\cite{Francisco}. Here, we apply the technique at the boundary of a turbulent flow generated in a square box by an impeller stirring the fluids. Although the measurement is made on a small region near the boundary, we show that the dissipation remains proportional to the injected power and follows the turbulent scaling $ε\propto \mathrm{Re}^3$, with Re being the Reynolds number ranging from $1.5 \times 10^4$ to $6 \times 10^5$. With this flow, there is no need for logarithmic corrections to reproduce the dissipation near the flat boundary. In addition, our setup allows us to measure the spatio-temporal fluctuations of the dissipation near the boundary. These fluctuations are quite large (the relative fluctuations are about 50\%) and are well described by a log-normal distribution, as expected for the dissipation rate in the bulk of homogeneous and isotropic turbulence (HIT) but Power Density Spectr
In this work, we study the realisation of unitarity-based cutting rules for primordial cosmological correlators computed within the Schwinger-Keldysh path integral formalism. While cutting rules have been previously derived for wavefunction coefficients, here we examine them directly at the level of cosmological observables expressed diagrammatically. The resulting rules closely resemble those familiar from flat-space scattering amplitudes, but with an additional subtlety: in order to express the discontinuity of a correlator as the product of lower-order correlators, one must introduce a specific combinations of diagrams which do not appear in the computation of observables themselves. We explicitly verify these rules for several classes of correlators, both at tree level and with loop corrections, arising from theories involving different types of interactions.
We show that the polyhedron defined as the convex hull of the lattice points above the hyperbola $\left\{xy = n\right\}$ has between $Ω(n^{1/3})$ and $O(n^{1/3} \log n)$ vertices. The same bounds apply to any hyperbola with rational slopes except that instead of $n$ we have $n/Δ$ in the lower bound and by $\max\left\{Δ, n/Δ\right\}$ in the upper bound, where $Δ\in \mathbb{Z}_{>0}$ is the discriminant. We also give an algorithm that enumerates the vertices of these convex hulls in logarithmic time per vertex. One motivation for such an algorithm is the deterministic factorization of integers.
We show that the shifted Lonely Runner Conjecture (sLRC) holds for 5 runners. We also determine that there are exactly 3 primitive tight instances of the conjecture, only two of which are tight for the non-shifted conjecture (LRC). Our proof is computational, relying on a rephrasing of the sLRC in terms of covering radii of certain zonotopes (Henze and Malikiosis, 2017), and on an upper bound for the (integer) velocities to be checked (Malikiosis, Santos and Schymura, 2024+). As a tool for the proof, we devise an algorithm for bounding the covering radius of rational lattice polytopes, based on constructing dyadic fundamental domains.
This paper presents experimental validation of a time-domain load parameter determination method for single-phase circuits. The verification is performed in a state-of-the-art smart grid laboratory equipped with power hardware and real-time emulators. The proposed method enables the identification of circuit parameters using only instantaneous voltage and current measurements at the point of common coupling. The experimental setup includes a range of test cases covering linear and non-sinusoidal single-phase conditions. Voltage and current waveforms are acquired, preprocessed, and used to calculate the relevant circuit parameters. The experimental results demonstrate a high degree of accuracy and robustness, with minimal percentage errors across all test cases. The identified parameters show excellent agreement with the theoretical expectations, confirming the validity and applicability of the proposed method to identify the load of single-phase systems. This validation highlights the potential of the method for improved monitoring, control, and protection of smart grids, paving the way for future extensions to three-phase systems and real-time implementations.