\label{sec:numerical-results}
We began by simulating the full two parameter space by using the flatPERM algorithm \cite{Prellberg_2004} which is a flat histogram version of the Pruned and Enriched Rosenbluth Method (PERM) developed in \cite{Grassberger_1997}. For the PERM algorithm, at each iteration a polymer configurations is generated kinetically (which is to say that each growth step is selected at random from all possible growth steps) along with a weight factor to correct the sample bias. At each growth step, configuration with very high weight relative to other configurations of the same size are enriched (duplicated) while configuration with low weight or that cannot be grown any further are pruned (discarded). Despite introducing a correlation between each iteration, this simple mechanism greatly improves the algorithm efficiency. A single iteration is then concluded when all configurations have been pruned and the total number of samples generated during each iteration depends on the specifities of the problem at hand and on the details of the enriching/pruning strategy. FlatPERM extends this method by cleverly choosing the enrichment and pruning steps to generate for each polymer size \(n\) a quasi-flat historgram in some choosen micro-canonical quantities \(\mathbf{k}=(k_1,k_2,\dotsc,k_{\ell})\) and producing an estimate \(W_{n,\mathbf{k}}\) of the total weight of the walks of length \(n\) at fixed values of \(\mathbf{k}\). From the total weight one can access physical quantities over a broad range of temperatures through a simple weighted average, e.g. \[\begin{aligned} {\langle}\mathcal O {\rangle}_n(\tau) = \frac{\sum_{\mathbf{k}} \mathcal O_{n,\mathbf{k}}\, \left(\prod_j \tau_j^{k_j}\right) \, W_{n,\mathbf{k}}}{\sum_\mathbf{k} \left(\prod_j \tau_j^{k_j}\right) \, W_{n,\mathbf{k}}}.\end{aligned}\] The quantities \(k_j\) may be any subset of the physical parameters of the model. To study the full two parameter phase space we set \((k_1, k_2) = (m_t, m_c)\) and \((\tau_1, \tau_2) = (\tau, \omega)\).
We have first simulated the model using the full two-parameter flatPERM algorithm up to length \(n = 256\), running \(4.4 \cdot 10^6\) iterations, and collecting \(2.3 \cdot 10^{11}\) samples at the maximum length. To obtain a landscape of possible phase transitions, we plot the largest eigenvalue of the matrix of second derivatives of the free energy with respect to \(\tau\) and \(\omega\) (measuring the strength of the fluctuations and covariance in \(m_t\) and \(m_c\)) at length \(n=256\) on Figure \ref{fig:eigenvalue_plot}.