I am wondering if there exist CPTP maps $T$ such that the purity of a quantum state $\rho$ can increase, i.e.
$$ \text{tr} ( T ( \rho )^2 ) \geq \text{tr} ( \rho ^2). $$
If so, what are the conditions on $T$ and/or $\rho$ for this to be possible?
I am wondering if there exist CPTP maps $T$ such that the purity of a quantum state $\rho$ can increase, i.e.
$$ \text{tr} ( T ( \rho )^2 ) \geq \text{tr} ( \rho ^2). $$
If so, what are the conditions on $T$ and/or $\rho$ for this to be possible?
Yes, some quantum channels can increase purity. For example the preparation channel $$ T(X) = \mathrm{Tr}[X] |\psi\rangle \langle \psi| $$ that can be thought of as throwing away your system and preparing $|\psi\rangle$ will always produce a pure state. Therefore $\mathrm{Tr}[T(\rho)^2] = 1$ for any $\rho$.
A channel $T$ is called unital if $T(I) = I$ (i.e. preserves the identity matrix). If $T$ is a unital channel then it is possible to show that $\mathrm{Tr}[T(X)^2] \leq \mathrm{Tr}[X^2]$. Thus any channel that increases purity must be non-unital. It turns out non-unitality is necessary and sufficient see Theorem II.4 in Contractivity of positive and trace preserving maps under Lp norms.
Here's the simplest argument I could come up with for the statements Rammus made about unital channels. We actually don't need complete positivity, just positivity is enough.
Suppose that $T$ is a positive, unital, and trace-preserving map and let $\rho$ be a density operator input to $T$. Consider a spectral decomposition of $\rho$. $$ \rho = \sum_{k = 1}^{n} p_k \vert\psi_k\rangle\langle\psi_k\vert $$
Define $\sigma_k = T(\vert\psi_k\rangle\langle\psi_k\vert)$ for $k = 1,\ldots,n$. We can use our assumptions on $T$ as follows.
We can now obtain the desired inequality by Cauchy-Schwarz.
\begin{align*} \operatorname{Tr}(T(\rho)^2) & =\sum_{j,k=1}^n p_j p_k \operatorname{Tr}(\sigma_j\sigma_k)\\ & \leq\sqrt{\sum_{j,k=1}^n p_j^2 \operatorname{Tr}(\sigma_j\sigma_k)} \sqrt{\sum_{j,k=1}^n p_k^2 \operatorname{Tr}(\sigma_j\sigma_k)}\\ & = \sum_{k=1}^n p_k^2\\ & = \operatorname{Tr}(\rho^2) \end{align*} For the inequality, we're making use of the fact that $\operatorname{Tr}(\sigma_j\sigma_k)$ is always a nonnegative real number, which follows from the fact that $\sigma_1,\ldots\sigma_n$ are positive semidefinite, and for the equality that follows we're using both the equality $\sigma_1+\cdots+\sigma_n = \mathbb{I}$ and the fact that $\sigma_1,\ldots,\sigma_n$ all have unit trace.
On the other hand, if $T$ isn't unital, then the inequality above will be violated for the completely mixed state $\rho = \mathbb{I}/n$. This is because $\sigma = T(\rho)$ is a density operator that isn't completely mixed, so $$ 0 < \operatorname{Tr}((\sigma - \rho)^2) = \operatorname{Tr}(\sigma^2) - \frac{1}{n} = \operatorname{Tr}(T(\rho)^2) - \operatorname{Tr}(\rho^2). $$ In short, the completely mixed state is uniquely qualified as the most impure a state can be, so you necessarily get more pure if you move away from it.
Another common example of a CPTP channel which can increase purity (excluded from that link of canonical channels) is amplitude damping, with Kraus operators (parameterised by probability $\lambda$) $$ K_1 = \begin{pmatrix} 1 & \\ & \sqrt{1-\lambda} \end{pmatrix}, \;\;\;\;\;\;\;\; K_2 = \begin{pmatrix} & \sqrt{\lambda} \\ & \end{pmatrix}, $$ which act on a general one-qubit state $\rho$ as $$ \rho \rightarrow K_1 \rho K_1^\dagger + K_2 \rho K_2^{\dagger}. $$ This channel represents a simple dissipative process toward a minimum observable state, above assumed as $|0\rangle$. It is easy to see that as $\lambda$ approaches $1$ and ergo the probability of the dissipative transition $|1\rangle \rightarrow |0\rangle$ becomes certain, the maps become $$ K_1 \rightarrow |0\rangle\langle 0|, \;\;\;\; K_2 \rightarrow |1\rangle\langle 0|. $$ The maps take both $|0\rangle$ and $|1\rangle$ states and output a completely pure $|0\rangle$ state, even if the input $\rho$ is maximally mixed.