Hi Simone,
Based on my understanding, it appears that you are seeking to apply orthogonal rotation to new data and obtain the corresponding factor scores using the "factoran()" method.
To achieve this, you can utilize the "Rotate" parameter and set its value to "orthomax" to specify the orthogonal rotation method. Additionally, you can utilize the "Scores" parameter to compute the factor scores. There are two possible values you can specify for the "Scores" parameter, namely "wls" (weighted least-squares) and "regression".
For more detailed information and guidance on this topic, I recommend referring to the following page: Estimate Factor Scores from New Data.
I hope this information assists you in successfully obtaining the factor scores for your matrix.
Thanks,
Sachin.