Model-based Cognitive Neuroscience: Multifield Mechanistic Integration in Practice

Autonomist accounts of cognitive science suggest that cognitive model building and theory construction (can or should) proceed independently of findings in neuroscience. Common functionalist justifications of autonomy rely on there being relatively few constraints between neural structure and cognitive function (e.g., Weiskopf, 2011). In contrast, an integrative mechanistic perspective stresses the mutual constraining of structure and function (e.g., Piccinini & Craver, 2011; Povich, 2015). In this paper, I show how model-based cognitive neuroscience (MBCN) epitomizes the integrative mechanistic perspective and concentrates the most revolutionary elements of the cognitive neuroscience revolution (Boone & Piccinini, 2016). I also show how the prominent subset account of functional realization supports the integrative mechanistic perspective I take on MBCN and use it to clarify the intralevel and interlevel components of integration

Human Time-Frequency Acuity Beats the Fourier Uncertainty Principle

The time-frequency uncertainty principle states that the product of the temporal and frequency extents of a signal cannot be smaller than $1/(4\pi)$. We study human ability to simultaneously judge the frequency and the timing of a sound. Our subjects often exceeded the uncertainty limit, sometimes by more than tenfold, mostly through remarkable timing acuity. Our results establish a lower bound for the nonlinearity and complexity of the algorithms employed by our brains in parsing transient sounds, rule out simple "linear filter" models of early auditory processing, and highlight timing acuity as a central feature in auditory object processing.Comment: 4 pages, 2 figures; Accepted at PR

Phase transition in the Jarzynski estimator of free energy differences

