In this thesis, we extend the numerical S-matrix bootstrap program to 1+1d theories with a boundary, where we bootstrap the 1-to-1 reflection matrix (R-matrix). We review the constraints that a physical R-matrix must obey, namely unitarity, analyticiy and crossing symmetry. We then carve out the allowed space of 2d R-matrices with the O(N) nonlinear sigma model and the periodic Yang Baxter solution in the bulk. We find a variety of integrable R-matrices along the boundary of the allowed space both with and without free parameters. The integrable models without a free parameter appear at vertices of the allowed space, while those with a free parameter occupy the whole boundary. We also introduce the extended analyticity constraint where we increase the domain of analyticity beyond the physical region. In some cases, the allowed space of R-matrices shrinks drastically and we observe new vertices which correspond to integrable theories. We also find a new integrable R-matrix through our numerics, which we later obtained by solving the boundary Yang--Baxter equation. Finally, we derive the dual to the extended analyticity problem and find that the formalism allows for R-matrices which do not saturate unitarity to lie on the boundary of the allowed region.

<p>Entanglement entropy (EE) is a quantum information theoretic measure that quantifies the correlations between a region and its surroundings. We study this quantity in the following two setups : </p> <ul> <li>We look at the dynamics of a free minimally coupled, massless scalar field in a deSitter expansion, where the expansion stops after some time (i.e. we quench the expansion) and transitions to flat spacetime. We study the evolution of entanglement entropy (EE) and the Rényi entropy of a spatial region during the expansion and, more interestingly, after the expansion stops, calculating its time evolution numerically. The EE increases during the expansion but the growth is much more rapid after the expansion ends, finally saturating at late times, with saturation values obeying a volume law. The final state of the subregion is a partially thermalized state, reminiscent of a Gibbs ensemble. We comment on application of our results to the question of when and how cosmological perturbations decohere.</li> <li>We study the EE in a theory that is holographically dual to a BTZ black hole geometry in the presence of a scalar field, using the Ryu-Takayangi (RT) formula. Gaberdiel and Gopakumar had conjectured that the theory of N free fermions in 1+1 dimensions, for large N, is dual to a higher spin gravity theory with two scalar fields in 2+1 dimensions. So, we choose our boundary theory to be the theory of N free Dirac fermions with a uniformly winding mass, m e^iqx, in two spacetime dimensions (which describes for instance a superconducting current in an N-channel wire). However, to O(m²), thermodynamic quantities can be computed using Einstein gravity. We aim to check if the same holds true for entanglement entropy (EE). Doing calculations on both sides of the duality, we find that general relativity does indeed correctly account for EE of single intervals to O(m²).</li> </ul>