The transition between a regime in which thermodynamic relations apply only to ensembles of small systems coupled to a large environment and a regime in which they can be used to characterize individual macroscopic systems is analyzed in terms of the change in behavior of the Jarzynski estimator of equilibrium free energy differences from nonequilibrium work measurements. Given a fixed number of measurements, the Jarzynski estimator is unbiased for sufficiently small systems. In these systems, the directionality of time is poorly defined and configurations that dominate the empirical average, but which are in fact typical of the reverse process, are sufficiently well sampled. As the system size increases the arrow of time becomes better defined. The dominant atypical fluctuations become rare and eventually cannot be sampled with the limited resources that are available. Asymptotically, only typical work values are measured. The Jarzynski estimator becomes maximally biased and approaches the exponential of minus the average work, which is the result that is expected from standard macroscopic thermodynamics. In the proper scaling limit, this regime change can be described in terms of a phase transition in variants of the random energy model (REM). This correspondence is explicitly demonstrated in several examples of physical interest: near-equilibrium processes in which the work distribution is Gaussian, the sudden compression of an ideal gas and adiabatic quasi-static volume changes in a dilute real gas.Comment: 29 pages, 5 figures, accepted for publication in Physical Review E (2012

Aperture synthesis for gravitational-wave data analysis: Deterministic Sources

Gravitational wave detectors now under construction are sensitive to the phase of the incident gravitational waves. Correspondingly, the signals from the different detectors can be combined, in the analysis, to simulate a single detector of greater amplitude and directional sensitivity: in short, aperture synthesis. Here we consider the problem of aperture synthesis in the special case of a search for a source whose waveform is known in detail: \textit{e.g.,} compact binary inspiral. We derive the likelihood function for joint output of several detectors as a function of the parameters that describe the signal and find the optimal matched filter for the detection of the known signal. Our results allow for the presence of noise that is correlated between the several detectors. While their derivation is specialized to the case of Gaussian noise we show that the results obtained are, in fact, appropriate in a well-defined, information-theoretic sense even when the noise is non-Gaussian in character. The analysis described here stands in distinction to ``coincidence analyses'', wherein the data from each of several detectors is studied in isolation to produce a list of candidate events, which are then compared to search for coincidences that might indicate common origin in a gravitational wave signal. We compare these two analyses --- optimal filtering and coincidence --- in a series of numerical examples, showing that the optimal filtering analysis always yields a greater detection efficiency for given false alarm rate, even when the detector noise is strongly non-Gaussian.Comment: 39 pages, 4 figures, submitted to Phys. Rev.

Three-dimensional sound propagation models using the parabolic-equation approximation and the split-step Fourier method

Author Posting. © IMACS, 2012. This article is posted here by permission of World Scientific Publishing for personal use, not for redistribution. The definitive version was published in Journal of Computational Acoustics 21 (2013): 1250018, doi:10.1142/S0218396X1250018X.The split-step Fourier method is used in three-dimensional parabolic-equation (PE) models to compute underwater sound propagation in one direction (i.e. forward). The method is implemented in both Cartesian (x, y, z) and cylindrical (r, θ, z) coordinate systems, with forward defined as along x and radial coordinate r, respectively. The Cartesian model has uniform resolution throughout the domain, and has errors that increase with azimuthal angle from the x axis. The cylindrical model has consistent validity in each azimuthal direction, but a fixed cylindrical grid of radials cannot produce uniform resolution. Two different methods to achieve more uniform resolution in the cylindrical PE model are presented. One of the methods is to increase the grid points in azimuth, as a function of r, according to nonaliased sampling theory. The other is to make use of a fixed arc-length grid. In addition, a point-source starter is derived for the three-dimensional Cartesian PE model. Results from idealized seamount and slope calculations are shown to compare and verify the performance of the three methods.This work was sponsored by the Office of Naval Research under the grants N00014-10-1-0040 and N00014-11-1-0701

Private quantum decoupling and secure disposal of information

Given a bipartite system, correlations between its subsystems can be understood as information that each one carries about the other. In order to give a model-independent description of secure information disposal, we propose the paradigm of private quantum decoupling, corresponding to locally reducing correlations in a given bipartite quantum state without transferring them to the environment. In this framework, the concept of private local randomness naturally arises as a resource, and total correlations get divided into eliminable and ineliminable ones. We prove upper and lower bounds on the amount of ineliminable correlations present in an arbitrary bipartite state, and show that, in tripartite pure states, ineliminable correlations satisfy a monogamy constraint, making apparent their quantum nature. A relation with entanglement theory is provided by showing that ineliminable correlations constitute an entanglement parameter. In the limit of infinitely many copies of the initial state provided, we compute the regularized ineliminable correlations to be measured by the coherent information, which is thus equipped with a new operational interpretation. In particular, our results imply that two subsystems can be privately decoupled if their joint state is separable.Comment: Child of 0807.3594 v2: minor changes v3: presentation improved, one figure added v4: extended version with a lot of discussions and examples v5: published versio

Notes on the integration of numerical relativity waveforms

A primary goal of numerical relativity is to provide estimates of the wave strain, $h$, from strong gravitational wave sources, to be used in detector templates. The simulations, however, typically measure waves in terms of the Weyl curvature component, $\psi_4$. Assuming Bondi gauge, transforming to the strain $h$ reduces to integration of $\psi_4$ twice in time. Integrations performed in either the time or frequency domain, however, lead to secular non-linear drifts in the resulting strain $h$. These non-linear drifts are not explained by the two unknown integration constants which can at most result in linear drifts. We identify a number of fundamental difficulties which can arise from integrating finite length, discretely sampled and noisy data streams. These issues are an artifact of post-processing data. They are independent of the characteristics of the original simulation, such as gauge or numerical method used. We suggest, however, a simple procedure for integrating numerical waveforms in the frequency domain, which is effective at strongly reducing spurious secular non-linear drifts in the resulting strain.Comment: 23 pages, 10 figures, matches final published versio

Floquet-Markov description of the parametrically driven, dissipative harmonic quantum oscillator

Using the parametrically driven harmonic oscillator as a working example, we study two different Markovian approaches to the quantum dynamics of a periodically driven system with dissipation. In the simpler approach, the driving enters the master equation for the reduced density operator only in the Hamiltonian term. An improved master equation is achieved by treating the entire driven system within the Floquet formalism and coupling it to the reservoir as a whole. The different ensuing evolution equations are compared in various representations, particularly as Fokker-Planck equations for the Wigner function. On all levels of approximation, these evolution equations retain the periodicity of the driving, so that their solutions have Floquet form and represent eigenfunctions of a non-unitary propagator over a single period of the driving. We discuss asymptotic states in the long-time limit as well as the conservative and the high-temperature limits. Numerical results obtained within the different Markov approximations are compared with the exact path-integral solution. The application of the improved Floquet-Markov scheme becomes increasingly important when considering stronger driving and lower temperatures.Comment: 29 pages, 7 figure

Observing binary inspiral in gravitational radiation: One interferometer

We investigate the sensitivity of individual LIGO/VIRGO-like interferometers and the precision with which they can determine the characteristics of an inspiralling binary system. Since the two interferometers of the LIGO detector share nearly the same orientation, their joint sensitivity is similar to that of a single, more sensitive interferometer. We express our results for a single interferometer of both initial and advanced LIGO design, and also for the LIGO detector in the limit that its two interferometers share exactly the same orientation. We approximate the evolution of a binary system as driven exclusively by leading order quadrupole gravitational radiation. To assess the sensitivity, we calculate the rate at which sources are expected to be observed, the range to which they are observable, and the precision with which characteristic quantities describing the observed binary system can be determined. Assuming a conservative rate density for coalescing neutron star binary systems we expect that the advanced LIGO detector will observe approximately 69~yr${}^{-1}$ with an amplitude SNR greater than 8. Of these, approximately 7~yr${}^{-1}$ will be from binaries at distances greater than 950~Mpc. We explore the sensitivity of these results to a tunable parameter in the interferometer design (the recycling frequency). The optimum choice of the parameter is dependent on the goal of the observations, e.g., maximizing the rate of detections or maximizing the precision of measurement. We determine the optimum parameter values for these two cases.Comment: 40 pages (plus 7 figures), LaTeX/REVTEX3.0, NU-GR-

Positive Quantum Brownian Evolution

Using the independent oscillator model with an arbitrary system potential, we derive a quantum Brownian equation assuming a correlated total initial state. Although not of Lindblad form, the equation preserves positivity of the density operator on a restricted set of initial states